F/8  17/8 


/  AD-A097  187  ADAPTIVE  OPTICS  ASSOCIATES  INC  CAM8RID6E  HA 

ADVANCED  WAVEFRONT  SENSOR  CONCEPTS. <U> 

JAN  81  J  K  BOWKER.  J  FEINLEIB.  L  SCHNUTZ  F30602-80-C-0063 

RADC-TR-80-368  NL 


UNCLASSIFIED 


AD  A  097  1  8  7 


<f  LEVEL* 

KABC-TR-tO-369  p 

Hwi  Tdwld  ispsit 
Jammy  19S1 


ADVANCED  WAVEFRONT 
SENSOR  CONCEPTS 

Adaptive  Optic*  Associates,  Inc. 


J.  K.  Bowkor  L  Schmuti 

J.  Foinloib  M.  Sham 

S.  J.  Tubbs 


1  SWtOVtO  FOS  rustic  IftUSI;  wstmsution  UNLIMITED  I 


DT1C 

SELECTE1% 

APR  1  1981 1  ■ 


ROME  AIR  DEVELOPMENT  CENTER 
Air  Foret  Systems  Command 
Griffiss  Air  Fore#  Ban,  Now  York  13441 

81  4  01  023 


This  report  her  boon  reviewed  by  the  SADC  Public  Affaire  Office  (PA)  and 
la  releaaable  to  the  National  Technical  Information  Service  (HTIS) .  At  NTIS 
it  will  be  releaaable  to  the  general  public,  including  foreign  nation*. 

RADC-TR-80-368  has  been  reviewed  and  is  approved  for  publication. 


APPROVED: 


Patrick!  j.martone,  i/Lt, 


Project  Engineer 


USAF 


APPROVED: 

FRANK  J.  REHM 
Technical  Director 
Surveillance  Division 


FOR  THE 

J^jcm  p.  huss 

Acting  Chief,  Plans  Office 


SUBJECT  TO  EXPORT  CONTROL  LANS 


If  your  address'  has  changed  or  if  you  wish  to  be  reaoved  fron  the. SADC 
■ailing  Hat,  or  if  the  addressee  is  no  longer  eaployed  by  your  organisation, 
please  notify  RADC.(OCSE)  Griff lee  AFB  NT  13441.  This  will  assist  .us  in 
maintaining  a  current  sailing  list. 


Do  .not  return  this  copy.  Retain  or  destroy 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  FACE  fWhtn  Omt»  Erw.r.ai 


M ^REPORT  documentation  page 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

A— B*BO«tLm*IE*R  _ _  |]  OOVT  ACCESSION  no 

RADC*TR-80-368  1 A 

S  RECIPIENT'S  CATALOG  NUMBER 

a.  nni  <Lr#  BRiwwy — 

Q 

^ADVANCED  JJAVEFR0NT  SENSOR  CONCEPTS  g  ' 

S  type  OF  REPORT  a.  PERIOD. GOVE  AE0 

;  FinalTechnical  Rep**t# 
Jan*Bl-  Sep  8 j6± 

A.  RERFORMIHO  ORO.  REPORT  HUMBER - - 

.  . 

; J./Felnleib 
i L. /Schmutz 


j  Shao 

s/ J.yTubbs 


/^/F3jtap-80-C-^jJ63 


WTTation  name  ANO  AOORESS 

Adaptive  Optics  Associates,  Inc. 


12  Prentiss  Street  4o  4  4  M  /) 

Cambridge  MA  03140 _ 


}J^~ 


10.  pnosram  element,  fnoject.  task 
AREA  A  WORK  unit  NUMBERS 


*^702F  — n  . 

45061424 


n.  controlling  optics  namc  and  aoork&s 

Rome  Air  Development  Center  (OCSE) 


n 


»a.  rcrort  oat* 

JairfjjrTfcBl 


Griffiss  AFB  NY  13441 


IX  NUMBER  OF  RACE* 

183 


11 


/ 


I*  MONlVoRINO  AOEnCv  name  a  AOORESS/lf  FKIamhI  Inm  Conmlllni  Olll cm) 

Same 


is.  security  class.  iKn 
UNCLASSIFIED 


V*.  01 Rl BU flON  STAtImEnT  /•/  Ihlt  Kmpatu 

Approved  for  public  release;  distribution  unlimited. 


'Sa.  DECLASSIFICATION/  DOWNGRADING 
SCHEDULE 

N/A 


IT.  DISTRIBUTION  STATEMENT  Im I  !*•  aAamci  mtlmrwd  in  BlacA  20,  II  dllltfnl  I w  R. port) 

Same 


I*.  SUPPLEMENTARY  notes 

RADC  Project  Engineer: 


Patrick  J.  Martone,  lLt,  USAF  (OCSE) 


K|V  WOROt  ON  NVNIO  0I#0  it  IIKOIMry  iRNfffr  6r  bloc*  wuwQOf) 


Adaptive  Optics 
Wavefront  Sensors 
Hartmann  Sensor 
Atmospheric  Turbulence 


10.  ABSTRACT  (CohiInm  oft  rovoroo  »*4o  it  mcoooon  iRNflfp  *y  Mac*  nwRr) 

The  objective  of  this  effort  was  to  Investigate  advances  in  wavefront 
sensor  technology  for  use  in  future  adaptive  optics  systems.  The 
improvements  sought  over  presently  developed  adaptive  optics  systems  are 
preliminary  in  the  areas  of  componenet  efficiency  and  sub-system 
complexity,  since  this  would  have  potential  for  major  operational 
simplicity  and  reduced  cost. 

( Con  t ' J ) 


UNCLASSIFIED 


Ml MjO 


SECURITY  CL  AAfIFiCATlON  OF  This  PACE  ,RRfh  Dm*  Zmam « 


00  •  J  AN  TJ  1473  EOlTlON  of  I  NOV  *t  It  OBSOLETE 


UNCLASSIFIED 


SlCUWlTY  CLASSIFICATION  or  THU  EAgCQWHw  0««  Bntmd) 


^The  concept  under  study_in  this  effort  is  the  Large  Scale  Integrated 
Imaging  Irradiance  (LSl^)  sensor.  This  is  based  on  the  Integrated 
Imaging  Irradiance  (1~)  sensor  (U.S;  Patent  No.  4,141,652).  The  I  sensorl 
is  based  on  the  Hartmann  concept,  but  does  not  suffer  from  the  optical 
alignment  problems  inherent  to  previously  developed  Hartmann  sensors. 
Another  of  the  major  advantages  of  the  LSi^-1  sensor  over  other  Hartmann^ 
sensors  (including  the  previously  developed  I^-sensor)  is  that  the  LSf—i 
sensor  utilizes  the  full  aperture  of  the  system  in  dividing  the  input 
photon  flux  into  the  respective  quadrants  which  results  in  better 
diffraction  -  limited  resolution.  The  LSIliii  sensor  is  based  on  an  optical 
system  that  is  independent  of  the  number  of  sub-apertures  selected.-,. 

The  major  tasks  of  this  study  included  the  conceptual  design  of  a  small 
version  (approximately  20  channels),  a  large  scale  (approximately  200 
channels)  wavefront  sensor,  and  the  evaluation  of  the  expected  performance 
complexity,  and  cost  of  these  prototype  models.  ^During  the  course  of  the 
study,  major  design  modifications  were  introduced (which  appeared  to 
greatly  improve  the  performance  of  the  LSI^  concept  of  the  1^  sensor,  and 
finally  showed  a  significant  performance  improvement  over  previously 
developed  wavefront  sensors  including  the  shearing  interferometer.  At  the 
same  time,  these  design  changes  appear  to  have  greatly  reduced  the 
complexity  of  the  electronic  processing  required,  simplified  the  optical 
system  to  a  few  critical  elements  and  allowed  large  diode  array  detectors 
such  as  photomultipliers  to  be  utilized  without  a  significant  performance 
loss. 


The  result  of  these  design  advancements  is  a  new  wavefront  sensor  (LSf^1) 
which  could  provide  cost  and  performance  advances  for  future  adaptive 
optics  systems. 


Accession  For 
GRAJcI 


NTIS 
DTIC  TAB 
Unannounced 
Justification 


□ 

□ 


Rv  .  _ 

Distribution/ 
Availability  Coded 
.Avail  and/or 
Dlst  j  Special 


A 


unclassified 


;L  istinc  ATiO*  O*  -'•AC.K'Whan  Dmtm 


SUMMARY 


This  final  report  describes  the  results  of  a  nine  month 
effort  by  Adaptive  Optics  Associates,  Inc.  (AOA)  on  Rome  Air 
Development  Center  (RADC)  Contract  F30602-80-C-0063  for  an 
"Advanced  Wavefront  Sensor  Concepts"  (AWSC)  study.  The 
objective  of  this  study  was  to  investigate  advances  in  wavefront 
sensor  technology  for  use  in  future  adaptive  optics  systems,  and, 
in  particular,  to  advance  the  technology  over  the  Compensated 
Imaging  System  (CIS)1  presently  being  developed  for  open  space 
surveillance  application.  The  improvements  sought  are  primarily 
in  the  area  of  component  efficiency  and  subsystem  complexity, 
since  this  would  have  potential  for  major  operational  simplicity 
and  reduced  cost.  The  technical  performance  projected  for  CIS 
is  near  optimum  and  only  small  improvements,  if  any,  were  contem¬ 
plated  at  the  start  of  this  program. 

To  meet  the  objectives  of  this  program,  AOA  proposed  to 

analyze  and  evaluate  the  design  and  performance  of  a  new  wave- 

o 

front  sensor  concept  based  on  a  proprietary  sensor,  the  Inte- 

3 

grated  Imaging  Irradiance  (I  )  sensor,  developed  and  tested 

3  4 

under  other  Air  Force  programs  ’  .  The  concept  under  study 

3  3 

here  will  be  referred  to  as  the  LSI  ,  or  Large  Scale  I  , sensor 
throughout  this  report. 

The  major  tasks  of  this  study  included  the  conceptual  design 
of  a  small  version  (approximately  20  channels)  and  a  large  scale 
(approximately  200  channels)  wavefront  sensor,  and  the  evaluation 
of  the  expected  performance,  complexity  and  cost  of  these  proto- 

3 

type  models.  The  study  was  undertaken  using  the  baseline  LSI 
proposed.  During  the  course  of  the  study,  major  design  modifica¬ 
tions  were  introduced  which  appeared  to  greatly  improve  the  per- 

3  3 

formance  of  the  LSI  concept  over  the  I  sensor,  and  finally 
showed  a  significant  performance  improvement  over  that  projected 


1 


for  the  CIS.  At  the  same  time,  these  design  changes  greatly 
reduced  the  complexity  of  the  electronic  processing  required, 
simplified  the  optical  system  to  a  few  critical  elements  and 
allowed  large  diode  array  detectors,  instead  of  individual 
quantum  detectors  such  as  photomultipliers,  to  be  utilized  with¬ 
out  a  significant  performance  loss. 

As  a  result  of  these  design  advancements,  the  AWSC  investi¬ 
gation  met  all  program  objectives,  and  this  final  report 
describes  in  detail  the  design  and  analysis  of  a  new  wavefront 
sensor  which  could  provide  cost  and  performance  advances  for 
future  CIS  systems. 


REFERENCES 


1  References  will  be  made  available  to  qualified  military  and 
government  agencies  on  request  from  RADC  (OCSE)  Gnffiss 
AFB  NY  13441. 


2.  U.  S.  Patent  No.  4,141,652. 


3.  AFWL-TK- 78-65. 


4. 


Reference 

government 


will  be  made  available  to  qualified  military  and 
agencies  on  request  from  RADC  (OCSE)  Griff iss 


AFB  NY  13441. 


5.  Patent  pending. 


2 


TABLE  OF  CONTENTS 


Page 


1.0  INTRODUCTION  AND  BACKGROUND  .  9 

1.1  Compensated  Imaging  Wavefront 

Sensing  Problem  .  9 

3 

1.2  I  Sensor  Concept .  11 

3 

1.3  LSI  Sensor  Concept .  17 

1.4  Design  Goals .  18 

1.5  Parallel/Serial  System  Design  ...  19 

2.0  THE  LSI3  SENSOR  CONCEPT .  25 

3 

2.1  LSI  Optical  Configuration .  25 

2.2  Signal  Processing  .  29 

2.3  Detector  Selection  .  31 

2.4  Signal  Processing  Electronics  ...  34 

3.0  SIGNAL  PROCESSING  AND  PERFORMANCE  ....  36 

3.1  Introduction .  36 

3.2  LSI3,  MSI3  and  SI .  38 

3.3  Performance  Evaluation  .  42 

3.4  Optimal  Correlator  &  Object 

Size  Determination . 45 

4.0  OPTICAL  SYSTEM  DESIGN  .  54 

4.1  Reflective  Systems  .  65 

5.0  SIGNAL  PROCESSING  ELECTRONICS  .  69 

5.1  Processor  Architecture  .  . . 69 

5.2  Front  End  Processor  . . 69 

5.3  Programmable  Algorithm  Processor  .  .  77 

5.4  System  Supervisor  Processor  ....  88 

6.0  DETECTOR  STUDY . 92 

6.1  Direct  Mode  Systems . 99 


3 


Section 


Page 


6.2  Geometries  and  Processing 

Requirements  .  106 

6.3  Signal  Processing  .  119 

6.4  Input  to  Arrays:  Techniques  ....  142 

6.5  Electron  Bombardment  Mode .  143 

6.6  Photon  Coupled  Systems  .  150 

6.7  Electron  Coupled  Mode .  159 

6.8  Recommended  Detector  Design  .  164 

7.0  THE  21 -ELEMENT  BREADBOARD  SYSTEM  .  177 

8.0  CONCLUSION .  182 

LIST  OF  ILLUSTRATIONS 

Figure  Page 

3 

1.  Simplified  I  sensor  optics  .......  13 

3 

2.  Response  of  I  sensor  to  two  different 

image  plane  intensity  distributions  .  .  16 

3.  Simplified  optical  layout  .  20,21 

4.  Detector  arrays  used  as  pupil  dividers  .  .  23 

3 

5.  Hartmann,  SI  and  LSI  sensors .  26 

6.  Measurement  of  wavefront  tilt .  28 

7.  Optical  schematic  .  37 

8.  Phase  variance  vs.  range .  47 

9.  Relative  SNR  vs.  nutation  radius .  48 

10.  Relative  SNR  vs.  nutation  radius .  49 


4 


Figure  Page 

11.  Tilt  noise  as  function  of 

number  of  detected  photons  .  50 

12.  Phase  variance  vs.  range  for  several 

wavefront  sensors  ....  .  51 

13.  Double  pupil  image  optics  .  56 

14.  Single  pupil  image  optics  .  57 

15.  Single  pupil  image  optics  with 

field  lens . 59 

16.  Design . 61 

17.  Design,  reflective  system  .  66 

3 

18.  LSI  signal  processing . 70 

19.  Front  end  processor . 72 

3 

20.  LSI  sensor  analog  electronics  ....  73 

21.  Solid  state  line  scanner 

preamplifier  timing . 74 

22.  SCS  algorithm  computation . 78 

23.  SCS  pipeline . 79 

24.  AM  2903  bit  slice  circuit . 82 

25.  PAP  element  . 83 

26.  PAP  microcode  for  SCS  algorithm  ....  84 

27.  PAP  microcode  for  SCS  algorithm 

and  derotation . 86 

28.  PAP  architecture . 87 

29.  System  Supervisor  Processor  (SSP)  ...  89 

30.  Direct  detection  coupling  .  94 

31.  Electron  bombardment  coupled  .  95 

32.  Photon  coupled  devices.  ...  .  96 

33.  Anode  arrays . 97 

34.  PMT  arrays  for  breadboard . 98 

i 


5 


Page 


Figure 

35.  Solar  spectrum  and  detector 

sensitivities  .  100 

36.  Buried  channel  CCD  structure  .  109 

37.  Surface  channel  CCD  structure  ....  110 

38.  Array  array  subset  readout  technique  .  Ill 

39.  Structure  of  CID . 115 

40.  Schematic  of  Reticon  self-scanning 

diode  array  . 116 

41.  Fairchild  linear  image  sensor  ....  118 

42.  Rectangular  object  .  132 

43.  SSD  weighting  function . 133 

44.  Photosite  structure  .  147 

45.  Phosphor  spectral  output  and  fiber 

optics  plate  transmittance  ....  152 

46.  Phosphor  decay  curves  .  157 

47.  Preferred  EB  mode  detector, 

internal  optics  .  165 

48.  Aberrations  of  spherical  and 

ellipsoid  mirrors  .  170 

49.  EB  mode  detector,  external  optics  ,  .  166 

50.  EB  mode  detector  with  clam¬ 

shell  optics  .  167,168 

LIST  OF  TABLES 

Table  Page 

1.  Types  of  masks  and  modulation 

functions .  39 

2.  Demodulation  functions  .  41 

3.  Size  dependence  of  various 

algorithms .  52 


6 


Table 


Page 


4.  Integration  result,  total 

photo-electron  count  .  101 

5.  Bi-alkali  cathode  counts  and  rates  .  .  .  105 

6.  Surface/buried  channel  comparisons  .  .  .  107 

7.  Tradeoffs  between  arrays  .  107 

8.  Optimum  weighting  functions  .  135 

9.  CCD  weights  . 135 

10.  Reticon  nl  elements . 140 

11.  Recommended  detectors  .  173 

12.  Relationship  of  breadboard  to  AWS  ....  181 


7 


EVALUATION 


The  value  and  significance  of  this  effort  includes  many  design  advance¬ 
ments  incorporated  in  a  new  wavefront  sensor  (Large  Scale  Integrated 

3 

Imaging  Irradiance  LSI  )  which  could  provide  cost  and  performance  advances 
for  future  adaptive  optics  systems. 

The  results  of  this  effort  are  applicable  to  Electro-Optics  Advanced 
Systems  Analysis. 

Applications  of  this  Advanced  Wavefront  Sensor  include  future 
Compensated  Imaging  Systems  for  Space  Defense  applications,  and  laser  beam 
control  systems  such  as  the  Stragegic  Laser  Communications  Uplink  program. 


PATRICK  J.  MARTONE,  1/Lt,  USAF 
Project  Engineer 
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1.0  INTRODUCTION  AND  BACKGROUND 


1.1  Compensated  Imaging  Wavefront  Sensing  Problem 

The  development  of  real-time  wavefront  sensing  has 
produced  three  techniques  of  measurement:  interferometric, 
successfully  embodied  in  the  shearing  interferometer  which 
measures  phase  shifts  in  the  pupil;  geometric,  embodied  in  the 
Hartmann  test,  which  measures  wavefront  tilts;  and  spread  func¬ 
tion  analysis,  which  obtains  the  pupil  phase  distortion  by 
Fourier  transformations.  The  first  of  these  has  been  completely 
developed  by  Itek  under  contract  to  RADC  which  produced  the 
Real  Time  Atmospheric  Correction  (RTAC)^,  a  21-channel  system. 
Later,  after  a  design  competition  with  Perkin  Elmer,  which 
espoused  a  Hartmann  system  under  another  CIS  program,  Itek 
developed  a  250-channel  digital  sensor  for  CIS.  Because  of  these 
contracts,  other  techniques  have  not  been  developed  extensively. 

The  design  contest  between  Itek  Corporation  and  Perkin 
Elmer  which  knocked  the  Hartmann  test  out  was  won  by  Itek,  in 
part  because  the  Hartmann  test  was  optically  and  electronically 
too  sensitive  to  misalignment  and  drift  —  so  sensitive,  in  fact, 

that  a  system  could  not  be  demonstrated.  The  appreciation  of 

3 

these  difficulties  was  the  basis  for  the  invention  of  the  I 

2 

sensor,  patented  and  developed  by  AOA  . 

The  CIS  wavefront  sensing  problem  is  severe  because  the 
targets  are  dim,  large,  complex  and  moving  rapidly.  In  conse¬ 
quence,  the  sensor  measurement  must  be  insensitive  to  object 
shape,  must  be  able  to  operate  with  as  few  as  one  photon/cycle 
on  the  average,  and  must  provide  a  correction  bandwidth  as  great 
as  1  kHz. 


1.  Reference  will  made  available  to  qualified  military  and 
government  agencies  on  request  from  RADC  (OCSE)  Griffiss 
AFB  NY  13441. 

2.  U.S.  Patent  No.  4,141,652. 
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Target  size  and  shape  effect  have  been  crucial  in 

determining  whether  a  sensor  will  work  or  not,  and  these  have 

generally  made  adaptive  function  necessary  in  the  design  of 

sensors.  For  example,  if  the  target  subtends  more  resolution 

elements  over  its  area  than  the  system  has  subapertures,  then 

the  target  will  be  resolved  by  the  subapertures.  The  shearing 

interferometer  signal  visibility  will  drop  and  n  phase  error 

occurs  due  to  interference  effects.  The  classic  Hartmann  test 

will  lose  sensitivity,  often  dropping  to  zero,  if  the  target 

has  multiple  bright  spots.  The  shearing  interferometer  (SI) 

can  adapt  to  these  difficulties  by  decreasing  the  shear  and 

reducing  sensitivity.  The  Hartmann  test  cannot  adapt  and  is 

vulnerable  to  object-induced  failure. 

3 

The  I  sensor,  although  based  on  the  Hartmann  test, 
uses  an  AC  algorithm  and,  in  essence,  uses  edge  information  as 
well  as  the  broad  radiant  distribution.  Thus,  response  becomes 
less  dependent  on  object  size  and  signal  strength.  The  adaptive 
features  of  the  proposed  sensor  are  used  to  optimize  sensitivity 
and  to  prevent  performance  anomalies  such  as  those  embodied  by 
the  dead  zone  of  the  dumbbell  response  of  the  Hartmann  or  by  the 

visibility  nulls  with  2tt  ambiguity  and  n  phase  flips  of  the  SI. 

3  3 

The  LSI  sensor  is  a  major  advance  of  the  I  sensor 

for  applications  requiring  a  large  number  of  subaperture  measure¬ 
ments.  It  embodies  ail  of  the  optical  efficiencies  of  the  Hart- 

3 

mann  sensor  and  the  I  insensitivity  to  target  shaoes. 

The  major  advantages  are  as  follows. 

-  Fundamental  performance  improvement  over  the  shearing 
interferometer  by  up  to  a  factor  of  two  for  simple  targets  of  all  sizes 
and  ranges. 

-  An  optical  system  that  is  independent  of  the  number 

3 

of  subapertures  selected,  unlike  the  I  sensor. 


10 
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!■ 
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-  An  optical  system  that  is  independent  of  the 

3 

detector  configuration,  unlike  the  I  sensor. 

-  An  optical  system  design  that  can  be  used  with  pulsed 
or  CW  actively  illuminated  targets  as  well  as  passive  illuminated 
targets. 

-  A  versatile  optimum  processing  algorithm  that  can 
be  adapted  to  maximize  performance  for  unusual  target  shapes. 

-  An  electronic  processor  architecture  that  requires 
only  four  pipeline  serial  processors  to  service  256  subapertures,  as 
opposed  to  1024  parallel  processors. 

-  A  selection  of  three  different  intensified  serial 
detector  array  configurations  which  would  maintain  an  optical  efficiency 
exceeding  that  projected  for  the  current  CIS  and  that  are  compatible 
with  the  pipeline  processor. 

-  A  system  design  that  incorporates  all  the  performance 
advantages  in  a  small  lightweight  optical  package  using  few  components 
and  simple  to  align. 

-  A  system  design  with  low  component  count,  compact 
electronics  processing  package  that  promises  high  reliability. 

These  features  will  be  discussed  in  detail  in  the  following 

sections. 

3 

1.2  I  Sensor  Concept 

q  12  3 

The  1°  sensor  *  ’  is  a  wavefront  sensing  concept  based  on 


q 

1.  Julius  Feinleib  et  al .  ,  "The  I  Sensor:  A  New  Method  for 
Real-Time  Wavefront  Mensuration,"  Proc.  SPIE  Tech.  Symp.  East, 
Washington,  D.C.(1979). 

2.  J.  K.  Bowker  et  al . ,  "I3  Sensor  Wavefront  Performance," 

SPIE  Vol  190  TT973). 

3 

3.  l.  E.  Schmutz  et  al . ,  "Experimental  Performance  of  the  I 
Wavefront  Sensor  for  Closed-Loop  Adaptive  Optics,"  Proc.  SPIE 
Tech.  Symp.  East,  Washington  D.C.  (1980). 

i  1 


the  Hartmann  technique.  In  the  Hartmann  test,  the  shape  of  an 
optical  wavefront  is  estimated  by  measuring  the  mean  tilts 
at  each  of  a  set  of  small  subapertures  in  the  pupil  aperture, 
and  then  fitting  a  smooth  curve  to  the  tilt  samples.  Dynamic 
Hartmann  sensors,  designed  to  measure  wavefront  changes  at  rates 
comparable  to  atmospheric  turbulence  variations,  have  been 
attempted  but  have  not  succeeded  due  to  severe  problems  of 
optical  alignment  stability,  and  to  a  fairly  high  sensitivity  to 
image  shape  in  cases  where  the  source  object  is  resolved  by  the 
subaperture . 

3 

Ref.  1  (p.ll)  is  an  introduction  to  the  I  sensor  concept 
and  describes  some  of  its  performance  advantages.  Figure  1  is 
a  schematic  of  the  concept  and  shows  the  changes  from  a  basic 
Hartmann  technique.  The  two  major  innovations  which  make  the 

3 

I  sensor  a  practical,  high  sensitivity  wavefront  measuring 
technique  are: 

(1)  the  addition  to  the  input  wave  of  a  local  plane 
wave  which  acts  as  an  internal  optical  tilt  reference,  and 

(2)  a  conical  scan,  or  nutation  modulation  of  the 
combined  signal  and  reference  beam. 

1.2.1  Optical  reference.  The  use  of  a  local  reference 
provides  enormous  flexibility  in  the  design  of  a  practical 
measuring  instrument.  Once  combined,  the  signal  and  reference 
beams  are  equivalently  affected  by  internal  optical  misalignments, 
which  might  be  due  to  thermally  or  vibrationally  induced  compo¬ 
nent  movement,  or  to  finite  specification  tolerances  of  the 
optical  elements.  In  the  subaperture  image  plane,  tilts  in  the 
input  wavefront  are  given  by  the  difference  between  the  reference 

and  signal  focal  spot  positions.  The  differential  character  of 
3 

the  I  measurement  is  useful  with  any  kind  of  position  sensing 
detector,  be  it  a  quad  cell,  lateral  cell,  CCD  or  diode  array, 
video  imager,  etc. 


12 


SIGNAL 

WAVEFRONT 


Figure  1 
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Simplified  I  sensor  optics 


1.2.2  Nutation.  When  coupled  with  the  quadrant 
detector  configuration,  nutation  modulation  permits  the  use 
of  some  extremely  powerful  and  efficient  detection  and  processing 
techniques.  Several  performance  advantages  can  be  obtained, 
which  include  optimal  SNR  performance,  AC  carrier  generation,  and 
the  ability  to  discriminate  between  signal  and  background  radia¬ 
tion  in  the  complex  scenes. 

Nutation  modulation  impresses  an  AC  carrier  on 
the  optical  input  power  by  utilising  spatial  variations  on  the 
image  itself.  It  differs  from  both  the  AC  shearing  interfero¬ 
meter  and  the  chopped  Hartmann  type  of  system  in  that  the  input 
photon  flux  is  never  interrupted  by  a  chopped  or  rotating  grating. 
In  those  systems,  uniform  illumination  from  the  object  background 
is  modulated  along  with  the  object  signal  power  which  is  limited 
by  field  stops  of  the  detector. 

However,  when  nutation  is  used,  uniform  background 
power  is  not  modulated  at  all,  since  the  relative  power  falling 
on  each  quadrant  of  the  detector  is  unchanged  during  the  course  of 
a  nutation  cycle.  In  addition,  locally  generated  radiation, 
either  from  scattering  or  thermal  sources,  does  not  appear  at  the 
AC  carrier  frequency.  This  has  been  found  experimentally  to  have 
very  significant  advantages,  particularly  when  using  infrared 
( IH)  sensitive  detectors. 

Another  more  subtle  consequence  of  nutation 

3 

modulation  is  the  ability  of  the  I  sensor  to  operate  with  equal 
facility  using  objects  of  either  positive  or  negative  contrast. 
That  is,  the  object  may  be  a  bright  spot  against  a  dark  back¬ 
ground  or  a  dark  spot  in  a  bright  background.  In  both  cases, 
correct  signs  for  the  wavefront  tilt  response  are  retained  (a 
crucial  result  for  closed-loop  operation). 


H 


r 
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1.2.3  Image  insensitivity.  The  unique  features  of 
“ 

nutation  modulation  extend  also  to  the  way  in  which  the  I 
sensor  responds  to  complex  shapes  resolved  by  the  subaperture, 
which  is  a  problem  quite  relevant  to  imaging  systems.  One 
manifestation  of  the  problem  occurs  in  the  change  of  the  transfer 
function  of  the  sensor  for  different  image  shapes.  In  a  closed- 
loop  system  this  amounts  to  changes  in  loop  gain  for  different 
input  scenes,  complicating  the  problem  of  obtaining  unconditional 
loop  stability.  For  the  shearing  interferometer,  image  periodi¬ 
cities  related  to  the  grating  period  cause  significant  signal 
variations,  while  in  simple  quad  Hartmann  systems,  a  "dumbbell" 
image,  or  pair  of  bright  spots  located  in  separate  quadrants, 

results  in  near-zero  signal  regardless  of  spot  intensity. 

3 

The  I  sensor  responds  to  the  slope  of  the  image 

edges,  rather  than  to  the  centroid  of  the  image  power.  In 

a  properly  designed  instrument,  the  edge  slope  of  the  subaperture 

image  is  dominated  by  the  point  spread  function  of  the  subaper- 

3 

ture.  That  is,  the  I  sensor  transfer  function  is  determined  by 
the  pupil  of  the  instrument,  rather  than  by  the  spatial  extent 
of  the  input  image. 

3 

Figure  2  shows  the  response  of  the  I  sensor  to 
two  drastically  differing  image  distributions.  These  results 
are  taken  from  extensive  simulation,  verified  by  systematic  ex¬ 
perimentation,  already  completed  at  AOA.  The  dashed  line  shows 
the  sensor  response  due  to  displacement  of  a  Gaussian-shaped  spot 
along  a  quadrant  detector  axis.  The  quadrant  detector  is  2.25 
spot  radii  square,  with  a  nutation  radius  of  one  spot  radius. 

The  broken  line  is  the  response  to  a  pair  of  spots  oriented  at 
45°  to  the  quad  detector  axes,  and  displaced  along  the  x-axis. 

The  spot  shapes  have  the  same  Gaussian  shapes,  and  the  system 
response  is  the  same  for  both  images  until  detector  edge  effects 
become  important  at  large  displacements.  Similar  analysis  for 
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Responsivity  of  the  I  sensor  to  two  different 
image  plane  intensity  distributions.  The  results 
are  obtained  using  a  high  level  computer  simyla- 
tion  developed  by  AOA  which  duplicates  the  1 
processing  steps  as  well  as  a  variety  of  possible 
input  optical  distributions.  The  2  distributions 
are  a  single  circularly  symmetric  peak  with  Gaus¬ 
sian  edges,  chosen  to  emulate  an  unresolved  image 
point,  and  a  double-lobed  "dumbbell"  distribution 
having  3  times  the  lateral  extent  along  the  dumb¬ 
bell  axis  and  the  single  spot  but  with  the  same 
Gaussian  edge  slopes.  The  dumbbell  is  oriented 
45°  to  the  quad  cell  axis  with  one  lobe-  in  quad¬ 
rant  1,  the  other  in  quadrant  3.  The  displacements 
occur  in  the  +x  direction  and  the  outputs  are  sim 
to  be  essentially  identical  for  both  distributions. 
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other  even  more  complex  shapes  show  similar  results:  the  I 
sensor  response  has  a  very  high  degree  of  insensitivity  to  image 
shape . 


3 

1.2.4  Summary.  The  I  sensor  has  a  large  variety  of 
operating  characteristics  that  make  it  a  very  competitive 
wavefront  sensing  technique.  There  include  the  following. 

(1)  Fundamental  signal-to-noise  ratio  (SNR) 
performance  comparable  to  that  of  the  AC  shearing  interferometer 
and  the  DC  quad  Hartmann;  all  three  approach  theoretically 
optimal  performance  limits. 

3 

(2)  The  internal  reference  in  the  I  sensor  allows 
great  flexibility  in  the  detailed  optical  design  of  a  specific 
instrument . 

(3)  With  the  use  of  proper  electronic  processing, 
nutation  modulation  provides  an  AC  carrier  in  a  highly  efficient 
way,  yields  operation  with  both  positive  and  negative  contrast 
objects,  and  gives  a  very  high  degree  of  system  gain  insensi¬ 
tivity  to  object  shape  and  extent.  In  addition,  the  nutation 
modulation  technique  does  not  chop  a  uniform  background,  which 
considerably  reduces  the  dynamic  range  problems  at  the  front  end 
processing  stages. 

3 

1.3  LSI  Sensor  Concept 

The  LS>  ,  which  is  the  baseline  system  described  in 

3 

this  study  for  a  visible  light  sensor,  is  based  on  the  I  sensor 
concept  studied  extensively  for  use  in  IR  adaptive  optics  systems 
where  a  relatively  small  number  of  subaperture  measurements  are 
required  for  compensation  of  atmospheric  turbulence  effects.  Al- 

3 

though  the  LSI  sensor  has  a  completely  different  and  more 

efficient  optical  configuration,  the  processing  algorithms  and 

3 

sensor  behavior  is  much  like  that  of  the  I  sensor.  Therefore, 
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many  of  the  sensor  characteristics  proven  to  be  of  importance 

3  3 

for  the  I  sensor  are  projected  to  carry  over  into  the  LSI 

performance.  These  characteristics  obtained  by  analysis, 

simulation  and  experimental  measurements  are  summarized  briefly 

here  and  are  described  in  detail  in  the  final  report  of 

Reference  4.  The  important  advantages  of  this  system  are  the 

3 

motivation  for  incorporating  them  into  the  LSI  sensor  design. 

1.4  Design  Goals 

As  discussed  in  Section  1.1,  any  alternative  technolo¬ 
gy  to  that  currently  employed  in  the  CIS  must  retain  high  speed 
and  must  approach  quantum  limited  sensitivity.  However,  if 
really  large  (e.g.,  greater  than  200  subapertures)  compensation 
systems  are  to  be  built,  some  alternative  to  fully  parallel 
sensing  and  processing  must  be  found. 

The  cost  of  expanding  the  current  CIS  in  a  parallel 
fashion  is,  no  doubt,  well  studied.  We  have  made  similar 
estimates  of  the  cost  and  complexity  of  constructing  a 

3 

completely  parallel  I  sensor,  using  the  most  current  electronic 
technology.  Assuming  either  a  set  of  small  PMT's  at  the  input 
or  an  intensified  quadrant  detector,  and  employing  some  ad¬ 
vanced  processing  components  such  as  the  Intel  2920  signal 
processing  chip,  the  detection  and  processing  stages  would  cost 
approximately  $1,500  per  channel.  The  channels  would  contain 
no  analog  circuitry  requiring  any  adjustment,  and  would  be 
relatively  compact  (perhaps  one  3"  X  5"  PC  board  per  processor). 
The  cost  for  a  large  system  is  thus  $300K  for  detection  and 
processing  alone. 

There  is  another  major  cost  component  for  large  systems, 
and  that  is  in  the  optics  required  for  pupil  division  into  the 
measurement  subapertures.  For  a  Hartmann-type  system,  even  a 


cleverly  designed  system  becomes  intricate  and  somewhat 
sensitive  (even  using  a  referencing  technique).  In  a  parallel 
shearing  interferometer,  the  problem  is  one  or  precise  regis¬ 
tration  of  800  individual  detectors  in  the  conjugate  pupil 
plane . 

From  these  considerations,  our  design  goals  have 
centered  on  the  following. 

(1)  Eliminating  the  complex  pupil  division  optics 
which  becomes  increasingly  cumbersome  with  scale. 

(2)  Employing  a  monolithic  detection  array  to  reduce 
the  problems  of  size,  cost,  and  alignment  sensitivities. 

(3)  Attempt  a  parallel/serial  hybrid  processing 
scheme  which  eliminates  the  need  for  a  strict  processor  per 
channel  correlation. 

In  the  following  section,  a  wavefront  sensor  design 
which  meets  each  of  these  criteria  is  presented.  The  detection 
system  is  based  on  the  use  of  intensified  diode  arrays  which 
also  act  as  pupil  dividers  for  the  system.  Furthermore,  the 
characteristics  of  the  array  elements  are  used  as  part  of  the 
processing  stage,  so  that  only  a  single,  fairly  simple  processor 
would  be  required  to  operate  a  full  256  subaperture  system.  A 
fivefold  to  tenfold  cost  reduction  should  be  realized  with  this 
approach,  with  complexity  at  the  256  subaperture  level  less  than 
that  for  a  21  subaperture  CIS. 

1.5  Parallel/Serial  System  Design 

Figure  3  is  a  simplified  layout  showing  the  basic 
features  of  the  proposed  optical  system.  At  the  input  to  the 
system,  the  input  signal  wavefront  is  combined  with  the  locally 
generated  reference  plane  wave,  giving  the  highly  desirable 
optical  and  electronic  self-calibration  capability.  The 
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Figure  3.  Simplified  optical  layout 


CAPTION  FOR  FIGURE  3 


Figure  3.  Simplified  optical  layout  for  the  proposed 
wavefront  sensing  system.  The  conventional 
Hartmann-type  pupil  divider  and  discrete 
quadrant  detector  array  has  been  replaced 
by  a  quadrant  divider  and  four  monolithic 
detector  arrays.  The  input  wavefront  is 
combined  with  a  reference  wave  and  passed 
through  a  pair  of  galvanometer  mirrors, 
which  impart  the  circular  nutation.  The 
primary  focus  occurs  at  the  tip  of  the 
prismatic  quadrant  divider,  which  reflects 
the  beam  into  four  detector  channels. 

Final  lenses  project  the  input  pupil  onto 
the  detector  array,  which  also  acts  as  a 
pupil  divider.  The  number  of  subapertures 
is  determined  solely  by  the  number  of 
detectors  in  each  array;  otherwise  the 
optical  arrangement  is  the  same  for  any 
size  system. 


21 


reference  is  modulated  in  the  time  domain  so  that  it  can  be 


discriminated  from  the  signal.  In  order  that  shot  noise  from 
the  reference  not  degrade  the  system  performance,  the  reference 
is  switched  on  for  one  nutation  cycle  out  of  every  ten  or  more, 
so  that  most  time  is  spent  observing  signal  power.  In  the  sys¬ 
tems  which  have  already  been  constructed  at  AOA,  a  ratio  of 
about  one  reference  cycle  for  every  32  signal  cycles  has  been 
used  with  good  result. 

The  circular  tilt  or  nutation  modulation  is  provided  by 
a  pair  of  sinusoidally  driven  resonant  galvanometer  mirrors 
locked  in  phase  quadrature.  AOA  has  developed  and  operated 
such  systems  at  10  kHz,  which  provides  ample  sampling  rate  to 
give  the  required  1  kHz  information  bandwidth. 

The  innovations  in  this  design  occur  in  the  next  stage. 
Instead  of  placing  a  quadrant  detector  in  the  image  plane,  a 
prismatic  quadrant  divider  is  placed  at  the  focus  of  the  full 
aperture.  This  element  acts  to  simultaneously  divide  all  sub¬ 
aperture  areas  into  four  channels  for  detection  and  processing. 

After  the  optical  divider,  the  four  quadrant  beams  are 
collimated  and  fall  onto  four  separate  detector  arrays.  Because 
the  arrays  are  in  a  plane  conjugate  to  the  input  pupil,  the 
individual  detector  elements  behave  as  subaperture  masks.  This 
is  illustrated  more  directly  in  Figure  4.  By  projecting  the 
detector  array  pattern  back  through  the  optical  system,  it  is 
seen  that  each  detector  maps  into  an  input  subaperture 

For  every  subaperture  there  are  four  conjugate  detec¬ 
tors,  one  in  each  quadrant  plane.  The  outputs  of  these  detec¬ 
tors  then  act  exactly  as  the  outputs  of  a  single  quad  cell  in 
the  focal  plane  of  the  selected  subaperture  would  in  a  conven¬ 
tional  Hartmann  system.  Therefore,  any  processing  techniques 
which  work  in  the  Hartmann  system  are  applicable  in  this 
configuration  also. 
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Figure  4.  Detector  arrays  used  as  pupil  dividers.  In  each 

quadrant  channel  a  single  detector  element  is  pro¬ 
jected  onto  the  input  pupil  plane  and  defines  a 
subaperture.  The  4  detector  elements,  one  from 
each  channel,  which  map  onto  the  same  input  sub¬ 
aperture,  form  a  quadrant  detector  which  may  then 
be  used  to  determine  the  wavefront  tilt  in  that 
subaperture . 
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Thus,  the  I  sensor  processing  algorithms  with 

3 

nutation  modulation  are  immediately  applicable  to  the  LSI 

3 

processor.  The  processing  algorithms  proposed  for  the  LSI 

3 

sensor  are  optimized  versions  of  the  infrared  I  sensor  algo¬ 
rithms  for  use  in  quantum  limited  systems. 


I 


2,0  THE  LSr3  SENSOR  CONCEPT 
3 

The  LSI  wavefront  sensor  promises  significantly  improved 
performance  and  lower  cost  than  first  generation  wavefront  sen¬ 
sors.  This  has  been  made  possible  by  several  basic  advances  in 
wavefront  sensor  technology  that  will  be  described  in  this  ohaDter. 

The  advances  in  performance  are  due  to  the  optical  config¬ 
uration  of  the  wavefront  sensor  which  combines  the  advantages  of 
3 

the  Hartmann,  I  and  shearing  interferometer  wavefront  sensors, 
and  the  signal  processing  algorithm  using  modern  estimation 

3 

theory.  Quantitatively  the  LSI  sensor  is  approximately  a  fac¬ 
tor  of  two  more  sensitive  than  first  generation  photon  noise 
limited  wavefront  sensors,  like  the  shearing  interferometer  and 

3 

I  sensor,  when  operated  under  identical  conditions.  This  fac¬ 
tor  of  two  is  explained  in  greater  detail  later. 

3 

2.1  LSI  Optical  configuration 

3 

The  LSI  wavefront  sensor  has  some  of  the  properties 

3 

of  a  Hartmann  sensor  and  the  I  sensor  modification  of  it,  and 
some  of  the  properties  of  a  shearing  interferometer.  Figures 

5a,  5b,  and  5c  are  optical  schematics  of  these  wavefront 

3  3 

sensors.  The  LSI  is  like  a  Hartmann  and  I  sensor  in  that  a 
quadrant  divider  is  used  to  measure  the  centroid  of  the  target 
to  determine  subaperture  wavefront  tilt.  It  is  similar  to  the 
shearing  interferometer  in  that  the  subapertures  are  defined  at 

3 

a  conjugate  pupil  plane.  Where  the  I  sensor  required  a  refer¬ 
ence  beam  to  avoid  the  catastrophic  problem  of  aligning  n  quad 
cells  with  their  respective  subaperture  optical  axes  as  in  the 

3 

traditional  Hartmann  sensor,  the  LSI  sensor  uses  a  single  4- 
sided  pyramid  divider  instead  of  n  quad  cells,  and  avoids  that 
alignment  problem  entirely.  By  using  a  quadrant  scheme  for 

3 

centroid  location  as  in  the  I  sensor,  the  same  photons  are 


25 


1 


PUPIL  PLANE  DIVISION  INTO  SUBAPERTURES 


T> 

Q> 


QUAD  CELL  WAVEFRONT  TILT  DETECTOR 


(a)  Hartmann  sensor. 
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Figure  5.  Hartmann,  SI  and  LSI  sensors.. 
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used  for  x  tilt  and  y  tilt  measurements  resulting  in  greater 

sensitivity  than  the  shearing  interferometer;  this  is  explained 

3 

in  detail  later  in  this  chapter.  In  short,  the  LSI  sensor  of¬ 
fers  the  advantages  of  both  the  Hartmann  and  shearing  interfero¬ 
meter  sensors  without  the  disadvantages  that  are  peculiar  to 
either. 

Light  from  the  target  goes  through  a  nutator,  PI,  and  is 
focused  onto  the  pyramid.  From  the  pyramid,  the  light  is  recol¬ 
limated  by  a  lens.  The  detector  array,  which  defines  the  sub¬ 
apertures,  is  placed  behind  the  collimator.  In  the  case  of  Fig¬ 
ures  5b  and  5c,  the  pupil  division  is  performed  at  a  plane 
that  is  a  conjugate  of  the  PI  plane  in  front  of  the  first  lens. 

The  nutator  is  shown  schematically  in  Figure  5c  as  a  ro¬ 
tating  wedge  of  glass  (in  practice  galvanometer  nutators  are 

3 

used, as  in  the  I  sensor).  The  effect  of  the  nutator  is  to  move 
the  image  of  the  target  in  a  circular  pattern  at  plane  P2 ,  where 
the  four-sided  pyramid  is  located.  This  nutation  modulates  the 
intensity  of  the  light  arriving  at  the  detectors.  The  output 
of  the  detectors  is  then  demodulated  to  give  the  subaperture 
wavefront  tilt.  The  nutation  in  a  practical  wavefront  sensor 
is  performed  by  two  resonant  galvanometer  driven  mirrors  (see 
Section  6)  located  near  a  conjugate  image  of  the  primary.  The 

3 

nutation  radius,  the  LSI  sensor's  version  of  shear,  is  elec¬ 
tronically  controlled  by  varying  the  drive  voltage  to  the  gal¬ 
vanometers  . 

A  simplified  description  of  the  operation  of  the  wavefront 
sensor  may  be  understood  by  masking  the  PI  plane  in  Figure  5c 
to  a  single  subaperture  and  ignoring  the  nutator.  Figure  6 
shows  how  wavefront  tilt  is  measured.  Light  from  a  single  sub¬ 
aperture  is  focused  on  the  tip  of  the  pyramid;  light  from  the 
central  subaperture  will  only  fall  on  the  central  detector  in 
each  of  the  four  detector  arrays.  Also  note  that  if  the  sub- 
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Figure  6.  Measurement  of  wavefront  tilt. 


aperture  wavefront  is  tilted,  its  image  will  be  shifted  at  the 
plane  of  the  pyramid.  As  a  result,  more  light  will  be  detected 
by  the  detector  on  the  right  side  than  the  left,  as  shown  in 
Figure  6. 

3  3 

The  difference  between  the  LSI  and  the  I  sensor  is  clear. 
3 

In  the  I  sensor,  as  in  the  Hartmann  sensor  of  Figure  5a,  the 

wavefront  is  first  divided  into  n  subapertures  by  an  optical 

divider  and  each  subaperture  tilt  measured  by  four  detector  ele- 

3 

ments  in  a  quad  cell  configuration.  In  the  LSI  sensor,  the 

quad  cell  measurement  is  provided  by  a  single  4-sided  optical 

divider  and  then  the  wavefront  is  divided  into  subapertures  for 

detection  by  4  detector  elements  for  each  subaperture.  In  the 
3 

LSI  configuration,  the  detector  elements  are  not  clumped  to¬ 
gether  as  quad  cells;  therefore  a  uniformly  spaced  detector 

array  may  be  utilized.  However,  the  signal  processing  algo- 

3  3 

rithms  are  similar  in  the  I  and  LSI  since  quad  cell  centroid- 
ing  and  nutation  modulation  are  used  for  the  basic  subaperture 

3 

tilt  measurements.  This  is  the  principle  behind  the  LSI  sen¬ 
sor  . 


2.2  Signal  Processing 

Without  nutation,  the  signal  processing  algorithm  is 
simply  the  Hartmann  quad  cell  algorithm. 


Tilti  =  DI,i  +  DIII , i  ~  DII,i  ~  DIV , i 
DI,i  +  DII , i  +  DII I , i  +  DIV,i 


(1  ) 


where  tilt  i  is  the  tilt  in  the  i'th  subaperture,  D  is  the 

**  9  * 

output  of  the  yth  detector  in  the  xth  array;  a  similar  formula 
is  used  for  tilt  in  the  orthogonal  direction. 
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With  nutation,  the  detector  outputs  are  a  function  of  nuta¬ 
tion  radius  and  phase.  The  outputs  of  four  detectors  over  a 
complete  nutation  cycle  are  needed  to  derive  subaperture  tilt. 
The  basic  concept  is  the  same  as  equation  1.  The  numerator  of 
equation  1  is  a  term  that  is  linearly  related  to  the  wavefront 
tilt  for  small  tilts  and  linearly  related  to  target  brightness. 
The  denominator  is  independent  of  tilt  and  linearly  related  to 
target  brightness. 

The  ratio  is  then  linearly  related  to  tilt  and  independent 
of  target  brightness. 

The  formula  for  demodulation  with  nutation  is 

2  4 

/  ZD  (0)NT(6)d0 
o  i  =  i'J  J 

tnt  =  - (2) 

/  ZD  (0)P  (6)de 
o j=lJ  J 

where  Dj(0)  is  the  Jth  detector  output  as  a  function  of  nuta¬ 
tion  phase,  Pj(9)  is  a  correlation  function  that  extracts  total 
power,  Nj(0)  is  a  correlation  function  that  extracts  tilt  and 
power . 

A  number  of  different  P j ' s  and  Nj's  were  investigated. 

These  properties  are  discussed  in  Section  6.  It  was  decided 
that,  because  a  number  of  demodulation  algorithms  worked,  it 
might  be  profitable  to  use  modern  estimation  theory  to  find  an 
optimal  algorithm,  given  target  size  and  shape  and  a  statistical 
description  of  the  noise.  Photon  noise  limited  and  detector 
noise  limited  conditions  were  investigated. 

The  major  result  is  that  with  the  optimized  algorithms 
the  LSI3  sensor,  for  most  target  scenarios,  is  predicted  to 
have  a  factor  of  two  higher  phase  accuracy  for  a  given  target 
signal  power  and  equivalent  detectors  than  the  theoretical 
performance  of  the  shearing  interferometer  CIS.  The  details 

of  these  calculations  are  given  in  Section  3. 
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2.3  Detector  Selection 

The  AWS  optical  configuration,  in  which  the  quadrant 

or  reticle  elements  precede  the  detectors,  provides  two  unique 

virtues  for  the  design  of  the  detectors:  the  optical  system 

is  independent  of  the  detector  configuration  except  for  the 

final  collecting  optics;  the  optical  system  is  independent  of 

the  number  of  subapertures  into  which  the  pupil  is  divided. 

As  a  result,  both  the  detector  and  optical  design  can  be 

optimized  separately.  This  is  not  true  with  a  quad  Hartmann 

3 

type  design  of  the  I  sensor,  in  which  the  detector  arrange¬ 
ments  and  subapertures  dictate  the  optical  arrangement.  The 

interface  between  the  sensor  optics  and  detectors  is  an  f 

imaging  system  which  divides  the  pupil  images  into  subaper¬ 
tures.  The  final  AWS  optical  design  presented  in  Section  4 
creates  16  mm  diameter  images  of  the  pupil. 

This  pupil  is  imaged  directly  on  arrays  which  then 
divide  the  pupil.  The  arrays  that  can  be  used  are 

(1)  area  CCD's  or  self  scanned  diode  area 
arrays  preceded  by  intensifiers  for  direct  detection  and  sub¬ 
aperture  division,  or 

(2)  lens  arrays  which  condense  the  ‘ ight  in 
the  subaperture  to  an  array  of  small  points  which  can  be  used 
as  input  to  optical  fibers,  or  which  can  be  made  to  match  de¬ 
tectors  directly  or  through  an  intensifying  device. 

The  fibers  can  be  spread  out  to  large  detectors  or 
rearranged  to  convert  from  area  to  a  linear  format  to  permit 
the  use  of  small  linear  CCD  or  self  scanned  diode  arrays. 

The  results  of  the  study  produced  three  recommended 
detector  systems  for  a  256  element  array,  and  a  photomultiplier 
array  for  a  reduced  scale  21  subaperture  breadboard.  These 
systems  and  their  spec i f i cat i ons  in  order  of  preference  are 
as  follows. 

I 
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(1)  EB  mode  detector,  internal  optics 

(a)  Characteristics  (see  Figure  47) 

-  Intensification  with  a  256  element 
linear  self  scanned  diode  array 

-  Optical  input;  lenticular  array, 
fiber  optic  area  linear  converter 

-  Internal  ellipsoid  focusing  fiber 
optic  output  on  opaque  GaAs  cathode 

-  Electron  optics;  4:1  reduction  and 
15  -  20  kV 

-  Reticon  256  element  array;  10  MHz 
sample  rate,  10  KHz  frame  rate  for 
processed  tilts 

(b)  Performance 

-  Highest  possible  efficiency  in  using 
solar  spectrum  (16.8%  of  all  solar 
photons ) 

(2)  EB  mode  detector,  external  optics 

(a)  Characteristics  (see  Figures  40a,b,e) 

-  Same  as  (1),  but  uses  external  op¬ 
tics  for  input  to  opaque  cathode 

(b)  Performance 

-  Less  than  (1)  because'  optics  are 
1  ess  e f f ici en t 


-  Straightforward  intensifier  design, 

2  to  3  coupled  stages  with  P-24 
phosphors 

-  Transparent  GaAs  photocathodes 

-  Fiber  optically  coupled  CCD 
buried  channel  array 

-  Limited  array  size  by  selective 
gating;  samples  taken  in  32  x  36 
element  blocks 

-  All  elements  commercially  available, 
no  special  design  required 

(b)  Performance 

-  15%  of  solar  flux  utilized,  more 
efficient  input  compensates  for  lower 
cathode  sensitivity 

-  Poorer  accuracy  because  of  large 
noise  factor  (1.2) 

(4)  Breadboard  (21  channel) 

(a)  Characteristics 

-  PMR's  1/2"  diameter 

-  Lenticular  array  input  to  fiber 
optics  which  spread  out  to  tubes 

-  Photon  counting  for  low  noise  fac¬ 
tor 

(b)  Performance 

-  Bi  alkali  photocathode  comparable 
with  CIS,  potentially  5%  of  solar 
photons  collected 
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The  reasons  given  for  these  choices  and  a  large 
variety  of  alternate  schemes  are  given  in  Section  6  and  the 
above  performance  data  is  summarized  in  Table  11. 

The  major  advantages  of  the  designs  presented  here 
are  as  follows. 

(1)  Intensified  detectors  provide  photon 
limited  signals  for  the  highest  theoretical  accuracy. 

(2)  The  ability  of  the  LSI^  sensor  to  use  the 
wide  solar  spectra  without  loss  of  accuracy  due  to  a  loss  of 
visibility  permits  the  use  of  the  new  wide  band  GaAs  cathodes. 

3 

The  variance  for  the  best  LSI  over  that  of  the  best  shearing 
interferometer  is  a  factor  of  2.76  smaller. 

(3)  The  combined  improvement  over  theoretical 
CIS  due  to  the  algorithm  and  the  efficiency  of  detection  is 
equivalent  to  1.63  magnitudes  increased  sensitivity. 

2.4  Signal  Processing  Electronics 

O 

The  LSI°  optical  system  requires  four  separate  detec¬ 
tor  arrays  for  each  subaperture  output  as  a  quadrant  output 
signal.  In  each  detector  array  there  will  be  at  least  one 
detector  element  per  subaperture.  For  the  large  scale  sub¬ 
aperture  division,  up  to  256,  the  detector  systems  selected  all 
use  integrated  circuit  detector  arrays. 

These  detector  arrays  accumulate  photoe 1 oct rons  and 
this  charge  accumulation  replaces  analog  circuit  integrators 
for  signal  processing.  Furthermore,  the  integrated  arrays 
take  the  subaperture  optical  data  and  transform  it  into  10  MHz 
serial  electronic  data.  This  charge  accumulation  and  parallel- 
serial  conversion  reduces  the  analog  circuitry  requirements 
from  four  preamp  I i f i er- i ntegra tors  per  subaperture  (  1024  per 
system)  to  1  preamplifier  per  array  (4  per  system)  for  up  to 
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a  256  subaperture  system.  These  4  preamplifiers  are  each 
followed  by  a  10  MHz  pipelined  digital  processor  and  the 
4  parallel  outputs  are  then  processed  by  a  single  processor 
programmed  with  the  optimum  algorithm.  The  output  from  the 
processor  iz  a  serial  stream  of  tip/tilt  error  signals  from 
each  subaperture. 

3 

The  LSI  sensor  electronic  processor  design  is 
intended  to  minimize  operator  intervention.  Novel  features 

3 

of  the  LSI  are  exploited  to  permit  automatic  system  calibra¬ 
tion,  extensive  self  test,  and  automatic  dynamic  alignment 
during  system  operation.  These  extraordinary  capabilities 
are  in  addition  to  a  signal  processor  architecture  which 
totally  eliminates  subaperture  level  electronic  components 
and  adjustment.  This  proprietary  architecture  results  in  a 
considerable  reduction  in  circuitry  compared  to  the  tradition¬ 
al  processor  per  subaperture  approach.  The  details  of  the 
design  are  described  in  Section  5. 


3,0  SIGNAL  PROCESSING  AND  PERFORMANCE  EVALUATION 
3 . 1  Introduct ion 

This  section  describes  the  various  algorithms  for  de¬ 
riving  subaperture  tilts  from  the  output  of  the  detectors.  Sev¬ 
eral  algorithms  were  developed  in  the  course  of  this  program. 

3 

These  algorithms,  in  addition  to  the  original  I  sensor  algo¬ 
rithm  (the  SCS  algorithm)  were  examined  as  part  of  this  pro¬ 
gram.  In  order  to  facilitate  the  comparison  of  these  algo¬ 
rithms  with  each  other  and  with  the  AC  shearing  interferometer, 
a  uniform  mathematical  notation  is  adopted. 

3 

The  optical  layout  of  the  I  sensor  has  been  described 
in  previous  sections.  Basically,  a  divider  is  placed  in  a  plane 
where  an  image  of  the  target  would  be  formed.  After  the  di¬ 
vider,  the  light  is  recollimated  and  a  detector  array  is  placed 
in  a  plane  that  is  conjugate  to  the  pupil  of  the  optical  system. 
For  each  subaperture  there  are  2  to  4  detectors,  one  for  each 
facet  of  the  divider.  The  outputs  of  the  detectors  for  each 

subaperture  are  demodulated  to  find  the  centroid  of  the  image 

3 

in  the  divider  plane.  The  difference  between  the  various  1 
algorithms  and  the  shearing  interferometer  is  that  each  algo¬ 
rithm  is  a  different  mathematical  description  of  the  centroid. 

Because  of  the  similarity  of  the  various  wavefront 
sensors,  a  single  model  for  all  of  them  is  easily  constructed. 
Figure  7  is  an  optical  schematic  of  a  generalized  wavefront 

sensor.  For  a  shearing  interferometer,  the  mask  is  a  Ronchi 

3 

grating.  For  an  LSI  sensor,  each  detector  array  has  a  differ¬ 
ent  mask.  The  mask  transmits  the  light  in  one  quadrant  of  the 
plane.  The  target  is  imaged  onto  the  reticle.  The  different 
type  sensors  have  different  reticles  and  different  demodulation 
algorithms . 

The  image  on  the  reticle  is  a  convolution  of  the  true 
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Figure  7.  Optical  schematic. 
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object,  the  diffraction  limited  point  spread  function  of  the 

telescope,  and  the  atmospheric  point  spread  function.  We  will 

3 

now  describe  the  AC  shearing  interferometer,  the  LSI  sensor, 

3 

and  the  MSI  sensor. 


3  3 

3.2  LSI  ,  MSI  and  Shearing  Interferometer 


For  all  these  types  of  wavefront  sensors,  the  output 
of  the  detector  is  eq .  3. 


Detector  output  =  / dx  / dy  M( x-f , y-g)T( x , y ) 


where  T(x,y)  is  the  effective  irradiance  distribution  of  the 
target  (convolution  of  true  object  and  atmospheric  point  spread 
function  (PSF)  and  diffraction  PSF) ;  M(x,y)  is  the  mask  func¬ 
tion;  f,g  are  the  modulation  functions. 

Table  1  lists  the  types  of  masks  and  modulation 
functions  for  the  three  types  of  wavefront  sensors. 

Equation  3  is  different  from  the  traditional  formula 
for  a  shearing  interferometer  derived  using  Fourier  optics.  In 
effect,  it  says  that  the  difference  between  the  various  sensors 
are  the  masks  and  modulation  patterns,  and  that  these  are  the 
only  differences,  since  the  subapertures  are  divided  in  the  same 
way.  Hence,  difference  in  the  signal-to  noise  ratios  (SNR)  of 
the  various  sensors  may  be  accurately  derived  from  equation  3. 

3 

In  an  LSI  system,  there  are  four  detectors  per  sub¬ 
aperture.  Because  of  nutation,  the  outputs  of  these  detectors 
change  with  time.  This  is  also  true  for  the  shearing  interfero¬ 
meter.  In  addition  to  different  mask  functions  and  modulation 
functions,  the  demodulation  algorithms  for  the  shearing  inter- 

3 

ferometer  and  I  sensors  are  different.  For  a  shearing  inter¬ 
ferometer,  the  centroid  is  at  zero  phase  for  a  single  spatial 
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TABLE  1.  MASKS  AND  MODULATION  FUNCTIONS. 


co 
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frequency  component  of  the  target.  Phase  is  defined  by  equation 

4. 


A 

x 


4  2 Ti¬ 


tan 


l  /  S.O. d( 
-1  i=l  °  11 


4  2 ?? 


L  f  C.0.d9 
i=l  °  11 


(4) 


where  CK  represent  the  outputs  of  the  4  detectors,  S ^  is  a  cor¬ 
relator  function  which  extracts  the  sine  of  the  phase  from  the 
detector  output,  and  the  cosine  extraction. 

For  a  conventional  radial  grating  interferometer,  sep¬ 
arate  detectors  are  used  for  x  tilt  and  y  tilt  and  hence  for  <p 
two  of  the  correlation  functions  are  0  for  x  tilt  and  the  other 
two  are  0  for  y  tilt. 

In  an  I  sensor,  the  tilt  is  defined  by  equation  5. 

4  2t? 

.1  iHw6 


tilt  = 


i=r 


4  2tt 

l  27?  ^  D  j°  j  dG 
o  1  1 


(5) 


where  CK  is  the  output  of  the  detectors,  ;  and  are  correla¬ 
tion  functions.  As  explained  in  Section  2,  the  numerator  is  a 
term  that  is  proportional  to  the  tilt  and  target  brightness, 
while  the  denominator  is  a  term  that  is  proportional  to  target 
brightness . 

Table  2  shows  the  various  demodulation  functions 
for  the  three  wavefront  sensors. 
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TABLE  2.  DEMODULATION  FUNCTIONS. 


Sensor 


Numerator 
funct ion 
for  x  tilt 


Denominator 

function 


Shearing  Inter.  sin  9 

(Ronchi  grating) 

MSI3  sin  6 


cos  6 

sin  26  both  x  &  v 


3 

LSI  There  are  4  different  demodulation 

algorithms  for  the  LSI3  optical 

configuration 

See  text  for  details 


DC  (x  tilt) 

SCS  (original  I3) 


N] 

Nr 

N; 


N . 


=  Sin  6 
-  Sin  6 


N .  =  +1 

l 


D^  =  Sin  6 
D2  *  N. 

D_  =  -  Sin  0 

»!  ■  -3 


MS  (This  algorithm  is  equivalent 

Nx  =  Sin 


N4  "  -N1 


3 

to  the  MSI  sensor) 

6  Dj  =  D3  =  Sin  26 

D  =  D  =  -  Sin  26 
2  4 


ACDC 


N1  1  -  cos  0  Dj^  =  1 

N0  =  -1  +  cos  6  D_  =  1 

A  - 1 

N.  =  N.  D  =  1 

4  1  4 


These  correlation  functions  are  optimal  in  the  sense 
that  for  a  Gaussian  object  and  nutation  radius  equal 
to  the  width  of  the  object. 

The  noise  in  the  derived  tilt  is  a  minimum  for  a  given 
level  of  photon  noise. 


'■<■3  Performance  Evaluation 

Once  we  have  a  description  of  the  various  wavefront 
sensors  (equations  3,  4,  and  5  ),  it  is  possible  to  evalu¬ 

ate  the  relative  SNR’s  of  the  sensors. 

In  order  to  calculate  the  noise  equivalent  tilt  of  a 
sensor,  we  must  determine  the  responsivity  of  the  sensor,  the 
change  in  the  derived  tilt  (equation  ■}  or  )  for  a  given 
change  in  the  centroid  of  the  target,  as  well  as  the  expected 
fluctuations  in  the  derived  tilt  due  to  measurement  noise. 

When  the  target  is  centered,  the  numerator  terms  of  equations 
■1  and  3  are  zero.  For  small  changes  in  the  target  position 
the  numerator  term  will  change  significantly,  while  in  first 
order,  the  denominator  term  will  not  change.  As  a  result,  all 
of  the  following  calculations  will  concentrate  on  the  numerator 
term.  For  practical  purposes,  the  SNR  of  the  numerator  is  the 
SNR  of  the  ratio.  This  is  shown  below. 

aRpX  =  (l/D)ONpX)  -  (N/D2)OD/^X)  (6) 

or  (1/D)(3N/3X) 

dd 

where  R  is  the  ratio  of  two  numbers,  N  or  0  (small),  ^  or  0 
(small),  and  x  is  the  displacement  of  the  target  centroid. 

Equation  (5  shows  that  the  responsivity  of  the  ratio 
R  is  dependent  primarily  on  the  first  term  in  equation  (> ,  the 
second  term  being  very  small.  Similarly,  the  fluctuations  in 
R  in  equation  7  due  to  measurement  noise  are  dependent  pri- 
marily  on  the  numerator. 

aR  /an  =  (1/D)  (>N/>n)  -  (N/D2 ) O D/3 n )  (71 

cL 

where  n  is  a  random  variable  of  which  N  and  D  are  functions. 

Since  N  is  small,  only  the  first  term  in  equation  7 
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is  important.  The  variance  of  R&  due  to  fluctuations  in  n  are 
then  equation  8. 


a 


P. 


1 

D 


f^ON/an)2  <(An)2> 
Ln 


1/2 


(8) 


where  An  are  the  fluctuations  due  to  noise  and  the  summation  is 
over  the  random  variables.  Consequently,  the  SNR  of  Ra  ratio  is 
equal  to  the  SNR  of  the  numerator. 

The  numerators  of  equations  4  and  5  have  the  form 
of  equation  9. 


4  2  TT 

Num  =  l  (1/2tt)  /  d9  C  (0)0. (G)  (9) 

r.=l  o  1  1 


where  C.  is  the  correlation  function  for  the  ith  detector,  0.  is 
the  output  of  the  ith  detector. 

In  the  photon  noise  limited  case,  equation  10  is 
true;  in  a  detector  noise  limited  case,  H  is  true. 

a02(6)  =  0(6)  (10) 


where  the  variance  of  0(6)  is  0(6)  for  photon  noise. 
2 

a  (9)  =  constant 
o 


(ID 
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for  detector  noise  limited  operation. 

Using  standard  error  propagation  formulae,  the  SNR  is 
Equation  12  for  photon  noise. 


(SNR)  = 


[l  57/  ci(e)R1(6)dej 


/'  57/  c  2(8)o  <?)df 

O 


(  1  2  ) 


where  CL(6)  is  defined  by  equation  3  ,  and  (0)  is  defined  by 
equation  13. 


OO  00 

R(9)  =  fdx  fdv  ?tjx-  f(  *  )  ,  y-g(  0  )~]  7T(x,v)  •  x“  (13) 

— 00  _oo  *”  -* 


The  difference  between  3  and  12  is  that  T(x,y)  the 
effective  target  irradiance  distribution  in  equation  3  is  re¬ 
placed  by  the  slope  of  the  irradiance  distribution.  The  numera¬ 
tor  of  equation  12  is  the  square  of  the  responsivi ty ,  while  the 
denominator  is  the  variance  due  to  photon  noise.  In  the  follow¬ 
ing  graphs  and  tables,  the  data  was  generated  from 
equations  3  $■  13.  The  integrals  were  evaluated  numerically. 

The  assumptions  for  the  calculations  are  listed  in  the 

figures.  A  Hartmann  sensor  is  also  evaluated.  The  main  diflev- 

3 

ence  between  the  Hartmann  sensor  and  the  LSI  sensor  is  that  t he 
effective  target  irradiance  distribution  in  the  Hartmann  system 
is  the  convolution  of  the  subaperture  diffraction  pattern  and 
the  true  target  while  for  the  LSI^  and  shearing  interferometer, 
the  effective  shape  is  the  convolution  of  the  full  aperture  cl  i  f  — 


M 


1 


I 


fraction  pattern  and  the  true  target. 

The  atmosphere  was  ignored  in  these  calculations 
for  two  reasons.  One  was  to  simplify  the  calculations.  The 
second  was  to  make  the  comparison  with  other  published  data* 
which  compared  the  shearing  interferometer  to  the  Hartmann 
sensor. 

The  effect  of  the  atmosphere  is  to  limit  the 

us>able  shear  on  the  shearing  interferometer  to  approximately 

r  ,  the  atmospheric  phase  coherence  length.  The  result  of 

including  the  atmosphere  in  calculating  the  data  for  the 

3 

figures  would  be  to  reduce  the  advantage  of  the  LSI  and  SI 
over  the  Hartmann  sensor  for  very  small  targets  or  very 
distant  targets.  We  repeat  that  the  results  in  this  report 
do  not  conflict  with  other  published  data*. 

3.4  Optimal  Correlator  and  Object  Size  Determination 

While  evaluating  the  relative  SNR  of  several  I3 
demodulation  algorithms,  it  was  decided  that  it  would  be 
profitable  to  analytically  determine  an  optimal  correlation 
function  for  equations  4  and  5.  Here  we  use  the  word  optimal 
correlator  to  mean  a  function  C^(0)  that  will  result  in  a 
minimum  variance  tilt  estimate  given  a  known  noise  in  the 
detector  output  0(6). 

For  photon  noise  limited  operation,  the  optimal 
correlation  is  found  by  setting  the  variational  derivative 
of  equation  12  to  zero  and  solving  for  the  correlation 
function  CL (8).  The  answer  is  eouation  14. 

R,(0) 

Optimal  C.(0)  =  -  (14) 

1  0.(0) 


Reference  will  be  made  available  to  qualified  military  and 
government  agencies  on  request  from  RADC  (OCSE)  Grifliss 
AFB  NY  13441. 
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where  R  is  defined  by  equation 

d . 


Id  and  0  is  delined  by  equation 


In  all  the  figures  in  this  sect  ion*  the  graphs  labeled 
optimal  1^  use  the  correlation  function  defined  by  equation 
M.  Note  that  the  optimal  correlator  change's  with  object 
size,  shape  as  well  as  nutation  radius. 


Another  topic  that  was  studied  as  part  of  the'  AWS  pro¬ 
gram  was  the  measurement  of  object,  size  by  proper  demodulation 
of  the  outputs  of  the  wavefront  sensor  detectors.  The'  effective 
object  size  of  a  combination  of  atmospheric  seeing  blur  and  true 
object  size  is  an  important  parameter  in  a  shearing  interfer¬ 
ometer  wavefront  sensor  because  the  highest  tilt  accuracy  is 
achieved  when  the  shear  is  matched  to  the  object  size,  when  the 
grating  period  was  twice  the  object  size.  While  Figure  !’ 
shows  that  for  a  Gaussian  object,  the  accuracy  of  a  tilt  meas¬ 
urement  is  not  a  strong  function  of  nutation  radius  when  an  op¬ 
timal  correlator  is  used,  the  measurement  of  object  size  is 
still  important.  The  reason  is  that  some  of  the  I  cubed  demod¬ 
ulation  algorithms  are  object  size  dependent,  in  the  sense  that 
the  derived  lilt  in  radians  is  a  function  of  the  angular  size 
of  the  object.  Table  d  lists  the  size  dependence'  of  various 
algorithms . 


*  F  i  gures  7  1  li  r<  High  1  2 


Ki^tirc  S. 


NUTATION  RADIUS 


Fipcure  10.  Relative  SNR  vs.  nutation  radius  for  a 
square  object. 
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V  RADIANS) 


CIRCULAR  OBJECT,  4  ;j RAD  DIA. 
ATMOSPHERE  PSF  IGNORED 
DIFFRACTION  PSF  -<  4  -RAD 


N  PHOTON'S 


Till  noise  ns  a  function  of  number  of 
delected  photons. 
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F  i  iru  re  11 


As  an  illustration,  we  note  that  in  the  classical  quad 
cell  Hartmann  system  CDC  numerator  and  denominator)  the  ratio  is 
proportional  to  1/W.  In  other  words,  if  the  wavefront  tilt  is  1 
yrad  and  the  width  of  the  object  is  10  prad,  the  output  of  the 
DC  algorithm  would  be  ^  2/10;  if  the  object  size  was  100  prad,  a 
1  ;  rad  tilt  would  give  an  output  of  2/100.  On  the  other  hand, 
use  of  the  DC  numerator  and  SCS  denominator  would  produce  a  num¬ 
ber  that  was  independent  of  object  size,  depending  only  on  nuta¬ 
tion  radius.  This  is  the  case  with  the  shearing  interferometer 
where  the  output  represents  a  given  wavefront  tilt  that  de¬ 
pends  on  the  shear  but  not  on  object  size.  In  general,  the  size 
of  the  object  will  not  change  on  a  fast  line  scale,  e.g.  10kHz. 
As  a  result,  the  object  size  may  be  measured  with  very  high  ac¬ 
curacy  using  several  tens  of  milliseconds  of  data.  The  pre¬ 
ferred  method  of  measuring  size  is  to  use  the  SCS  denominator 
and  the  MS  denominator. 

Object  size  may  in  theory  be  measured  with  a  shearing 
interferometer  by  monitoring  fringe  visibility.  This  technique, 
however,  has  a  drawback.  The  fringe  visibility  formula  that  is 
usually  presented  in  descriptions  of  the  shearing  interferometer 
ignores  the  effect  of  scintillation.  Because  of  this  fringe 
visibility,  measurements  are  unreliable.  Currie  et  al  .*  have 
solved  the  problem  by  measuring  the  apparent  visibility  as  well 
as  the  scintillation  and  correcting  the  visibility  results  in 
the  data  analysis.  The  problem  is  due  to  the  fact  that  fringe 
visibility  is  in  effect  the  ratio  of  the  AC  optical  power  to 
the  DC  optical  power.  Modulating  at  10kHz  there  are  no  scintil¬ 
lation  effects  at  10kHz  but  there  are  at  DC. 

3 

With  an  I  sensor,  which  compares  the  AC  power  at 
10kHz  (SCS  denominator)  to  AC  power  at  20  kHz  (MS  denominator), 
scintillation  is  not  a  significant  source  of  error. 


♦Currie  et  al.,  Astrophys.  .J.  1 87 ,  131(  1974). 


4.0  OPTICAL  SYSTEM  DESIGN 


The  design  of  the  AWS  optical  system  is  less  demanding  than 
might  be  expected  because  the  instrument  uses  a  reference  beam  to 
define  a  perfect  wavefront.  By  measuring  the  difference  between 
the  reference  and  target  point  positions,  the  errors  of  the  com¬ 
mon  mode  optics  are  removed  from  the  measurement  of  the  wave- 
front.  In  consequence,  it  would  seem  that  any  optical  system 
would  be  adequate.  However,  the  effect  of  an  inferior  optical 
system  will  be  to  reduce  the  performance  in  detecting  and  cor¬ 
recting  the  wavefront  for  small  objects.  For  a  2  arc  second  ob¬ 
ject  (ot0),  the  response  will  be  reduced  to  90%  if  the  standard 
deviation  of  the  spread  function  a  is  0.86  arc  seconds. 

If  we  assume  Gaussian  spread  functions  and  compute  the 

2  2 

change  in  peak  intensity  to  be  i(P)  =  (cr0/aSF)  ,  where  aQ  is  the 

diffraction  limited  spread  function  which  is  given  approximately 
2  2 

by  (0 . 375^1/D)  ,  and  where  ogF  is  the  variance  of  the  allowed 
spread  function  for  90%  response  to  a  target,  we  can  use  the  ex¬ 
pression  for  peak  intensity  or  the  Strehl  Ratio,  to  yield  a  re¬ 
quired  rms  phase  error  for  the  sensor  optical  system.  Thus 

i(t>)  =  (o0/ogF)2  =  ^  1  -  ^(2tt/X)2(AX)2]  j  =  e~(  2tt  AX  /  X  ) 

(15) 

or 

4 

Jo(2ttAX/X)  for  large  arguments 

for  X  =  0.64  pm  (because  of  the  use  of  GaAs  cathodes)  at  D  = 

1.6  m,  i(P)  =  0.0304,  and  the  rms  phase  error  AX  allowed  is 
X/3.3.  The  allowed  peak  to  peak  error  is  about  one  wave. 

This  is  not  a  difficult  requirement  for  an  f/16  optical  sys¬ 
tem,  and  it  will  permit  some  latitude  in  handling  difficulties 
due  to  nutation,  pupil  imagery  and  chromatic  aberrations.  The 
optical  systems  that  have  been  designed  began  as  refracting  sys- 
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terns.  The  evolution  of  the  detector  design  has  extended  the 
range  of  colors  to  such  an  extent  that  a  reflecting  system  became 
desirable. 

The  first  problem  is  to  design  a  system  which  prevents  sub¬ 
aperture  smearing. 

The  nutation  device  is  the  major  geometrical  constraint  on 
the  system.  In  order  to  operate  at  10  kHz,  the  beryllium  mir¬ 
rors  are  limited  to  dimensions  of  approximately  5  x  7  mm.  For  a 
system  to  have  no  nutation-induced  subaperture  smearing,  the  tel¬ 
escope  primary  should  be  imaged  onto  the  nutator  mirrors  and  then 
onto  the  detector  array,  while  the  object  should  be  imaged  onto 
the  divider  plane.  These  operations  are  accomplished  by  the  con¬ 
figuration  in  Figure  13.  In  this  type  of  system,  the  maximum 
number  of  resolved  subapertures  is  diffraction  limited.  If  f/4 
lenses  are  used  with  5  mm  mirrors,  the  maximum  number  of  subaper¬ 
tures  is  approximately  (5  mm/1.2  if)2  -  1.6  x  107. 


This  is  not  an  important  constraint,  therefore,  and  we  may 
consider  simpler  optics, 

The  configuration  of  Figure  14  shows  the  nutation  mirrors 
displaced  from  conjugate  pupil  planes,  eliminating  a  lens.  The 
primary  is  considered  to  focus  at  LI  with  a  speed  of  f/20.  The 
focal  length  of  L3  is  chosen  to  match  the  detector  size  to  pupil 
size . 

For  a  typical  system,  assume  the  following. 

Cl)  A  10  cm  subaperture  size,  which  gives  an  angular  nuta¬ 
tion  radius  of  approximately  1  arc  sec  at  the  pupil.  The  lateral 
shift  at  the  conjugate  pupil  plane  between  nutators  is 


F^/FgdsinCl  arcsec) 


(  16) 
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or,  at  the  pupil  plane,  we  want 


Ax  =  d( sin  1")(F4/F3)2  <<  10  cm  (17) 

(2)  If  the  telescope  has  a  1,5  meter  primary  @  f/20  for  a 
5  mm  nutator  mirror  f3  must  be  10  cm.  For  d  -  1  cm  the  nutation 
shift,  or  pupil  smear,  will  be 

Ax  (3000/2  x  106)  cm  4.5  dt  (18) 


or  about  5%  of  a  subaperture  diameter.  For  LI  of  f/10,  diffrac¬ 
tion  smearing  is  very  small,  approximately  5  micrometers. 

This  configuration  requires  a  nutation  angular  range  of 

15  waves  *  5  mm 

or  about  1.5  millirads,  which  is  easily  achieved  with  the  current 
nutator . 

It  is  clear  that  the  simple  optical  design  of  Figure  14  will 
perform  amply  well  in  the  approximately  200  subaperture  prototype 
system,  with  negligible  degradation  due  to  subaperture  smearing. 
In  order  to  easily  match  the  detector  size  to  the  pupil  (a  prob¬ 
lem  of  magnification),  a  field  element  might  be  added  at  the  di¬ 
vider  plane,  as  in  Figure  15. 

The  later  design  was  much  larger  than  it  had  to  be  because 
the  reimaging  lens  could  be  placed  between  the  galvanometers. 

Further,  there  are  needless  design  problems  created  when  the 
system  pupil  in  a  1:1  imaging  system  is  asymmetrically  posi- 


tioned.  This  pupil  location  is  the  telecentric  position  shown  in 
Figure  15  and  all  principal  rays  are  collimated  on  the  pyramid 
side . 

All  aberrations  are  present  in  such  an  arrangement,  and  the 
design  must  correct  all  of  them,  particularly  the  worst  of¬ 
fenders,  lateral  color  and  coma. 

If  the  pupil  is  moved  to  a  symmetric  location,  there  are 
specific  gains,  since  such  systems  at  unit  magnification  are  free 
of  coma  distortion  and  chromatic  difference  of  magnification. 
Since  the  wide  optical  bandwidth  is  the  most  important  require¬ 
ment  in  a  photon  limited  environment,  and  since  it  is  difficult 
to  remove  all  chromatic  aberrations  with  the  limited  dispersions 
of  glasses  available,  obtaining  freedom  from  one  chromatic  aber¬ 
ration  is  of  utmost  value  is  securing  an  adequate  design.  In 
consequence,  the  symmetric  design  of  Figure  16  was  adapted.  The 
galvanometers  are  placed  on  either  side  of  a  simple  symmetric 
triplet.  The  actual  separation  of  the  galvanometers  and  the 
fraction  of  subaperture  smearing  that  will  exist  is  dictated  by 
the  size  of  the  coils  and  the  closeness  of  the  mirrors  to  the 
pupil  image  plane.  The  positioning  of  a  lens  between  them  will 
actually  decrease  the  optical  distance  because  of  the  generally 
finite  separation  of  principle  planes,  which,  in  the  case  of  one 
design  described  by  Conrady,  for  the  present  geometry  would  be 
6,7  mm,  reducing  the  distance  of  the  galvanometer  mirror  from  the 
pupil  image  from  10  mm  to  6.7  mm  and  reducing  the  subaperture 
smearing  from  10%  to  6.7%. 

The  actual  lens  has  not  been  designed  at  this  time,  but  the 
basic  form  is  adequate.  It  will  probably  be  necessary  to  intro- 


duce  a  third  glass  type  and  two  more  elements  to  make  an  apo- 
chromat  to  secure  the  degree  of  chromatic  correction  desired. 

The  optical  layout  that  results  is  quite  compact.  Starting 
at  the  upper  left  at  the  telescope  prime  focus,  the  beam  diverg¬ 
es  and  is  then  collimated  for  passage  through  a  correction  system 
that  includes  a  tip/tilt  corrector  and  a  deformable  mirror.  The 
beam  is  then  refocused  into  the  wavefront  sensor.  Before  it  is 
brought  to  focus,  however,  the  central  subapertures,  which  are 
occluded  by  the  telescope  secondary,  are  tagged  by  introducing  a 
beam  through  a  centrally  located  prism  from  a  modulated  subaper¬ 
ture  reference  emitter.  This  signal  will  be  used  to  align  all 
four  arrays  by  the  use  of  coils  on  the  intensif iers .  This  elec¬ 
tronic  adjustment  will  compensate  for  the  small  changes  in  the 
intensifier  which  can  be  significant  if  they  are  only  a  few  mi¬ 
crons  in  error.  The  alignment  of  the  optical  system,  to  ensure 
registration  at  the  fiber  optic  input,  will  require  mechanical 
stability  to  within  1.3  min  of  arc  (0.0004/1).  This  is  not  dif¬ 
ficult  to  obtain. 

Following  the  reference  subaperture  input,  a  field  lens 
close  to  the  correction  system  focus  will  image  the  deformable 
mirror,  which  is  at  a  pupil  image,  onto  the  principal  plan  of  the 
relay  triplet . 

Alternately,  and  preferably  except  for  cost,  the  use  of  an 
anamorphic  field  lens  to  place  the  vertical  and  horizontal  images 
on  the  respective  galvanometers  would  eliminate  subaperture 
smearing  completely.  The  respective  focal  distances  are  93.3  mm 
and  106.7  mm  for  horizontal  and  vertical  galvanometers.  The  re¬ 
quired  anamorphic  power  is  1.3  diopters,  which  is  an  easily  ob¬ 
tained  amount  of  cylinder. 

After  this  field  lens,  an  almost  completely  reflective  beam 
splitter  is  shown  when  the  reference  beam  is  introduced.  The 
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wavefront  is  reflected  because  of  its  wide  bandwidth.  The  refer¬ 
ence  wavelength  will  be  derived  from  a  modulated  laser,  Being 
monochromatic,  the  dispersion  of  the  beam  splitter  is  of  no  con¬ 
sequence.  A  corrector  plate  is  included  to  correct  the  astigma¬ 
tism.  The  residual  negative  spherical  will  be  corrected  by  an 
additional  element. 

Since  the  reference  beam  is  generated  at  a  pinhole  close  to 
the  beam  splitter,  the  uncompensatable  noncommon  path  is  small, 
and  the  stability  of  the  system  is  assured. 

The  combined  beams  proceed  to  the  galvanometer  and  the  re¬ 
lay  lens,  and  are  then  focused  on  the  quad  pyramid.  This  pyramid 
is  a  critical  element  because  its  edges  must  be  sharp.  With  the 
f/16  system  illustrated  here,  the  radius  at  an  edge  must  not  be 
greater  than  2  ym.  There  is  no  way  to  obtain  such  an  edge  by 
polishing,  but  diamond  machining  and  the  use  of  cleaved  crystals 
seem  to  provide  reasonable  means  of  manufacture. 

A  large  number  of  tetragonal  system  crystals  (but  not  dia¬ 
mond)  provide  natural  pyramid  corners.  Among  these  are  the  holo- 
symmetric  class,  which  yield  some  four  faceted  corners  and  in¬ 
clude  rutile,  anatase,  zircon  and  others.  This  last  is  commonly 
synthetically  grown  and  used  in  jewelery.  The  pyramid  class, 
yielding  the  desired  crystal,  includes  wulfenite  (PbMoO^). 

The  size  of  the  quad  is  defined  by  the  largest  object  and 
the  amount  of  nutation  to  optimize  the  SNR.  For  100  urad  ob¬ 
jects,  the  smallest  pyramid  required  is  5.12  mm.  Gems  of  tbi‘ 
size  should  be  easy  to  acquire,  particularly  zircon,  which  ac¬ 
counts  for  the  biggest  rocks  in  the  cheap  discount  stores. 

The  directions  of  the  paths  reflected  from  the  fact's  will 
depend  on  the  crystal,  so  that  the  final  design  will  require 
further  investigation  into  the  fabrication  of  pyramids.  Contin¬ 
uing  along  the  path,  the  next  element  shown  is  a  pupil  imaging 
field  lens  which  forms  an  image  on  the  input  to  the  pupil  divid¬ 
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ng  optical  system,  If  the  preceding  fields  lens  were  anamor- 
hic,  then  lenses  should  be  anamorphic  as  well,  with  the  same 
.dd-on  cylinder  of  1.3  diopters.  The  images  of  the  pupils  in  the 
livider  plane  in  this  case  will  be  1.125  mm  and  0,875  mm  because 
>f  the  power  changes  required  to  transfer  the  images  add  up. 

There  is  no  problem  in  making  the  lenslet  arrays  in  the  pupil 
divider  rectangular. 

The  last  optical  elements  are  the  lenticular  lens  arrays. 

In  this  group,  the  first  element  is  a  simple  achromat ,  which  is 
used  to  collimate  the  principle  rays.  This  ensures  that  center 
to  center  distances  between  lenslets  and  fibers  are  all  the  same. 
The  double  array  system  shown  is  used  to  ensure  that  the  pupil  is 
imaged  on  the  fiber  end,  instead  of  the  image  of  the  object  on 
the  pyramid. 

This  ensures  that  the  distribution  of  light  on  the  fiber 
does  not  move  so  that  imperfections  in  the  fiber  end  do  not 
effect  the  nutated  output.  On  the  other  hand,  the  input  angles 
into  the  fiber  will  change.  Since  these  are  always  less  than  the 
NA  of  the  fiber,  and  since  mixing  occurs,  so  that  all  angles  up 
to  the  NA  are  present  at  the  fiber  output,  this  would  seem  to  be 
the  way  to  minimize  sensitivity  to  optical  variation.  It  may, 
however,  be  overkill,  particularly  since  the  use  of  the  reference 
beam  essentially  calibrates  out  the  effects  of  imperfections  in 
each  lens,  fiber  and  detector. 

The  remaining  element,  the  detector  system,  has  been  de¬ 
scribed  in  the  preceding  section.  The  simple  intensifier  de¬ 
sign  with  the  transparent  photocathode  with  a  wide  but  manageable 
spectral  sensitivity,  0.5  to  0.8  m,  is  appropriate. 


4,1  Reflective  Systems 

The  greatest  sensitivity  is  enjoyed  by  the  intensifier 
that  uses  the  internal  reflector  and  opaque  photocathode.  The 
spectral  sensitivity,  however,  is  so  wide  that  obtaining  adequate 
chromatic  correction  in  the  optical  system  may  not  be  possible. 
Because  of  this,  reflecting  systems  have  been  examined  and  one 
design  chosen.  This  is  shown  in  Figure  17. 

It  uses  the  aplanatic  properties  of  the  Schwarzschild 
telescope  to  ensure  freedom  from  spherical  aberration,  astigma¬ 
tism  and  coma.  Since  it  uses  concentric  spherical  surfaces,  fab¬ 
rication  will  be  simple  and  inexpensive.  The  prescription  of  the 
system  is  =  1.236F,  R2  =  3.236F,  and  d  =  2F .  The  design  has  a 
focal  length  F  of  80  mm. 

The  telescope  is  used  to  collimate  the  beam,  with  a  5mm 
diameter.  The  nutator  then  reverses  the  beam  direction  and 
brings  the  return  back  along  the  opposite  side  through  the  tele¬ 
scope  again. 

There  are  a  few  innovations  which  further  simplify  the 
design.  The  input  beam  with  the  tagged  central  subaperture  is 
reflected,  aligned  and  focused  on  a  negative  field  lens  which 
will  create  an  image  of  the  pupil  between  the  galvanometers.  The 
reference  beam  will  be  introduced  through  a  10  ym  diameter  hole 
in  this  mirror.  This  pin  hole  in  the  reflective  coating  will 
filter  the  reference  beam,  absolutely  define  where  it  is.  and  en¬ 
hance,  to  a  small  degree,  the  response  of  the  system  to  small 
targets.  The  reason  for  this  is  discussed  in  Section  6  . 

There  is  one  problem  in  the  design  that  has  not  been 
resolved.  The  aplanatic  solution  is  correct  for  a  location  of 
the  stop  at  the  center  of  curvature.  The  galvanometers  could  be 
placed  at  these  locations. 

The  field  lens  mirror  that  would  be  required  would  be 
an  anamorphic  toroid  which  for  an  80  mm  focal  length  Schwarzschild 
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design  would  have  a  focal  length  of  +  72.73  mm,  to  image  the 
pupil  which  is  800  mm  upstream,  on  the  first  galvanometer  and 
-123.5  mm  to  place  the  orthogonal  image  on  the  second  galvono- 
meter . 

The  result  yields  a  5  mm  pupil  image  on  the  first 
galvo,  but  because  of  the  telescope,  creates  a  39.8  mm  pupil 
image  on  the  second  which  is  too  great.  Thus  the  use  of  galvos 
in  the  expected  stop  locations  and  proper  imaging  through  them 
is  not  possible  because  of  the  size  limitation  on  high  speed 
resonant  galvanometers.  Fortunately,  the  field  of  view  is  small 
so  that  actual  deviations  from  the  rays  shown  in  the  figure  about 
the  stop  locations  are  small  and  any  astigmatism  that  might  re¬ 
sult  from  these  small  deviations  about  the  proper  alignment  off 
axis  should  be  small. 

Residual  aberrations  can  be  balanced  and,  providing  a 
means  for  rotating  the  wavefront  and  reversing  the  sign  of  the 
astigmatic  terms  would  be  desirable.  For  this  reason,  and  to 
keep  the  galvanometers  close  together,  a  single  mirror  plus  the 
two  galvanometers  are  located  in  the  collimated  space. 

The  last  fold  mirror  out  of  the  relay  telescope  is 
mounted  on  a  piezoelectrical ly  adjustable  mount.  This  is  used  to 
lock  the  reference  beam  on  the  pyramid.  The  focusing  beam  now 
proceeds  to  the  pyramid  where  it  is  reflected  from  the  four 
faces.  The  telescope  itself  forms  an  image  of  the  pupil  115.8  mm 
from  the  apex.  Consequently,  a  simple  negative  element,  f  = 

-41.6  mm,  at  80  mm  from  the  pyramid  will  magnify  the  pupil  7.17X 
onto  the  detector  lenticular  array  257  mm  away.  A  hyperbolic 
figure  is  probably  required  but  the  size  of  the  field  on  the  mir¬ 
ror  of  2 . 5  mm  is  so  small  that  a  spherical  surface  may  be  suffi¬ 
cient.  At  this  point  in  the  optical  path,  the  optics  serve  only 
to  get  the  signal  to  the  right  detector  and  performance  is  not 


critical.  The  arrangement  of  elements  in  the  layout  is  compli¬ 
cated  only  because  of  the  need  for  the  reflected  beams  to  miss 
occluding  structures. 

The  design  requires  a  minimum  number  of  elements,  two  good 
spherical  mirrors,  5  negative  spherical  mirrors,  and  a  few  fold 
mirrors.  It  is  completely  free  of  dispersive  elements.  Even 
the  beam  splitter  and  its  loss  has  been  replaced  by  a  signal  en¬ 
hancing  hole.  The  result  is  an  optical  system  that  is  very  ef¬ 
ficient.  If  yttria  protected  silver  mirrors  are  employed,  the 
mean  reflectance  from  0.44  to  0.84  pm  will  be  0.984.  With  twelve 
reflections,  the  transmission  will  be  82.4%  from  input  to  the  de¬ 
tector  . 
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5.0  SIGNAL  PROCESSING  ELECTRONICS 

_  g 

The  Large  Scale  I  Sensor  (LSI  )  electronic  processor  de¬ 
sign  is  intended  to  minimize  operator  intervention.  Novel  fea- 
3 

tures  of  the  LSI  are  exploited  to  permit  automatic  system  cal¬ 
ibration,  extensive  self  test,  and  automatic  dynamic  alignment 
during  system  operation.  These  extraordinary  capabilities  are 
in  addition  to  a  signal  processor  architecture  which  totally 
eliminates  subaperture  level  electronic  components  and  adjust¬ 
ment.  This  proprietary  architecture  results  in  a  considerable 
reduction  in  circuitry  compared  to  the  traditional  processor 
per  subaperture  approach. 

5. 1  Processor  Architecture 

3 

The  architecture  of  the  LSI  electronic  signal  proces¬ 
sor  is  based  on  the  fact  that  high  speed  digital  memory  is  the 
simplest,  most  reliable,  and  least  expensive  electronic  compo¬ 
nent  to  use  for  the  design,  fabrication,  operation  and  mainte¬ 
nance  of  a  large  scale  signal  processor.  Digital  memory  coupled 
with  high  speed  digital  arithmetic  elements  permits  Time  Domain 
Multiplexing  (TDM)  and  this  sharing  of  processing  resources  min¬ 
imizes  the  circuitry  required  for  signal  processing.  An  addi¬ 
tional  positive  feature  is  the  inherent  flexibility  of  program¬ 
mable  digital  arithmetic  elements, 

3 

There  are  three  major  subdivisions  to  the  LSI  elec¬ 
tronic  processor.  These  are  the  Front  End  Processor  (FEP),  the 
Programmable  Algorithm  Processor  (PAP),  and  the  System  Super¬ 
visor  Processor  (SSP).  Figure  18  illustrates  the  data  flow 
among  the  processors. 

5.2  Front  End  Processor 

The  FEP  digitizes,  normalizes  and  buffers  the  detect- 


ed  signals.  These  functions  are  performed  serially  at  the  de¬ 
tector  readout  rate  (10  MHz  for  256  subapertures).  Figure  19 
illustrates  the  pipeline  approach  to  the  FEP. 

The  processor  architecture  was  designed  so  that  any 
of  the  large  scale  (up  to  256  subaperture)  detectors  selected 
for  use  in  a  system  as  described  in  the  next  section,  can  be 
accommodated.  The  only  changes  will  be  in  the  detector 
preampl i f iers . 

Each  detector  array  requires  a  charge  sensitive  pre¬ 
amplifier.  For  the  Reticon  RL256C  Solid  State  Line  Scanner 
(SSLS)  the  maximum  charge  on  any  photodiode  is  approximately  4 
picocoulombs .  With  10  MHz  readout,  the  resulting  current  is  at 
most  40  nanoamps.  The  preamplifier  is  mounted  in  close  proximi¬ 
ty  to  the  diode  array  and  connected  via  carefully  shielded 
cables  to  eliminate  noise  currents  due  to  stray  magnetic  fields. 

The  SSLS  preamplifier  block  diagram  is  given  in  Figure 
20,  and  Figure  21  shows  the  timing  sequence  for  charge 
readout.  Each  photodiode  is  precharged  to  +5V  and  incident 
electrons  deplete  this  charge.  The  output  signal  is  the  current 
required  to  recharge  the  diode.  Hence,  the  current  to  voltage 
converter  output  is  negative  spikes,  as  shown  by  VI  of  Figure 
21  . 

These  spikes  are  integrated  and  buffered  by  the  sample 
;:d  hold  (S/H).  The  two  SSLS  outputs  (even  numbered  photodiodes 
on  one,  odd  numbered  photodiodes  on  the  other)  are  combined  at 
the  analog  multiplexer  and  sent  to  the  analog  to  digital  con¬ 
verter  (A/D). 

The  A/D  is  an  integrated  circuit,  8  bit  unit  of  com¬ 
mercial  manufacture  and  capable  of  30  MHz  conversion.  The  10 
MHz  readout  rate  is  well  within  the  A/D  capability.  After  this 
digitization  the  detected  signals  are  transmitted  via  80  MHz 
fiber  optic  cable  to  the  normalization  and  buffer  sections  of 
the  FEP .  Thus,  a  minimum  of  electronic  circuitry  is  mounted 
with  the  optical  package. 
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This  minimum  circuitry  consists  of  a  clock  driver 
circuit,  a  +5V  regulator,  two  preamplifiers,  two  integrators, 
two  S/K,  one  A/D,  and  one  fiber  optic  transmitter  per  detector 
array.  The  total  is  approximately  50  integrated  circuits  which, 
with  cables,  weigh  less  than  a  Kilogram. 

The  remainder  of  the  FEP  is  installed  in  the  same  card 
cage  as  the  PAP.  This  remainder  consists  of  2560  bytes  of  mem¬ 
ory  and  two  arithmetic  elements  per  detector  array.  The  memory 
is  organized  as  two  segments  of  256  bytes  each  (one  byte  per 
subaperture)  and  two  buffers  of  1024  bytes  each  (four  bytes  per 
subaperture) . 

The  first  256  byte  memory  is  connected  to  an  8  bit 
digital  adder.  As  each  digitized  photodiode  output  is  acquired 
ry  the  fiber  optic  receiver,  it  is  latched  and  the  dark  signal 
contribution  is  removed.  The  dark  signal  for  each  photodiode 
is  stored  in  a  separate  memory  location  of  the  256  byte  RAM. 
These  RAM  constants  are  determined  and  loaded  by  the  SSP  (see 
Section  5.4). 

After  dark  signal  subtraction,  the  new  8  bit  word  is 
latched  and  the  photodiode  output  is  corrected  for  nonuniform 
responsivity  of  the  detector  system.  This  correction  is  im¬ 
plemented  as  an  8  bit  multiplication  using  a  separate  normaliz¬ 
ing  factor  for  each  photodiode.  The  normalizing  factors  are 
stored  in  a  256  byte  RAM  and  again,  the  RAM  constants  are  de¬ 
termined  and  loaded  by  the  SSP  (see  Section  5.4). 

The  8  bit  output  of  the  multiplier  circuit  is  stored 
in  one  of  the  two  1024  byte  buffer  memories;  the  other  1024 
byte  memory  is  supplying  data  from  the  previous  nutation  cycle 
to  the  PAP.  The  1024  bytes  are  sufficient  for  4  samples  per 
nutation  cycle  per  photodiode  (256  photodiodes).  The  data  en¬ 
tering  the  buffer  memory  is 
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where  P(i,k)  is  the  kth  time  sample  (1  <_  k  £  4)  of  the  ith  pho¬ 
todiode,  D(i)  is  the  ith  dark  signal,  value  and  R(i)  is  the  ith 
responsivity  correction  constant.  Thus,  the  buffer  is  filled 
with  the  output  data  of  the  FEP . 

The  FEP  is  a  pipeline  processor  with  6  stages  to  the 
pipeline.  The  first  stage  is  the  charge  sensitive  amplifier  and 
the  integrator.  The  second  stage  is  the  S/H  unit  and  the  A/D. 
The  third  stage  is  the  fiber  optic  transmitter-receiver.  The 
fourth  stage  is  the  dark  current  removal  addition  circuit.  The 
fifth  stage  is  the  responsivity  normalizing  multiplication  cir¬ 
cuit.  The  sixth  and  last  stage  is  the  buffer  memory.  Data 
moves  from  one  stage  to  the  next  at  100  nS  intervals  and  hence 
each  photodiode  output  is  amplified,  digitized,  and  normalized 
600  ns  after  readout. 

The  FEP  pipeline  is  four  bytes  wide,  one  byte  for 
each  detector  array.  All  four  data  bytes  are  processed  simul¬ 
taneously  and  synchronously;  thus  a  single  RAM  address  generator 
and  a  single  SSLS  clock  generator  are  adequate  for  all  four  de¬ 
tector  arrays.  The  digital  portion  of  the  FEP  (i.e.,  that  por¬ 
tion  of  the  FEP  not  mounted  in  the  optics  package)  requires  ap¬ 
proximately  150  integrated  circuits,  all  operating  at  a  10  MHz 
clock  rate.  All  circuits  are  commercially  available. 


5.3  Programmable  Algorithm  Processor  (PAP) 

The  digitized  photodiode  signals  roust  be  processed  in 
order  to  extract  the  subaperture  tip/tilt  information.  There 
exist  various  algorithms  for  extracting  this  information  (see  Sec¬ 
tion  5.3.2) f  but  all  algorithms  are  similar  in  that  they  involve 
linear  combinations  of  the  16  detected  signals  (four  detectors/ 
subaperture  and  four  time  samples/detector)followed  by  two  ratio 
computations.  The  ratio  computation  (division)  normalizes  the 
tip/tilt  measurement  and  hence  the  final  tip/tilt  value  is  inde¬ 
pendent  of  optical  power.  Since  all  of  the  algorithms  possess 
similar  computation  requirements,  the  Synchronous  Centroid  Sam¬ 
pling  (SCS)  algorithm  is  used  as  the  processor  example.  Figure 
22  illustrates  the  computation  requirements  of  the  SCS  algo¬ 
rithm. 

Two  approaches  to  algorithm  implementation  were  con¬ 
sidered.  The  first  to  be  presented,  full  pipeline  processing, 
is  straightforward  but  relatively  inflexible.  The  preferred  ap¬ 
proach,  multiple  microprogrammed  processors,  possesses  more  com¬ 
plex  timing  but  this  is  more  than  offset  by  the  inherent  flexi¬ 
bility  of  microprogramming.  Either  approach  is  practical  with 
commercially  available  circuits. 

5.3.1  Pipeline  processor.  A  fully  pipelined  proces¬ 
sor  requires  a  separate  pipeline  stage  for  each  arithmetic  op¬ 
eration.  A  pipeline  stage  consists  of  a  primitive  arithmetic 
function  and  a  data  register  (latch).  Figure  23  illustrates  the 
pipeline  implementation  of  the  SCS  algorithm  computation. 

The  SCS  algorithm  requires  three  primitive  arith¬ 
metic  functions.  These  are  complement  (output  =  -input),  sum 
(output  =  inputl  +  input2),  and  a  divide  primitive  (shift/sub¬ 
tract  circuit).  Each  latch  acquires  a  new  data  word  at  100  ns 
(10  MHz  clock)  intervals  and  all  data  latches  are  clocked  simul- 
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Q1  ARRAY 


(n)  is  nth  (1  <  n  <  4)  sample  period  during  a  nutation  cycle 
Q  (n)  is  n1"*1  sample  from  mth  (1  <  m  <  4)  detector  array. 


Let  SI  =  Q1 ( 2 )  -  Q4( 2 ) 

52  =  Q4(4 )  -  01(4) 
T1  =  SI  +  S2 

53  =  02(2)  -  Q3( 2 ) 

54  =  03(4)  -  Q2(4) 
T2  =  S3  +  S4 

XN  =  T1  -  T2 
XD  =  T1  +  T2 

Therefore 

X  -  XN  i  XD 


1st  level 
1st  level 
2nd  level 
1st  level 
1st  level 
2nd  level 

3rd  level ,  x  numerator 
3rd  level,  x  denominator 


x  tilt 


Similarly 

55  =  Q 1 ( 1 )  -  Q2( 1 ) 

56  =  Q2( 3 )  -  Q1 ( 3) 

T3  =  S5  +  S6 

57  =  Q4( 1 )  -  Q3( 1 ) 

58  =  Q3( 3 )  -  Q4( 3) 

T4  =  S7  +  S8 

YN  =  T3  -  T4 
YD  =  T3  +  T4 
Y  =  YN/YD 

Figure  22.  SCS  algorithm  computat  on. 
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taneously.  Thus,  the  data  for  each  computational  step  ripples 
through  the  processor  as  if  the  data  were  water  flowing  in  a 
pipe.  The  time  required  for  data  to  pass  from  input  to  output 
is  100  ns  x  number  of  pipeline  stages,  which,  in  this  case,  is 
13.  The  first  subaperture  data  is  available  1300  nS  after  the 
nutation  cycle,  the  second  aperture  data  is  available  1700  ns 
after  the  nutation  cycle,  the  third  is  available  2100  n.s  after 
the  nutation  cycle,  etc. 

After  pipeline  processing,  the  data  is  stored  in 
one  of  two  512  byte  buffer  memories;  the  other  buffer  memory  is 
supplying  subaperture  tip/tilt  data  to  the  deformable  mirror 
control  circuitry.  These  two  buffers  interchange  their  roles  at 
the  end  of  each  nutation  cycle.  The  512  bytes  of  memory  are  or¬ 
ganized  as  256  bytes  of  subaperture  x  tilts  and  256  bytes  of 
subaperture  y  tilts. 

This  pipeline  processing  technique  is  the  most 
straightforward  approach  to  algorithm  implementation.  However, 
the  algorithm  is  defined  by  the  interconnections  of  the  primi¬ 
tive  arithmetic  elements,  i.e.,  by  the  wiring.  Even  minor 
changes  to  the  algorithm  would  require  hardware  revision  and 
possibly  refabrication.  For  instance,  spatial  or  temporal  av¬ 
eraging  of  the  subaperture  data  would  require  extensive  modifi¬ 
cations  and  additions  to  the  pipeline  processor.  Furthermore, 
any  averaging  implemented  by  hardwired  pipeline  processing  would 
involve  fixed  filter  characteristics.  Thus,  a  hardwired  filter 
could  not  be  changed  to  adapt  for  differing  signal  to  noise  ra¬ 
tios  or  changing  seeing  conditions. 


HO 


5,3,2  Microprogrammed  processor.  Flexibility  of  al- 

3 

gorithm  definition  is  a  desirable  feature  for  the  LSI  subaper¬ 
ture  signal  processor.  Commercial  bit-slice  circuitry  has  ad¬ 
vanced  to  the  stage  where  microprogrammed  arithmetic  elements 
provide  a  practical  mechanism  for  programmable  algorithm  defini¬ 
tion  . 

The  Programmable  Algorithm  Processor  (PAP)  is  de¬ 
signed  to  implement  any  subaperture  signal  processing  algorithm 
for  which  the  output  signals  are  a  linear  combination  of  the  in¬ 
put  signals. 

The  heart  of  the  PAP  is  the  Am29203  four-bit  ex¬ 
pandable  microprocessor  slice  (See  Figure  24).  This  device  is 
a  48  pin  integrated  circuit  containing  16  storage  locations  and 
an  Arithmetic  Logic  Unit  (ALU).  The  storage  locations  and  ALU 
are  equivalent  to  the  latches  and  primitive  arithmetic  functions 
of  the  pipelined  processor.  However,  the  PAP  algorithm  is  de¬ 
fined  by  a  sequence  of  32  bit  control  signals  where  the  pipe¬ 
lined  algorithm  was  defined  by  wires.  The  32  bit  control  sig¬ 
nals  (called  the  microcode)  are  stored  in  a  RAM  and  hence  the 
PAP  algorithm  can  be  modified  by  changing  the  contents  of  this 
RAM.  Figure  25  illustrates  the  block  diagram  of  a  micropro¬ 
grammed  PAP  element. 

To  illustrate  the  flexibility  of  microprogram¬ 
ming  approach,  consider  the  implementation  of  the  SCS  algo¬ 
rithm.  Figure  26  illustrates  the  microcode  required  to  dupli¬ 
cate  the  pipeline  processor  of  Figure  23.  Each  line  of  micro¬ 
code  requires  100  ns  for  execution.  The  entire  SCS  algorithm 
(less  the  division)  requires  16  operations  or  1600  ns  .  A  10 
kHz  nutation,  400  subaperture  system  requires  output  at  400  ns 
intervals,  so  a  full  256  subaperture  system  would  require  lour 
arithmetic  elements  as  shown  in  Figure  25. 
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UT  DATA  BUS  FROM  BUFFER  MEMORIES 


Input  R5  =  Ql(l)  -  Q2(l)  S5 

R7  =  Q4( 1 )  -  Q3( 1 )  S7 

R1  =  Q1 ( 2 )  -  Q4( 2)  SI 

R3  =  Q2( 2 )  -  03(2)  S3 

R6  =  02(3)  -  Ql(3)  S6 

R8  =  Q3( 3 )  -  Q4 ( 3 )  SB 

R2  =  04(4)  -  Q1 (4)  S2 

R4  =  Q3(  4  )  -  Q2(  4 )  S4 

(Release  input  bus) 

R9  =  R 1  +  R2  T1 

RIO  =  R3  +  R4  T2 

RU  =  R9  +  RIO  XN 

R 1 2  =  R9  -  RIO  XD 

R 1 3  =  R5  +  RC>  T3 

R 14  -  R7  +  R8  T4 

R 1 5  --  R 1 3  -  HI  !  YN 

RIO  =  R 1 3  +  R 1 4  YD 


Figure  20.  PAP  mirrococlr  for  S('S  algorithm. 
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Now  consider  an  algorithm  change,  specifically  a 
derotation  as  described  in  Section  6.  The  pipeline  processor 
would  require  refabrication  in  order  to  implement  this  function. 
Two  new  primitive  arithmetic  elements  (a  multiplier  circuit 
and  a  derotation  memory)  would  be  required.  The  microprogrammed 
processor  would  require  only  a  change  of  microcode,  i.e.,  a 
softward  change.  Figure  27  shows  the  new  microcode  required  to 
implement  the  SCS  algorithm  plus  derotation.  This  new  microcode 
requires  2400  ns,  so  six  arithmetic  elements  as  shown  in  Figure 
25  would  be  required. 

The  PAP  uses  multiple  identical  programmable 
processor  elements  in  place  of  arithmetic  hardware  dedicated  to 
each  subtask  of  the  algorithm  processing.  The  inherent  flexi¬ 
bility  of  programmable  processors  is  easily  exploited  because 
new  microcode  can  be  loaded  in  a  single  nutation  cycle.  The 
programmable  processors  can  employ  temporal  or  spatial  filter 
functions,  and  the  parameters  of  these  filters  can  be  changed 
easily  to  accommodate  the  signal  characteristics  (SNR,  seeing, 
turbulence  bandwidth,  etc.).  The  filters  could  even  be  changed 
instantly  during  a  data  run  simply  by  switching  to  a  different 
block  in  the  microcode  RAM. 

Figure  28  illustrates  the  PAP  architecture.  The 
FEP  buffer  memories  connect  to  a  32  bit  input  data  bus.  The 
processor  elements  access  this  bus  sequentially,  i.e.,  PI  ac¬ 
quires  the  data  from  subaperture  1.  While  PI  processes  this 
data,  P2  acquires  the  data  from  subaperture  2,  etc.  All  proces¬ 
sors  execute'  the  same  microcode,  but  each  processor  is  out  of 
phase  with  its  neighbors.  This  arrangement  is  effectively  a 
pipelined  instruction  (not  pipelined  data)  architecture  lie- 
cause  any  given  microcode  step  is  executed  first  by  1‘1,  then  by 
P2 ,  then  P3 ,  etc.  By  organizing  the  processors  in  this  manner. 


Qm(n)  is  n1'^1  time  sample  from  m^*1  detector  array. 
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Figure  27. 


there  is  never  a  condition  where  two  processors  simultaneously 
require  the  same  shared  resource  (input  bus,  output  bus,  etc.). 
Thus,  each  processor  element  is  invisible  to  the  other  processor 
elements . 

The  number  of  microcoded  processor  elements  re¬ 
quired  is  dependent  on  the  length  of  the  microcode,  i.e.,  the 
complexity  of  the  algorithm.  For  256  subapertures  and  10  kHz 
nutation,  the  requirement  is  one  processor  element  per  four  mi¬ 
crocode  steps.  Eight  processor  elements  identical  to  the  illus¬ 
tration  of  Figure  25  would  be  sufficient  for  32  microcode  steps 
per  subaperture.  Programmability  makes  each  element  extremely 
powerful  and  32  steps  are  more  than  adequate  for  all  algorithms 
under  consideration. 

5.4  System  Supervisor  Processor  (SSP) 

The  SSP  supplies  general  purpose  intelligence  to  the 

3 

LSI  Sensor.  The  heart  of  the  SSP  is  a  commercial  microproces- 

3 

sor ,  the  microNova  (see  Figure  29).  The  SSP  controls  the  LSI' 
processor  initialization,  self  test,  calibration,  nutator  mir¬ 
rors,  dynamic  subaperture  alignment,  and  autolock. 

3 

When  the  LSI  Sensor  is  turned  on.  the  SSP  will  down¬ 
load  the  microcode  to  the  PAP  and  will  download  the  calibration 
constants  to  the  FEP .  This  control  information  will  be  used  to 
execute  a  self  test  of  the  processor  functions.  The  full  aper¬ 
ture  reference  LED  will  be  pulses  and  all  detector  system  ele¬ 
ments  (thermoelectric  cooler,  magnetic  focus,  solid  state'  line 
scanner,  preamplifiers,  and  A/D  converters)  will  be  verified  op- 
rat  ion  al  . 

The  SSP  will  then  verify  the  PAP  functionality  by  ap¬ 
plying  arbitrary  subaperture  signals  to  the  input  bus  and  exam¬ 
ining  the  PAP  outputs.  This  test  of  the  FEP  and  PAP  will  be 


s,s 


r 


automatic  and  will  require  approximately  one  minute. 

If  desired,  a  detector  calibration  can  be  initiated. 
This  calibration  will  use  the  internal  reference  LED  as  an  op¬ 
tical  source  to  calculate  dark  signal  and  responsivity  for  each 
cell  of  the  four  detector  arrays.  The  new  calibration  constants 
will  then  be  sent  to  the  FEP .  A  full  calibration  will  require 

approximately  five  minutes. 

3 

The  LSI  Sensor  will  then  begin  supplying  subaperture 

tilt  data.  While  the  FEP  and  PAP  are  totally  occupied  in  this 

task,  the  SSP  is  free  to  monitor  the  subaperture  reference  LED 

signals  and  will  maintain  the  subaperture  alignment  using  the 

magnetic  deflection  elements  in  the  detector  arrays.  The  SSP 

will  also  monitor  the  full  aperture  reference  LED  signals  and 
3 

maintain  the  LSI  optical  axis  alignment  using  the  autolock  mir¬ 
ror.  The  subaperture  alignment  will  not  intrude  on  the  wave- 
front  sensing  computations;  however,  the  optical  axis  alignment 
will  involve  loss  of  data  for  a  full  nutation  cycle.  When  the 
full  aperture  reference  LED  is  pulsed  on,  the  input  optical  sig¬ 
nal  is  lost.  Thus,  a  nutation  cycle  is  stolen  for  the  dynamic 
alignment  function.  A  loss  of  one  cycle  in  32  would  be  adequate 
for  100  Hz  dynamic  alignment  control. 

The  SSP  will  also  monitor  and  control  the  nutation 


mirrors.  The  nutation  radius  can  be  varied  to  optimize  the  sig¬ 
nal  to  noise  ratio.  It  is  even  possible  to  use  different  am¬ 
plitudes  on  the  x  and  y  nutators  for  independent  optimization. 
Since  the  SSP  controls  the  PAP  microcode,  any  tip/tilt  scale 
changes  to  dynamic  nutation  control  are  easily  and  automatically 
removed  by  the  mierocoded  processor  elements. 

The  exact  definition  of  the  SSP  tasks  will  change  dur¬ 
ing  the  LSI'*  fabrication  and  installation.  Nutation  and  optical 
alignment,  parameters  are  very  dependent  on  the  telescope  opt  i  cs . 


> 


The  bandwidth  required  for  dynamic  subaperture  and  full  aper¬ 
ture  alignment  is  dependent  on  the  telescope  mechanics.  The 
self  test  and  calibration  procedures  may  require  expansion  to 
include  other  elements  of  the  overall  optical  system.  Since 
the  SSP  functions  are  defined  by  micronova  programs,  the  SSP 
control  parameters  can  be  easily  modified  in  software  without 

hardware  refabrication.  This  designed-in  flexibility  is  typical 
3 

of  the  LSI  signal  processing  system. 


91 


6 . 0  DETECTOR  STUDY 

The  four  detector  planes  in  the  LSI  sensor  are  not  limited 

to  a  specific  geometric  configuration  of  detectors,  as  are  the 

3 

quad  cell  arrays  for  the  I  sensor.  Therefore,  a  wide  choice  of 
possible  off-the-shelf  or  special  devices  may  be  considered  for 
the  detector  arrays.  The  detectors  must  meet  compensated 
imaging  requirements  and  provide  a  high  speed  serial  output  in 
order  to  achieve  a  configuration  with  a  minimum  number  of  com¬ 
ponents.  High  quantum  efficiency  and  a  photon  limited  signal- 
to-noise  ratio  can  be  obtained,  in  general,  only  by  the  use  of 
photomultipliers  or  similar  devices.  The  high  speed  serial 
output  sensors  that  are  available  and  which  meet  the  minimum 
hardware  requirements  are  the  integrating  sensor  arrays  such  as 
the  CCD's  and  self-scanned  diode  arrays. 

The  marriage  of  silicon  arrays  and  intensifying  devices 
would  meet  the  objectives  of  the  program.  This  is  the  subject 
of  this  study. 

The  variety  of  methods  of  intensification  and  transferring 
the  resulting  charge  that  is  generated  to  a  detector  can  be 
categorized  by  the  method  of  coupling  the  amplifier  to  the 
detector.  These  categories,  including  direct  (nonintensif ied) 
photon  detection,  are 

(1)  direct  mode:  nonintensif ied  direct  photon  detection 
(silicon  arrays); 

(2)  electron  bombardment  (EB)  mode:  the  photo-electron 
is  accelerated  20  kV  to  30  kV  and  strikes  silicon,  directly 
creating  electron  hole  pairs; 

(3)  photon  coupled  mode:  the  photo-electrons  are  amplified 
in  an  intensifier  and  strike  a  phosphor.  The  photons  emitted 
are  then  fiber  optically  coupled  to  the  array. 

(4)  Electron  coupled  mode:  the  detector  array  is  designed 
as  an  array  of  anodes  to  collect  the  current  generated  by  an 


intensif ier . 

(5)  Discrete  element  arrays:  this  stands  for  photomulti¬ 
pliers  and  represents  a  fall  back  position  that  can  be  used  for 
small  arrays  for  breadboard  purposes. 

These  categories,  which  will  be  discussed  in  detail  in 
the  following  sections,  are  illustrated  in  Figures  30  through  34. 
The  complications  are  a  consequence  of  applying  means  to  pre¬ 
serve  all  of  the  photoelectrons  that  will  be  generated  at 
a  cathode,  or  to  permit  the  use  of  existing  commerically  avail¬ 
able  devices.  The  conclusion  of  the  following  study  was  that 
the  configurations  of  Figures  31b  and  c  were  optimum  for  the 

3 

LSI  system,  and  that  photomultipliers.  Figure  34, would  provide 
a  cost  effective  solution  for  small  breadboard  systems. 

The  quantities  required  for  evaluating  the  various  methods 
are  simple. 

(1)  Detection  quantum  efficiency  (QE  =  n):  sensitivity 
of  cathode  or  fraction  of  photons  yielding  a  photoelectron. 

(2)  Collection  quantum  efficiency  (hc>:  fraction  of 
photoelectron  pulses  that  get  to  detector. 

(3)  Noise  factor:  a  consequence  of  variation  of  gain  in 
intensifying  devices. 

(4)  Detection  noise  (ne). 

The  phase  variance  of  all  of  the  wavefront  sensors  including 
those  not  studied  in  this  program  are  given  by 

2  _  f (object , sensor)  /Oo\ 

0  A  ^  r>  V*-'-/ 

*  (SNR)2 

2 

In  the  photon  counting  limit ,  (SNR)  =  number  of  detected 
photons,  N^,  which  is  given  by  *  nnc/F.  Thus,  the  phase 

variance  is  proportional  to 
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98 


where  ne  is  the  number  of  noise  electrons  and  g  is  the  gain. 

Going  back  to  the  five  sensing  systems,  we  find  a  basic 
difference  in  the  sensing  employed.  In  the  first  category, 
almost  all  photosensitive  methods  are  applicable.  Basically, 
however,  because  of  the  devices,  silicon  is  the  main  contender. 

In  all  of  the  other  categories,  the  same  photocathodes  are 
required.  Thus,  for  comparing  the  last  four  categories,  a 
simple  value  of  n  will  be  sufficient  and  it  could  be  ignored  or 
simply  set  to  1. 

Since  n  is  a  function  of  wavelength,  the  evaluation  of 
detectors  in  case  (1)  is  complicated  by  the  fact  that  we  must 
integrate  over  the  bandwidth  and  consider  the  effect  of 
wavelength  on  the  measurements. 

P.  i.  Direct  Mode  Systems 

We  are  concerned  here  with  the  efficiency  of  simple 
silicon  based  sensors  and  will  analyze  CCD's,  the  CID  and  the 
self-scanned  diode  array.  The  possible  value  of  these  devices 
over  current  photocathodes  can  be  seen  by  comparing  the  number 
of  photoelectrons  that  will  be  generated.  The  detector  sensiti¬ 
vities  and  the  solar  flux  as  seen  through  the  atmosphere  is 
shown  in  Figure  3b.  The  total  photoelectron  count  for  the 
different  sensors  integrated  over  the  spectrum  is  decidedly 
better  for  the  silicon  array.  The  result  of  the  integration  is 
given  in  Table  4. 

The  plain  silicon  devices  are  seen  to  be  five  to  ten 
times  more  productive  than  the  cathodes  that  provide  photoelec¬ 
trons  for  intensification. 
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SENSITIVITY  (A/W) 


Figure  35.  Solar  spectrun  and  detector  sensivities. 


TABLE  4.  INTEGRATION  RESULT,  TOTAL  PHOTO-ELECTRON  COUNT. 


The  accuracy  of  tilt  measurements  in  a  quantum  limited 
systems  is  dependent  on  object  image  size  and  on  the  resolution 
of  the  optics.  If  the  object  image  width,  including  diffraction, 
is  approximated  by 


a2.  =  ot2  ,  .  +  (X/D)2  f  p21 

image  obj  ' 

where  D  is  the  full  telescope  aperture,  and  the  power  received 
in  a  measurement  with  a  subaperture  L  is  given  by 

Boc  ,  . 2  L2  (23> 

obj 

where  B  is  the  brightness,  so  that  the  number  of  photons  is 

Ba2  L2  n /hv  (24) 


the  variance  on  a  tilt  measurement  is  proportional  to 


V  “  {%bj2  +  <X/D>2} 


(25) 


-  (hv/ABL2 ) { 1  +  ( A /D ) ( 1 / aobj ) ) 


We  find  that,  if  the  quantum  efficiency  is  constant  for  all 
targets  of  interest,  >  A/D,  the  accuracy  will  improve 

with  wavelength.  If  the  sensitivity  is  constant,  then  the 
number  of  photoelectrons  is  Ng  =  S-Pd  and  it  is  important 
to  measure  over  as  wide  a  spectral  band  as  possible. 

The  tilt  variance  induced  by  the  atmosphere  is  given 

by 

(1/2)  <■«.«>=  <ax2>  =  <ay2>  =  0.975(F/k)2  {D^  (  2a ) }  /  (  2a ) 2 

(26) 
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where  the  structure  function  for  the  aperture  is  D  (2a)  and  is 
given  by 

c 

2 . 9k2  r5  ^ 3  •  f  C  ( z )dz .  (27) 

J  nv 

o 

The  function  is  independent  of  wavelength,  thus  it 
does  not  matter  what  wavelength  is  used  for  a  wavefront  tilt 
measurement.  In  a  shearing  interferometer,  where  optical 
phase  is  measured  instead  of  tilt,  long  wavelength  measurements 
are  less  effective  than  short  because  the  tilt  is  given  by 
■1>*X/L  so  that  the  tilt  variance  is 

aa2  =  (X/L)2o02  -  (X/L)2(l/N)  (28) 

The  result  of  this  is  to  reduce  the  ef fecti veness  of  some  of 
the  wider  sensitivity  cathodes.  This  is  shown  in  the  second  row 
of  Table  4. 

If  tilt  and  phase  measuring  sensors  were  to  use  the 
same  cathode,  the  other  considerations  described  in  Section  3 
on  performance  would  apply.  The  increased  capability  of  a  tilt 
sensor  to  utilize  wide  optical  bandwidths  effectively  will 
produce  an  improvement  in  tilt  variance  by  factors  of  3  to  5 
if  opaque  cathodes  can  be  used. 

Another  possibility,  if  a  suitable  geometry  can  be 
devised,  is  to  use  opaque  cathodes  such  as  GaAsP  (Varian's  S- 
type  cathode)  which  has  a  quantum  efficiency  of  40"  at  450  y m 
(theoretical  maximum  for  a  cathode  is  50%).  The  relative  number 
of  photoelectrons  and  the  improvement  in  variance  are  given  in 
the  last  column  of  Table  4.  The  improvement  in  performance  is 
sufficiently  great  to  warrant  consideration  of  a  special  tube 
design . 
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The  silicon  devices  are  still  a  great  deal  more 
effective.  The  SNR  for  CCD's  can  be  quite  large  for  moderate 
signals  because  of  the  small  number  of  rms  noise  electrons 
that  can  be  expected.  Some  CCD's  have  as  few  as  30  noise 
electrons/bucket.  The  reticon  arrays  can  be  held  to  800. 

The  phase  variance  for  a  nonintensi f ied  diode  is 


when  is  the  number  of  photoelectrons  seen  by  the  diode. 


ne  are  the  noise  electrons. 


For  the  intensified  diode 


cathode 

when  N  .  is  the  number  of  cathode  photoelectrons.  Using  the 
ratios  Si/bi-alkali  of  11.5  and  5.28,  obtained  from  Table  4, 
we  can  find  the  minimum  count  and  rate  below  which  intensifica¬ 
tion  is  required  for  various  readout  noise  levels.  These  are 
given  in  Table  5. 

Since  the  flux  for  Cl  objects  and  similar  cathodes 
4  7 

provides  from  10  counts/sec  to  5  x  10  counts/sec,  we  find 
that  intensification  will  be  required. 

This  is  not  an  absolute  conclusion,  however.  At  the 
low  noise  levels,  the  CCD  is  very  competitive.  From  10  to  30 
noise  electrons  per  sample  (NES)  have  been  reported  in  buried 
channel  CCD  shift  registers.  In  the  Fairchild  256  linear  array, 
noise  levels  of  150  electrons  per  charge  packet  have  been 
obtained,  and,  at  low  temperatures  (160  K)  at  telephone  rate 
(15  kHz),  the  NES  was  about  25. 

Consequently,  although  the  present  conclusion  indicates 


TABLE  5 

Bi-alkali  cathode  counts  and  rates  to  be  exceeded  if 
nonintensif ied  arrays  are  to  be  superior  (25  ys 
integration  time). 


If  Si/BiA  ratio  is 


11.6 

Count  Rate 

Count 

5.28 

Rate 

Best  CCD 

5 

0.21 

8  kHz 

1.1 

44  kHz 

10 

0.83 

33  kHz 

4.4 

177  kHz 

30 

7.45 

298  kHz 

39.8 

1.59  MHz 

100 

82.8 

3.3  MHz 

443 

17.7  MHz 

Ret  icon 

300 

745 

30  MHz 

3,980 

159  MHz 

1,000 

8,282 

330  MHz 
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a  need  for  intensification,  the  matter  is  not 
and  the  future  development  of  CCD's  should  be 
The  comments  on  CCD  geometry  in  the  following 
relevant  in  choosing  a  detection  technique  as 
sensitivity . 


entirely  closed, 
watched  carefully, 
section  are  as 
the  question  of 


There  is  one  further  consideration;  cathodes  are 
improving  and  it  may  be  possible  to  obtain  50%  quantum 
efficiencies  in  the  future  as  the  development  of  III-V  cathode 
materials  proceeds.  If  this  happens,  cathodes  with  intensifica¬ 
tion  will  be  superior. 


6.2  Geometries  and  Processing  Requirement 

Before  discussing  the  coupling  techniques  for  the 
EB  mode  or  the  photon-coupled  mode,  a  discussion  of  the  geomet¬ 
ries  of  the  arrays  and  the  processing  required  is  desirable  in 
order  to  define  how  to  use  them. 

The  specifications  that  apply  to  the  devices  are 
simple.  Each  array  must  be  read  out  every  25  psecs  in  order 
to  provide  a  complete  data  set  at  a  10  kHz  rate.  The  number 
of  elements  in  the  array  then  are  limited  by  the  speed  of  the 
device. 

High  speed  devices  provide  10  MHz  video  (and  some 
as  high  as  25  MHz),  consequently,  256  element  arrays  can  be 
used . 

Four  different  kinds  of  arrays  are  to  be  considered. 
These  are  (1)  buried  channel  charge  coupled  devices,  BC-CCD; 

(2)  surface  channel  charged  coupled  devices,  SC-CCD;  (3)  the 
charge  injection  device,  CID;  and  (4)  the  self-scanned  diode 
array,  SSDA.  Some  of  the  characteristics  of  the  first  three 
of  these  are  summarized  in  Tables  6  and  7. 

The  devices  that  are  commercially  available  are  rela¬ 
tively  large  area  arrays,  100  x  100,  being  typically  the  smallest 
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TABLE  6.  COMPARISON  OF  SURFACE  CHANNEL  AND 


k 


BURIED  CHANNEL* 


SURFACE  CHANNEL  BURIED  CHANNEL 


Limiting  temporal  noise  Surface  trapping  Bulk  trapping 

>  100  electrons  >  10  electrons 

Limiting  pattern  noise  Nonuniform  Nonuniform  dark 

background  current,  limp, 

>  100  electrons  sensitive 

Detectivity  limit  >  100  electrons  >  10  electrons 

Speed  limit  <  10  MHz  25  MHz  t 


TABLE  7.  TRADEOFFS  BETWEEN  FRAME  TRANSFER  CCD, 

INTERLINE  TRANSFER  CCD,  AND  CID  AREA  ARRAYS* 

FRAME  TRANSFER  INTERLINE 

CCD  TRANSFER  CCD  CID 

Versatility  Front  or  back  Front  only  Random  read- 

illumination  illumination  out,  front 

i 1 lumination 

Sensitivity  High  silicon  Lower  quantum  High  quantum 

quantum  efficiency  efficiency 

ef  f i c iency 

Problems  High  speed  Complex  cell  Pattern  noise 

vertical 
transfer 


*D.F.  Barbe,  "Charge  Coupled  Device  and  Charge  Injection 
Device  Imaging,"  IEEE,  Vol  EI)23,  177  (  1976). 
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in  CCD  arrays,  and  32  x  32  in  the  SSD  arrays. 

Devices  tuat  have  as  few  as  256  elements  that  can  be 
read  out  at  the  required  rates  are  limited  to  linear  arrays. 

It  is  possible  to  obtain  custom  arrays  at  a  cost. 

We  have  not  pursued  this  except  for  an  anode  array  discussed 
below  in  Section  2  because  we  found  no  reason  to  do  so. 

The  structures  of  four  area  arrays  are  shown  in 
Figures  36,  37  and  38.  The  first  of  these  is  a  buried  channel 
CCD,  the  Fairchild  190  and  244  area  imager.  Since  the  photosite 
region  is  separated  from  the  transfer  lines,  the  device  is 
relatively  inefficient.  In  this  case,  the  image  area  is  14  x  16 
providing  a  sensitive  area  utilization  of  only  46^.  There  is  a 
second  problem;  only  every  other  photosite  is  read  out  at  any 
one  time.  The  array  is  interlaced,  all  of  the  even  numbered 
cells  "A"  are  read  out  in  one  frame,  and  then,  once  the  transport 
channels  have  been  cleared,  the  other  frame  would  be  read  out. 

The  result  of  such  a  readout  would  be  two  sots  of  data  which 
would  have  to  be  processed  differently.  One  set  of  data  would 
integrate  from  0°to  90°,  the  other  from  -45°  to  +45°,  then  90°  to 
180°,  and  45°  to  135°. 

The  geometry,  providing  two  cells  in  one  direction 
and  only  one  cell  in  the  other,  is  difficult  to  use  because  the 
subapertures  are  square.  Thus,  two  cells  and  a  blank  space 
would  cover  one  subaperture.  On  chip  addition  could  be  used 
to  consolidate  the  subapertures.  Instead  of  shifting  out  all 
of  frame  A  before  reading  'rame  B,  one  could  shift  one  pulse 
vertically,  then  transfer  frame  B,  thus  adding  A  and  B  into 
the  same  bucket.  This  technique  was  successfully  used  by  the 
author  in  a  CCD  wavefront  sensor  for  another  project. 

The  readout  of  a  32  x  16  array  could  1  hen  proceed  at 
a  10  MHz  rate  or  slightly  faster  to  make  up  for  the  transfer 
t imes . 
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Buried  channel  TCP  structure. 
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The  .subaperture  areas  would  be  30  x  36  am.  The  trans¬ 
fer  of  frames  A  and  B  would  occur  just  before  and  after  the 
desired  transfer  time  to  average  the  errors  out. 

The  buried  channel  device  with  its  capability  of 
adding  the  contents  of  a  photos i t e  to  a  transport  line  bucket  and 
then  to  read  out  at  high  speed  allows  the  array  to  be  used  as 
a  smaller  array.  If  only  one  corner  of  the  array  is  used, 
and  if  a  part  of  transfer  pulse's  are  given  after  256  vertical 
and  horizontal  transfers,  then  the  area  observed  will  be  the 
sum  of  the  sets  of  32  x  16  elements  that  can  be  extracted.  The? 
pulsing  sequence  is  shown  in  Figure  38  where  the  method  is  shown 
for  a  4  x  4  subaperture  set  out  of  a  20  x  8  array.  The  array  is 
masked  so  that  only  the  lower  left  hand  corner  can  receive  an 
exposure.  All  of  the  other  photosites  will  generate  a  dark 
current  which  would  add  to  the  readout  and  which  will  need  to 
be  suppressed  by  cooling.  The  readout  sequence  is  initiated  by 
shifting  all  of  the  charge  in  the  even  cells  into  the  vertical 
shift  buckets.  Then,  after  one  vertical  shift,  all  of  the  odd 
cell  charges  are  added  to  the  even  cell  charges.  The  next 
vertical  shift  brings  (he  line  of  the  first  sum  into  the  hori¬ 
zontal  buckets.  This  line  is  read  out  using  four  pulses.  At 
this  point,  the  fifth  horizontal  bucket  is  the  first  position 
and,  on  the  next  vertical  transfer,  will  receive  the  charge 
contents  of  vertical  cells  3  and  4.  This  signal  will  be  added 
to  the  dark  current  from  colls  !  and  2  in  column  f>.  If  this  is 
small,  the  addition  will  be  trivial  and  the  result  will  be  as  if 
a  4  x  4  array  had  been  used.  There  is  very  little  time  lost  in 
the  process.  Three  pulses  an*  required  for  cell  addition,  4  for 
vertical  and  16  for  horizontal  transfers,  a  total  of  23  lime 
increments  for  16  elements.  In  a  16  x  16  array,  we  st  ill  have 
3  for  addition,  16  for  vertical  and  256  horizontal  shifts, 
requiring  275  pulses  in  all.  The  result  is  an  dlociive  increase 
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in  the  cell  size  by  2X,  a  reduction  in  the  size  of  the  array  at 
a  trivial  cost  in  speed.  The  process  can  be  extended  to  enlarge 
the  added  pixel  size  further  to  yield  4x3  element  subapertures, 
for  example,  but  at  a  cost  reduction  of  3  in  video  bandwidth 
(but  not  12,  one  dimension  being  compensated  but  not  the 
horizontal),  which  is  fixed  by  the  readout  rate. 

The  smallest  Fairchild  area  array  presently  available 
is  the  CCD  211  with  244  x  190  elements  with  18  pm  vertical  and 
30  ,jm  horizontal  spacings.  With  two-fold  additions,  the  cell 
size  would  be  36  x  30  ym.  The  collection  efficiency  is  47%. 

The  dark  current  in  each  output  bucket,  because  of  the 
repetition  addition  averages  over  about  85  elements,  is  the 
same  as  it  would  be  in  normal  video  use,  and  would  be  0.3%  of 
saturation . 

The  Fairchild  244  x  190  device  can  be  driven  at  15  MHz. 
The  second  configuration,  the  surface  channel  CCD,  is  that  of 
the  RCA  imager.  In  this  device,  there  seems  to  be  no  structure 
at  all  because  the  photosites  are  the  transfer  channels.  The 
device  is  100%  efficient  because  there  are  no  covering  elements. 

Alternate  photosites  "A"  accumulate  the  charge.  The 
intervening  photosites  "B"  are  biased  slightly  negative  so  that 
charges  generated  in  these  regions  will  end  up  in  the  adjacent 
site  A  or  A'.  Normally,  the  readout  is  interlaced,  A  then  B, 
but  there  is  no  need  to  do  this  here,  as  there  is  in  the  BCCCD. 

The  one  difficulty  in  the  method  of  transport  is  that 
if  the  device  is  not  blanked  during  frame  transfer  to  the  storage 
area,  charge  will  be  accumulated  during  transfer,  creating  a 
smear.  In  a  16  x  16  element  array,  the  smear  time  will  be  1/16 
of  the  frame  rate.  Such  smear  may  limit  phase  accuracy  in  one 
dimension  because  the  charges  accumulate  only  in  the  biased 
photos i to . 
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The  biggest  difficulty  with  the  device  is  that  it  is 
slow.  The  3  MHz  video  rate  would  restrict  the  device  to  75 
subapertures  or  less.  Finally,  tricks  such  as  those  used  on  the 
buried  channel  device  cannot  be  used. 

The  third  device,  the  CID  manufactured  by  General 
Electric,  and  the  fourth,  the  self-scanned  diode  array  made  by 
Reticon,  are  different  from  the  CCD's  in  that  the  photosensitive 
elements  must  be  addressed  sequentially.  The  entire  field  can¬ 
not  be  addressed  at  one  time,  as  it  can  with  the  CCD's.  The 
problem  with  sequential  samples  is  that  it  places  restrictions 
on  the  algorithms  that  can  be  used.  This  problem  will  be  dealt 
with  in  the  next  section. 

The  structure  of  the  CID  is  shown  in  Figure  39.  The 
design  of  a  high  speed  device  is  shown  here.  The  operation  of 
readout  consists  of  transferring  a  charge  from  one  capacitor, 
the  photosite,  to  another,  and  measuring  the  displacement 
current.  In  the  high  speed  device,  the  signal  charge  for  a  whole 
line  is  transferred  in  parallel  to  the  column  electrodes  which 
are  then  scanned  sequentially. 

There  is  a  problem  of  lag  due  to  residual  charge  if 
the  injection  is  insufficient.  This  takes  time  and  is  the 
limiting  factor  in  speed  of  operation  5  MHz  seems  to  be  the 
upper  limit. 

The  photosite  area  of  one  32  x  32  array  consists  of 
43  x  43  um  picture  elements.  This  is  a  slow  array,  however, 
because  of  the  capacitance.  There  is  not  enough  data  on  avail¬ 
able  arrays  to  suitably  characterize  the  device  for  the  AWSC 
application. 

The  Reticon  self-scanning  diode  array  is  schematically 
shown  in  Figure  40.  Like  the  CID,  each  diode  is  addressed  by  an 
x  and  y  gate.  The  charge  is  sensed  directly,  however,  and  the 
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Diagram  of  baric  X-Y  aecauinf  schema  for  a  CID  imager,  (a) 
Schematic  diagram  of  a  4  x  4  array,  (b)  Se ruing  rite  noe  aectioa 
showing  rilicon  turf ace  pottntialt  aad  location  of  Rond  charge  for 
various  operating  conditions. 


Figure  39.  CID  block  diagram. 


*  H.  K.  Burke  and  G.  J.  Michon,  "Chargo-In.iect  ion  Imaging: 
Operating  Techniques  and  Performance  Characteristics", 
Electron  Devices,  Vol .  ED-23,  No.  2  (Feb  76),  p.  190. 


1  15 


Figure  40.  SSDA  block  diagram,  RL-256 


♦Reticon  Corporation  data  sheet  27100 
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diode  is  completely  discharged.  All  of  the  photoelectrons  are 
collected  on  the  diodes  which,  in  the  32  x  32  array,  are  100  pm 
apart.  In  the  area  array,  the  sample  rate  is  5  MHz.  The  quan¬ 
tum  efficiency  is  high  because  there  are  not  surface  electrodes 
which  can  create  reflection  losses,  as  in  the  CID.  The  basic 
problem  with  the  area  array  Is  the  high  video  line  capacitance 
of  100  pf.  This  produces  a  larger  ICTC  noise  and  a  NES  of 
about  500  electrons. 

The  problem  with  all  of  the  area  devices  is  the  un¬ 
availability  of  small  arrays.  Only  the  buried  channel  CCD  can 
be  used  like  a  small  array.  Other  than  this  case,  a  custom 
design  at  an  appropriately  high  cost  would  be  required,  which 
may  not  be  justifiable  in  terms  of  the  few  devices  that  would 
actually  be  used  in  compensated  imaging  systems. 

There  are  two  linear  devices  that  have  256  elements 
that  can  be  read  out  at  rates  of  10  MHz.  These  are  the  Fair- 
child  CCD  110/110F  linear  image  sensor,  shown  in  Figure  41  and 
the  Ret  icon  RL  256S  SDDA's,  shown  in  Figure  40.  One  version  of 
the  Reticon  C  series  arrays  is  expressly  designed  for  electron 
bombardment . 

The  CCD  photoelement  is  small,  13  x  17  pm.  The 
Reticon  element  which  has  been  specially  designed  for  spectros¬ 
copy  is  25  pm  wide  center  to  center,  and  430  Pm  long. 

Both  devices  are  equally  efficient  optically,  with  all 
of  the  photoelectrons  collected  in  the  photosite  on  a  diode. 

Both  devices  use  two  video  lines  carrying  odd  and  even 
elements  separately.  The  big  difference  in  the  devices  is  that 
all  of  the  elements  are  shifted  into  the  transfer  line  simulta¬ 
neously  in  the  CCD  array. 

In  the  Reticon  device,  the  diodes  are  read  sequentially 
so  that  the  integration  period  of  each  starts  and  ends  at  differ¬ 
ent  point  in  a  nutation  cycle. 
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Figure'  41.  Fairchild  linear  CCD  block  diagram. 


6.3  Signal  Processing 

The  basic  array  options  are  the  two  linear  arrays,  the 
Fairchild  CCD  111,  which  can  be  photon  coupled,  and  the 
Ret  icon  array,  which  can  be  used  in  the  EB  mode,  and  one  area 
array.  The  Fairchild  CCD  211,  or  a  CCD  202,  a  100  x  100  element 
array  which  is  not  made  any  more,  can  be  used  with  the  small 
area  sampling  technique  described  above,  and  would  be  photon 
coupled . 

The  signal  processing  that  would  be  used  for  the  CCD's 
would  be  about  the  same,  whether  linear  or  square,  except  for 
the  transfer  line  clocking.  In  the  photon  coupling  mode,  there 
will  be  a  delay  due  to  phosphor  decay,  which  will  affect  frequency 
response.  But,  this  will  not  affect  the  use  of  the  array.  All 
readouts  will  be  done  together  in  synchronization  with  the  nuta¬ 
tion  cycle. 

The  signal  processing  for  the  Reticon  array  is  more 
complicated  because  of  the  fact  that  the  diodes  are  read  out 
sequentially.  This  is  the  major  concern.  Otherwise,  the 
signal  processing  is  quite  similar  to  the  CCD  signal  processing. 

6.3.1  CCD  signal  processor.  The  CCD  detector 
selection  is  governed  principally  by  the  data  readout  rate. 

For  10  kHz  nutation,  there  are  four  readouts  per  photocell  per 
nutation  cycle.  For  256  subapertures,  a  10.24  MHz  readout  is 
required.  This  constraint  immediately  eliminates  surface  channel 
CCD's  because  the  transfer  inefficiency  increases  dramatically 
above  5  MHz. 

The  buried  channel  CCD  (BC-CCD)  architecture  can 
easily  support  the  10  MHz  readout  rate.  The  potential  difficul¬ 
ties  in  BC-CCD's  include  dynamic  range,  dark  current,  cell  uni¬ 
formity,  readout  noise,  and  preamplifier  noise.  These  will  be 
dealt  with  separately  below,  using  the  Fairchild  CCD  111  device 
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characteristics.  It  will  be  shown  that  the  preamplifier  is  the 
principal  area  of  concern. 

Dark  si gnal .  There  are  random  electron  hole 
generations  in  the  doped  silicon  used  to  construct  CCD's. 

These  minute  charges,  called  dark  signal,  are  collected  in  the 
cells  exactly  as  photoelectrons  would  be  collected.  The  statis¬ 
tics  governing  this  charge  generation  indicate  that  total  dark 
signal  is  proportional  to  integration  time  and  is  strongly  depen¬ 
dent  on  temperature. 

For  a  10  kHz  nutation  rate,  the  integration  time 

is  25  ys.  For  a  CCD  111  at  25°C,  this  integration  time  generates 

-4 

a  maximum  dark  signal  of  3.2  x  10  of  V  ,  the  saturation  voltage. 

s 

This  dark  signal  will  halve  for  each  approximately  -18°C  in 
temperature . 

Cell  uniformity.  CCD’s  are  fabricated  on  rela- 

2 

tively  large  (approximately  1.2  mm  )  silicon  dies  using  multiple 
mask  exposures.  Large  area  diffusion  inhomogeneities  and  mask 
misalignments  lead  to  nonuniform  responsi vit ies  and  nonuniform 
dark  signals  for  the  various  CCD  cells  and  photosites.  Dark 
signal  nonuniformity  in  the  CCD  111  is  a  maximum  of  2  5%  of  the 
saturation  voltage.  This  nonuniformity  is  not  due  to  time 
domain  statistical  fluctuations  but  rather  to  device  inhomogene¬ 
ities.  Thus,  dark  signal  nonuniformity  is  removable  in  real  time 
as  a  separate  dark  signal  parameter  is  stored  for  each  cell. 

The  responsivi t ies  of  the  photosites  vary  by  a 
maximum  of  ±6%.  Again,  this  is  not  a  time-varying  parameter  and 
hence  is  removable  using  a  cell-by-cell  scaling  factor  in  the 
signal  processor.  Using  digital  scaling  with  8-bit  accuracy, 
responsivity  uniformity  can  be  reduced  to,  at  most,  ±  0.2%. 

It  should  be  noted  that  this  scaling  will  most 

3 

likely  not  be  necessary.  The  I  sensor  measures  wavefront  tilts 
relative  to  an  optical  reference.  The  optical  path  and  electronic 
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processing  are  identical  for  the  source  and  reference  beams  and 
hence  any  responsivity  scaling  factors  operate  equally  on  two 
measurements.  If  both  measurements  are  scaled  or  biased  equally, 
then  the  difference  signal  will  still  indicate  a  proper  null  for 
servo  control  systems.  Thus,  uncompensated  responsivity  nonuni¬ 
formity  should  have  no  effect  worse  than  approximately  10% 
uncertainty  in  the  absolute  magnitude  of  the  subaperture  tilts. 

Dynamic  range.  CCD  dynamic  range  is  specified 
as  the  ratio  of  the  saturation  voltage  to  the  peak-to-peak  noise 
equivalent  exposure.  This  specification  does  not  take  into 
account  the  dark  signal  or  the  dark  signal  nonuniformity.  For 

the  CCD  111,  the  typical  dynamic  range  is  500. 

3 

The  AWSC  I  sensor  will  be  designed  for  photon 
noise  limited  performance.  The  mi crochannel  amplifier  should 
amplify  the  incoming  signal  to  generate  at  least  twice  the 
peak-to-peak  noise  equivalent  exposure  for  each  event,  or  about 
5,000  electrons  per  photon.  This  is  8  -  12  times  the  rms  CCD 
noise  signal  and  hence  reliable  detection  of  a  single  photon 
is  possible. 

The  dark  current  and  dark  current  nonuniformity 
also  affect  dynamic  range.  The  dark  current  itself  is  negligible 
however,  the  worst  case  nonuniformity  of  2.5%  saturation  voltage 
results  in  a  dynamic  range  of  40.  This  worst  case  value  can  be 
made  a  negligible  1/2  the  noise  equivalent  exposure  by  cooling 
the  CCD  approximately  60°C  to  -35°C.  Note  that  typical  dark 
signal  nonuniformt iy  is  0.5%  of  saturation.  Cooling  the  CCD  by 
about  18°C  to  +7°C  will,  in  general,  be  adequate  to  remove  all 
dark  current  effects  from  the  calculations.  This  is  well  within 
the  capabilities  of  thermoelectric  cooling  mechanisms. 

By  removing  dark  current  effects,  a  dynamic  range 
of  0  -  250  photons  is  possible  within  the  CCD.  The  system 
dynamic  range  is  at  least  this  value  times  the  AGC  capabilities 


12] 


of  the  microchannel  amplifier.  With  readouts  of  40  kHz,  this 
yields  a  peak  flux  of  5  n.  The  AGC  range  needs  to  be  only  10:1. 

Readout  noise.  The  transfer  of  charge  from  one 
cell  to  another  is  not  100%  efficient  in  a  CCD.  An  electron 
left  behind  by  one  cell  will  most  likely  be  picked  up  by  the 
next  cell>-and  hence  transfer  noise  is  correlated  from  one  cell 
to  the  next  but  is  uncor related  from  one  readout  cycle  to  the 
next . 

Typical  charge  transfer  inefficiencies  for  BCCD's 
-5  -4 

are  10  to  10  .  The  total  transfer  is  then  S  •  N  •  i  (for 

N*  e  <<  1),  where  S  is  the  total  signal,  N  is  the  number  of 
transfers,  and  c  is  the  charge  transfer  inefficiency.  The  CCD 
111  uses  two  output  CCD's  for  the  256  photosites  and  hence 
N  =  130  transfers  (the  two  extra  transfers  are  to  move  the  charge 
from  the  photosites  to  the  CCD's  and  from  the  CCD's  to  the  out¬ 
put).  Thus,  0.0013  1.  N  0.013  and  for  a  worst  case  of  N't  = 

0.013  and  S  =  saturation  voltage  (  250  source  photons  by  system 
design),  the  transfer  noise  of  0.013  x  saturation  is  much  less 

_  i  /9 

than  the  photon  noise  of  (250)  ' “  =  0.063  x  saturation.  For 

the  photon  noise  limited  system,  transfer  noise  in  a  BCCD  is 
negligible. 

Another  source  of  readout  noise  is  reset  clock 
feedthrough  in  the  output  charge  detector-preamplifier.  The 
charge  detector  is  a  diode  which  is  precharged  to  a  fixed  voltage 
by  a  reset  clock  signal.  The  diode  potential  then  changes 
linearly  with  the  quantity  of  charge  delivered  from  a  ('CD  bucket. 
This  potential  is  applied  to  the  gate  of  an  MOS  output  transistor. 

The  desired  signal  is  the  difference  between  the 
MOS  transistor  output  during  precharge  and  the  MOS  transistor 
output  after  the  CCD  charge  has  been  accumulated.  Any  change  m 
reset  clock  duty  cycle  or  reset  voltage  will  direct  ly  al led  the 
output  during  precharge  and  hence  this  precharge  value  must  be 


sensed  dynamically. 

The  accepted  method  for  this  dynamic  sensing  is 
correlated  double  sampling  (CDS).  At  10  MHz,  CDS  has  the 
disadvantage  or  requiring  accurate  capacitor  charge-discharge 
in  25  ns.  The  Fairchild  CCD  111  removes  this  necessity  by 
supplying  an  on-chip  dummy  output  which  involves  a  precharge 
diode  and  MOS  output  transistor  identical  to  the  on-chip  CCD 
preamplifier.  Thus,  the  reset  signal  and  the  output  signal 
are  available  simultaneously  and  simple,  nonstorage  subtraction 
is  adequate  to  recover  the  signal.  This  output  structure  is  a 
strong  positive  feature  of  the  CCD  111. 

Preampl i f ier .  The  CCD  111  is  2,3  cm  long  by 
0.7-J  cm  wide.  This  results  in  the  two  5  MHz  transport  clocks  and 
the  10  MHz  reset  clock  being  physically  very  near  (approximately 
0.2  cm)  the  sensitive  output  and  dummy  output.  Capacitive  coup¬ 
ling  of  the  high  speed  clocks  to  the  output  is  a  sever  problem. 
However,  this  problem  can  be  reduced  with  good  engineering  and 
careful  layout;  no  technological  developments  are  required.  The 
newer  110F  has  this  clock  at  a  different  location. 

All  hardware,  even  breadboard,  will  require 
double-sided  printed  circuit  boards  incorporating  a  ground 
plane.  Isolation  etching  around  output  pins  will  be  required. 

The  IK  ohm  output  impudences  will  immediately  be  lowered  to 
approximately  10  ohm  using  matched  pair  emitter  followers. 

The  output  and  dummy  signals  can  then  be  sent  to  a  commerical 
quality  differential  amplifier.  The  approximately  400  pV  output 
noise;  voltage'  of  the  CCD  11  L  will  dominate  any  diffential  ampli¬ 
fier  noise'.  After  amplification,  the;  d  i  f  f  ere'nt  i  a  1  signal  can 
safe-ly  be’  sample;d  and  hc>ld  and  conve rteel  to  an  8-bit  digital 
s  i  gna  1  . 
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C.3.2  Reticon  signal  processing.  This  siit,  •  ■<  :  mu < 
be  considered  in  two  parts,  the  electronic  processing  r.  quit.  - 
ments  in  which  much  is  similar  to  the  CCD  processing  i-cqu  i  r<*m«  n  t  s , 
and  the  algorithm  dependent  requirements. 

The  readout  requirements  of  the  Reticon  256  cle¬ 
ment  self-scanned  diode  array  are  substantially  the  same  as  i he 
requirements  for  the  256  element  CCD.  The  dark  signal  died, 
cell  to  cell  nonuniformity,  and  dynamic  range  considerations  arc 
identical.  The  bandwidths  and  sensitivities  of  the  detector 
preamplifiers  are  also  comparable. 

The  major  difference  between  the  CCD  and  the 
diode  array  is  the  readout  noise  median  ism.  The  CCD  sutlers 
from  transfer  inefficiency  and  the  diode  array  does  not  How¬ 
ever,  the  high  output  capacitance  (approximately  25  pf)  of  the 
diode  array  will  be  a  problem  at  the  10  Mil/,  readout  frequency. 

The  25  pf  capacitance  must  be  totally  discharged  within  100  ns 
and  this  may  require  a  cooled  preamplifier.  The  preamp! i tier 
must  possess  an  effective  input  resistance  at  10  MHz  at  less 
than  100  ohms. 

The  effect  of  sequential  sampling  on  the 
algorithms.  Because  this  problem  was  analyzed  before  t  lie 
optimum  algorithm  had  been  developed,  only  the  IS  I'  and  MSI1' 
algorithms  were  studied.  The  analysis  assumed  a  square'  Mat 
distribution  object, 'a  nutation  radius  less  than  the  halt  width, 
and  a  large  quad  cell.  This  analysis,  which  showed  that  the 
sampling  techniques  created  an  effective  rotation  ol  the  coor¬ 
dinate  system,  also  showed  that  it  did  not  affect  the  errors. 

3* 

The  MSI  algorithm,  however,  was  shown  to  be  defective,  and  could 
not  be  used  if  the  sampling  techniques  were  employed.  The  cal¬ 
culation  for  the  optimum  algorithm  will  follow. 

In  the  Reticon  arrays,  all  diodes  un¬ 
connected  to  a  single  video  line  through  gates.  These  gates 


I 
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are  turned  on  sequentially,  reading  one  diode  charge  at  a  time. 
All  of  the  other  diodes  are  integrating  while  one  is  being 
output.  Thus,  the  output  of  the  nth  cell  with  a  light  input 
I(t,x)  where  x  is  the  location  along  the  array,  is 

n  A  t 

In(tQ)  =  /  I ( t , x ) d  t  (31) 

(r-N)At 


where  N  is  the  number  of  elements  in  the  array  and  At  =  l/f  is 

the  readout  interval.  The  output  is  obtained  at  t  =  nAt. 
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In  the  use  of  the  I  sensor,  four  readouts 
must  be  obtained  for  each  nutation  cycle. 

The  optical  input  to  the  detectors  is  given 


by 


4  p  Axsin6  +  4pAycos0  +  4p  sin2d  +  DC  terns 


(32) 


where 

0  =■  2tt t /T  (with  T  =  nutation  period) 


(33) 


I  will  assume  four  arrays  operating  simultaneously,  with  the 
outputs  of  cells  I  —  II  +  III  —  IV  summed  as  usual. 

Three  correlating  outputs  are  obtained  in 
time  sequence.  These  are  weighted  by 

+-~+,  +-+- 

For  the  n^  detector,  the  periods  are 


1  2.3 


(irn/2N)  ±  (  tt  /4  ) 

(  tt  n  /  2N  )  +  (ti  /2)  +  (r/4) 

(  ’  /2N)  +  (r)  ±  (r/4) 

(  ”  n  /2N  )  +  (  dr/ 2  )  +  (r/4) 


The  integrals  are 


4  2N+".i+4  0 

4  i  V>  p  j  (ixsinP  +  AveosO  +  o“s  i  n2*1  ) d’ 
,)  =  1  ’  ^  2 
2N  ,i“4 


For  different  G’s  in  sequence,  the  weights 
in  the  sum  are  given  by 


WEIGHT 


w  =  s 

.)  .) 


v  C 
.)  .1 


v;  =  2s 


With  the  S  weight  we  obtain 


4 ,  2 ,  H 


{,  £sin(r*,/2N)  +  cos (  ” - .  / 2\1]  +  ■>[[<•< 


' 2 N  )  -  m(- 


and  with  t  he  (’  we  1  gh t 


4  2vT? 


H 


s  1  i)  (  •■  ’  1  '  2  N  -  e  o  s  (  "  ’  2  \  )  I  ♦  .'.v  aisl-  2\ )  *  sin(-  2\ ) 


These  can  be  rewritten  as 

J*S  =  16p  ^ixsi  n(  rri /2N'  +  7/4)  +  Aycos(  -rr /2S  +  r/d")] 


lfin  U'  ' 


fC  -  16p I -Axcos( -n / 2N  +  -/4)  +  ivsin( t /2N  +  - /4 ) 


16- lx' 


As  before,  the  /S  yields  the  Ay  term  and  fC  the  Ax  term. 

The  denominators  must  now  be  obtained.  First, 
to  obtain  the  MS  denominator,  simply  apply  the  2S  weight,  which 
y i e 1 ds 

/ 2S  -  16  p  2sin(-»  (40) 


This  is  not  usable  at  any  angle,  but  at  0  or 
n  =  N/2.  At.  the  end  point  it  yields  a  zero.  The  rate  of 
change  of  the  null  is  twice  as  fast  as  the  change  in  the  other- 


terms  . 


The  denominators  in  the  standard  I1  formula¬ 


tion  are  better  behaved.  Adding  the  terms  I  +  11  -  111  -  IV  for 
y  denominator  and  using  the  S  weighting,  we  have 


— i  •'  +  — 

2N  .)  2 


AS  •  2 we 
.1 


s  i  n"d- 


2. •  2  ’  s  1  n  (  -  ’  1  2N  )  +  <-osi--/2N) 


1 


(H) 


8wpsi n(  irri/2N  +  n/4) 


There  is  a  cross  talk  induced  by  the  phase 
of  the  gates.  If  we  transform  the  deviations  into  a  coordinate 
system  for  each  pixel  determined  by  the  turn  off  angle  for  the 
gate 

X  =  tt n  / 2 N  +  tt/4  (44) 

and  the  turn  on  angle 

X  =  it/2  (45) 

then 

/ S  =  1  f>p  (  A x s i  n x  +  Avcosy)  (4b) 
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and 


fC  =  16p(-Axcosx  +  Avsinx) 


(47) 


The  displacements  Ax  and  are  obtained  by 
the  inverse  transformation.  Since  the  ratios 


Rx  =  (2/w)(Ax  +  Aycotx) 


and 


(48) 


Ry  =  (2/w)(Ay  -  AxcotX) 


(49) 


we  obtain 


x  =  sinx  *  (w/2) (R^sinX  -  RycosX) 


(50) 


Does  this  change  the  errors  in  the  measurement? 

Since 


<J>  =  x  / 1 , 

X  '  ’ 


(51  ) 


the  phase  variance  is 


2  ,  .2  .  2  ,  2.2  ,  2  2  ,  .  ro, 
a.  =  ( w  /  2  A  )  sin  X(°R  sin  X  +  °r  cos  *)  (o2) 

* x  x  y 


Since 


2 

a  num  x 
D  2 


and  o , 


2 

a  num  v  (53) 


D 
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then 


a.  2  =  (w/2A)2sin2x(o  2/D2)(l  +  R2 )  (54) 

If  we  normalise  the  oblect  brightness  and  integration 

time,  we  can  express  the  ob/joct  size  in  terms  of  photon  counts. 

o 

Let  this  value  be  N  =  vw  =  w  ,  if  square.  Since  (55) 

D  =  8wpsinx  =  8N(p/w)sinx  (56) 

The  variance  on  the  numerator  for  this  flux  is 


Since 


var  /  S(  I-II  +  III-IV)  =  fS2(a12  +  on2  +  om2  +  o^2) 

(57) 


o2  =  Nj  =  wv/4  +  ( w2/2 ) ( p/w)sin0  +  ( v/2 ) ( p / v ) cos0 

(58) 

2 

Because  S  =  1 ,  it  does  not  matter  where  the  integration  begins; 
consequent ly 

var (num)  =  2rM  =  (59) 

and 

N 

o  2  =  (w/2\)2sin2x  - 5 - 5- - o - 2~  (60) 

TA  (8/27r)z(NT‘i)(p.w)^sinzx 


at  the  minimum 

p/w  =  1/2 


o 


(w/2A  )2(  1/N)(2tt2/8) 


(61  ) 


2 


(62) 


and  we  find  that  the  accuracy  is  independent  of  where  the  sub¬ 
aperture  is  sampled. 

The  optimum  algorithm  with  CCD  and  reticon 
readout .  The  optimum  algorithm  weighting  functions  are  composed 
of  DC  and  AC  terms.  This  introduces  a  disparity  between  DC  and 
AC  results  because  the  axis  of  the  AC  measurement  may  be  rotated, 
but  the  DC  axis  would  not.  In  order  to  evaluate  this  effect,  a 
power  dependent  calculation  must  be  made, 

First,  the  CCD  readout  of  a  nutated  square  object 
using  an  optimum  algorithm  will  be  derived  and  then  the  algorithm 
will  be  modified  for  the  Reticon  and  the  results  will  be  compared. 

The  actual  photon  flux  on  the  detectors,  fol¬ 
lowing  the  notation  of  Figure  42  to  define  the  geometry  of  a 
rectangular  object  with  width  w  and  depth  v  with  a  brightness  B, 
so  that  the  total  number  of  photons  over  a  nutation  period  is 

2  IT 

Nt  =  nBvw/hv  =  /  ZN'.de  (63) 

is  given  by 

=  (r,B/hv2it)  {  ( w/2  +  Ax) (v/2  +  Ay)  + 

+  (w/2  +  *.x)(psin4)  +  ( v/2  +  Ayjpe.osO  + 

9 

+  ( p“si nOcosP )  f 

N2  =  (r;B/hv2r)  ‘  (v;/2  -  Ax)  (v/2  +  Ay)  + 

+  (w/2  -  Ax  )  ( :>s  i  nn  )  +  (v/2  +  Ay  )  (  pens"  ) 

2 

-  (p  sinOcosO)  ) 


1  3  I 
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Ng  =  (nB/hv2Tr)  f  (w/2  -  Ax)(v/2  -  Ay)  -  (w/2  -  Ax)(psin9) 

o 

-  (v/2  -  Ay)(pcos0)  +  (p  sin0cos9)  } 

=  (nB/hv2Ti)  {  (w/2  +  Ax)(v/w  -  Ay)  -  (w/2  +  Ax)(psin6) 

o 

+  (v/2  -  Ay)(pcos0)  -  (p  cos0sin8)  } 

(64) 


One  set  of  optimum  weighting  functions  for  the 
numerators  and  denominators  is  shown  in  Table  8. 


TABLE  8 


for  Ax 

for 

Ay 

s. 

=  1  + 

sin  0 

s. 

_ 

1 

+ 

cos 

lx 

iy 

S2x 

=  -1  - 

sin  0 

So 

2y 

1 

- 

cos 

S3x 

=  -1  + 

sin  9 

S3y 

-1 

+ 

cos 

S4x 

=  1  - 

sin  0 

S4y 

-1 

- 

cos 

and 

Dlx  = 

1  +  sin 

0 

lv 

= 

1 

+ 

cos 

D2x  = 

1  +  sin 

e 

°2y 

= 

1 

- 

cos 

D3x  = 

1  -  sin 

0 

D3y 

= 

1 

- 

cos 

D4x  = 

1  -  sin 

e 

D4y 

= 

1 

+ 

cos 

TABLE  9 
CCD  WEIGHTS 


Mean  angle 
limits  of 
integration 
weights 
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In  CCD  processing  the  samples  are  read  out 
whenever  the  nutation  is  at  a  45°  angle  and  a  different 
weighting  is  applied  to  each  depending  on  the  mean  angle  of 
the  integrated  sample.  Thus,  the  weights  are  as  shown  in 
Table  9. 

Applying  these  weights  to  the  detector 
outputs  and  integrating  over  a  nutation  cycle,  we  obtain  the 
x  nume  rators : 


2  TT 

/  l  N.S. 

0  ill  1  1 

=  (nB/hv)^2Axv 

+  (8/2it)/7Ax7| 

(65) 

and  the  denominator 

2"  4 

/  /  N.D. 

J  . L .  l  i 

0  1  =  1 

=  (pB/hv)(wv  + 

4vi1Z/2t.  +  wf  ) 

(  66  ) 

so  that  the  ratio  is 

r  =  2Ax/w  (67) 

Thus,  the  displacement  is 

Ax  =  rw/2X  ( fib ) 

and  the  phase  is  wave's  is 


A  $ 


rw/2 


(  C)!* ) 


The  variance  on  the  phase  is 


2  2  2 
a.  =  (w/2X)  a  * 

h 


where  the  variance  of  the  ratio  is  given  by 


a 


2 

r 


2  2 
a  /(Den) 
num  ' 


— ' 

2 

_ 


(1  + 


r2) 


(70) 


(71) 


The  variance  of  the  numerator  is 


2  IT 


num 


2  -  [  K,2s, 


.  N.  ~  i 
1  =  1  o  1 


where 


(72) 


o 


N. 

l 


2 


N. 

l 


N  =  nBvw/hv 
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Thus,  the  phase  variance  is 


When  this  calculation  is  applied  to  the  Reticon 
readout,  two  values  are  obtained  for  the  AC  and  DC  calculations. 
The  integrations  are  over  the  intervals 

TTH/2N)  +  (tt/4)  *  (I  -  1  )(ti/2)  (74) 

where  I  is  the  segment ,  n  the  number  of  the  diode  and  N  the  num¬ 
ber  of  diodes. 

If  this  is  done  using  the  weighting  functions 
described  above,  the  results  are  discordant. 

The  DC  term  is  not  rotated,  but  the  AC  term  is. 

In  the  case  of  a  square  object  used  here,  and  using  the  results 
of  the  preceding  discussion  substituted  into  the  optimum  algorithm 


i  :ss 


AC  term,  we  would  find 


(  2Axv  )  +  p  (  8/2/2tt  )L£3cs  i  n(  rr./2N  +  n/J_)  +  Avcosf  -r  /2N  +  W  lT] 


wv  +  ^4/2’/2tt  )  pws  i  n( -n  n /2N'  +  r/dfj 


(  75 ) 


The  result  is  now  complicated  not  only  by  p  but 
by  shape  factors  w  and  v,  which  in  general  are  unknown.  The 
rotation  cannot  be  removed  easily.  If  the  weighting  function 
were  to  contain  the  phase  of  the  readout  location,  the  derotation 
may  be  obtained  before  the  ratio  is  obtained.  A  similar  set  of 
functions  to  those  used  with  the  CCD  but  the  mean  angle  of  inte¬ 
gration  would  be  used  for  weighting. 

The  weights  are  no  longer  integers  and  a  table 
can  be  constructed  similar  to  the  one  given  for  the  CCD  (refer 
to  Table  10)  . 

The  integrals  for  each  period  for  the  ’x’  numera¬ 
tor  are 

for  mean  angle  0: 

sinu  Uxp^sina  +  dAyp/ScosoTI  (70) 


for  angle  90°: 


cosa 


^4Axp  J  2cosa  -  •'lAyp/2s  i  naj 


for  angle  1 80  : 


-s  i  nu  |dAX(- (  -  )  i  net  +  -lAyn  (  -  )/2eosrvJ 


(  77) 


(  78) 


139 


Ret  icon  nth  element  (x  weights) 


(  7!-*) 


and  for  angle  270°: 

-cosa  4Axr  (  -1  )  t’Jcosa  +  -lAvp»  7s  i  n  x 

and,  the  sum  is 

8/?Axp/2  (80) 

which  is  independent  of  a. 

A  similar  calculation  for  the  denominator  yields 
for  mean  angle  o: 

2 

sin“a  •  2wp  •  S’!  (81) 

for  tt/2 

2 

cos  a  •  2wp  •  i”!  (  82 ) 

for  tt 

9 

-sin  a*  2wp  *  /2(-l)  (82) 

and  3tt  /  2 

-cos  a  •  2wp(-l)/!T  (81) 

with  the  sum 

4/2wp  / 2tt  (,8.r. ) 


Adding  these  to  the  DC  terms,  we  obtain  exact ]y 
the  same  result  that  was  obtained  with  the  CCD  algorithm. 

In  conclusion,  the  use  of  a  Retic.on  with  sequen¬ 
tial  readout  will  not  affect  the  results  if  the  appropriate 
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weighting  functions  are  applied  for  each  diode.  What  this  implies 
is  a  look-up  table  and  two  multipliers  before  division,  Derota¬ 
tion  is  not  required. 

Finally,  since  the  readout  angle  or  axis  rotation 
angles  do  not  appear  in  either  numerator  or  denominator,  the 
variance  of  measurement  will  not  be  affected  by  readout  angle. 

0.4  Input  to  Arrays:  Techniques 

With  the  assurance  obtained  from  the  preceding  dis¬ 
cussions  that  buried  channel  two-dimensional  arrays,  linear  CCD’s 
and  self-scanned  diode  arrays  can  be  used,  methods  of  efficiently 
introducing  light  to  the  sensitive  elements  must  be  found.  Some 
of  these  are  shown  in  Figure  30.  The  buried  channel  CCD  has  dead 
columns  and  is  only  47%  efficient.  Two  strategies  are  possible. 
One,  shown  here,  is  to  use  cylindrical  condensing  elements, 
either  directly  to  the  array  or  to  a  cathode  which  is  then  imaged 
on  the  array.  The  cylindrical  lens  elements  are  only  30  micro¬ 
meters  wide  and  would  need  to  be  made  by  deposition  or  by  some 
photoresist  technique.  If  the  window  to  a  cathode  were  1  mm 
thick,  a  SiO£  deposit  only  0.16  micrometers  thick  would  be  suf¬ 
ficient  to  form  a  lens. 

The  other  technique  would  be  to  create  big  photons  by 
diffusion  of  an  amplified  signal  so  that  its  width  would  be  equal 
to  the  detector  subaperture  area.  This  would  occur,  for  example, 
in  photon  coupling  because  of  loss  of  resolution  in  the  phosphor. 
The  result  would  be  an  effective  overlap  of  subapertures  which, 
if  matched  to  a  deformable  mirror  influence  function,  would  be 
an  advantage*  in  the  accuracy  of  wavefront  reconstruction.  The 
disadvantage  would  be  an  increase  in  the  width  of  the  pulse 
height  distribution  and  consequently  an  increase  in  the  noise 
factor  of  10  to  25%.  The  problems  will  be  discussed  further 
under  photon  coupled  systems. 
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The  coupling  of  the  linear  arrays  to  a  two-dimensional 
subaperture  area  is  more  difficult.  Two  options,  both  of  which 
can  be  handled  by  fiber  optics  manufacturers,  are  shown  in 
Figure  30b.  The  first  of  these  is  a  drawn  fiber  bundle  made  of 
face  plate  material  with  7  micrometer  fibers  and  with  a  high 
transmission.  The  only  problem  is  that  no  one  will  use  the  fine 
fused  materials  which  would  be  much  smaller  than  the  subaperture 
size.  The  coarser  materials  that  manufacturers  will  use  have  50 
micrometers  or  greater  diameter  cores  and  substantial  cladding, 
at  least  equal  to  the  core.  This  would  create  a  very  ineffi¬ 
cient  bundle.  The  solution  is  to  go  one  step  further;  separate 
the  fibers  by  a  millimeter  and  provide  a  16x16  mm  lenticular 
lens  array  to  input  the  subaperture  field  to  the  fiber.  This 
can  be  done  efficiently.  The  author  has  made  one  such  array  by 
milling  spherical  indentations  in  copper  and  then  casting  a  posi¬ 
tive  in  plastic.  The  result  was  a  precision  array  of  foci  which 
were  nearly  diffraction  limited  at  f/10. 

The  other  end  of  the  fibers  would  be  arranged  in  a 
coherent  line  25.6  mm  long  and  the  assembly  encapsulated. 

Galileo  has  given  us  a  budgetary  quote  for  an  assembly  of  $250 
each,  plus  a  tool ing  charge  of  $250. 

We  feel  that  the  area  to  linear  conversion  is  tech¬ 
nically  reasonable  and  that  it  can  be  done  as  efficiently  as 
anti-reflection  coatings  will  permit. 

The  output  spacing  would  be  at  least  100  micrometers. 

6.5  Electron  Bombardment  Mode 

The  development  of  devices  which  utilize  the  electron 
hole  pairs  generated  by  absorbing  a  high  speed  electron  began 
with  beta  particle  silicon  and  germanium  detectors  in  the  early 
sixties.  The  idea  of  using  the  method  for  detecting  photo 
electrons  inside  an  image  intensifior  originated  with  Beaver  and 


McLlwain  of  U.C.,  San  Diego,  and  Choissoir  of  Electric  Vision  in 

1970 .  * 

The  first  device  consisted  of  a  cesium  antimonide 
cathode  at  one  end  of  a  vacuum  tube  and  a  38-element  diode  array 
at  the  other.  Focusing  was  magnetic.  A  potential  difference  of 
20  kV  was  applied  to  the  tube.  Since  an  electron  hole  pair  is 
generated  for  every  3.66  eV  loss  in  the  material,  the  potential 
number  of  pairs  generated  in  the  semi-conductor  will  be  5700  if 
the  20  kV  is  completely  absorbed  in  the  depletion  region.  Deduc¬ 
ting  a  3  kV  loss  in  a  1.7  micrometer  thick  dead  layer  on  the  diode 
the  tube  had  a  gain  of  4800.  Since  the  tube  had  a  noise  component 
equivalent  of  7  keV ,  single  photon  initiated  events  could  be 
detected . 

The  pulse  height  distribution  for  such  a  device  is 

quite  narrow,  and  the  noise  factor  will  be  low  --  about  F  =  1.016, 

because  the  variance  in  the  number  of  photoelectrons  generated 

is  E/E  =  4800. 
o 

Since  the  first  tube,  a  large  variety  of  configurations 
has  been  built,  mostly  with  arrays  of  discrete  diodes  -  up  to 
500  in  the  case  of  the  space  telescope  detector. 

There  has  been  a  large  number  of  attempts  to  introduce 
scanned  arrays  into  intensif iers .  So  far,  however,  there  seems 
to  have  been  only  one  successful  solution.  The  structure  of 
silicon  array  devices  is  complicated  and  includes  conductors, 
insulating  layers  and  active  elements.  The  diodes  in  the  early 
tubes  were  relatively  unaffected  by  the  bombardment,  since  there 
was  nothing  in  the  focal  plane  except  diodes  and  conductor  leads 
through  the  header. 


"A  Digital  Multichannel  Photometer”,  Edward  Beaver  and  Carl 
McLlwain,  Rev .  Sci .  Inst  .  ,  42,  1321. 
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The  complicated  devices  contain  switches,  electrodes, 
clock  circuitry,  etc.,  all  of  which  would  be  permanently  affected 
by  the  creation  of  electron  hole  pairs  inside  them. 

The  first  problem  that  was  solved  was  to  find  a  way  to 
prevent  contamination  of  the  silicon  devices  by  the  alkali  metal 
cathode  material.  This  was  first  successfully  accomplished  by 
Electron  Vision,  and  then  by  a  variety  of  firms,  including  ITT, 
The  methods  are  described  in  their  literature. 

The  most  important  question  of  all  is  whether  the 
arrays  can  survive  direct  electron  bombardment,  In  photon 
imagery,  the  photons  create  hole  pairs  only  in  absorbing  media 
such  as  silicon,  and  travel  through  transparent  insulators  and 
electrodes  without  effect.  The  electron  beams,  on  the  other 
hand,  crash  through  the  materials  leaving  a  trail  of  electron 
hole  pairs  everywhere. 

The  resulting  ionization  becomes  trapped  at  insulated 

interfaces,  causing  increases  in  leakage  currents,  creates  fields 

that  must  be  overcome  by  increasing  gate  voltages  which  must  get 

progressively  greater  until  the  limits  of  the  materials  are 

exceeded.  Experiments  with  buried  channel  devices  which  are 

5  9 

directly  exposed  show  that  at  10  rads /cm“  they  become  unusable. 

This  value,  translated  into  the  compensated  imaging  exposure 

levels  of  100  photons/c.ycle/subaperture,  would  be  reached  in  1.8 

minutes.  Most  of  this  failure  is  incurred  in  the  electrode 

transport  line  structure.  Things  may  be  better  if  these  channel 

are  protected  by  very  heavy  aluminum  barriers  10  micromotors  or 

more  thick.  The  life  in  this  case'  could  be  extended  to  3  hours. 

Because  of  this,  considerable  <‘ffort  has  been  expended 

to  create'  thinned  CTO's  which  can  bo  exposed  through  the  base. 

The  expectation  is  that  the  exposure  will  be  confined  to  the 

silicon  depletion  layer  below  the'  electrodes.  Leakage  currents 

still  increase,  however,  and  measurements  show  that  leakage 

<) 

reaches  (J>  of  saturation  at  1  x  10'  el  oc  t  rons/p  i  xe  1  . 
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not  bad  and  is  about  10  better  than  the  front  side  illumination. 

But,  the  device  in  the  CIS  environment  would  be  completely  dead 
4 

in  1.8  x  10  minutes  and  would  be  unusable  for  weak  targets 
(one  photon/cell)  after  three  hours  of  the  moderate  exposure. 

The  only  devices  that  have  survived  the  radiation  are 
the  discrete  diode  and  the  shielded  linear  Reticon  self-scanned 
diode  arrays.  This  is  because  there  are  no  electrodes  —  just  a 
conductor  --  to  lead  the  charge  out  to  the  protected  gates. 

The  reason  this  device  would  work  and  a  shielded  CCD 
would  not  is  due  to  the  structure  of  each  photosite,  diagrammed 
in  Figure  44 . 

In  the  Reticon,  there  is  a  simple  pin  junction  forming 
a  diode  which  is  protected  by  a  1  micrometer  thick  layer  of  SiO^. 
Above  this  is  a  mechanically  separate  metal  shield  which  is 
easily  aligned  because  the  diode  widths  are  425  to  600  micro¬ 
meters  . 

The  diffusion  region,  in  which  electron  hole  pairs  can 
be  created  and  whose  charge  will  be  accrued  to  the  junction,  is 
50  micrometers  deep.  This  permits  the  use  of  up  to  77  kV 
electrons.  The  Si02  layer  on  top  is  normally  3  micrometers, 
but  Reticon  will  process  the  arrays  to  provide  the  thinner  layer 
for  electron  detection. 

In  the  lateral  directions,  the  12.5  micrometer  diodes 
are  separated  by  N  type  silicon  so  that  electrons  which  are 
created  in  this  space  between  will  flow  to  both  junctions, 
dividing  the  current  between  them  in  proportion  to  the  distance 
to  them.  This  creates  an  overlapping  sensitivity  between  the 
diodes,  as  shown. 

This  spreading  can  provide  an  advantage  if  it  matches 
the  influence  function  of  a  preceding  adaptive  optical  element. 

!'n  fort  unat  e  1  y  .  the  spreading  is  one-dimensional  and  continuity 
from  element  to  element  at  the  edges  of  the  input  array  requires 
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that  the  fiber  arrangement  be  uncoiled  to  a  line  from  the  square 
in  a  serpentine  manner.  This  organization  is  exactly  what  hap¬ 
pens  when  a  flat  piece  of  paper  (like  a  wide  ribbon  of  fibers) 
is  folded  to  form  a  fan. 

The  spreading  unfortunately  will  increase  the  width  of 
the  pulse  height  distribution  and,  for  this  reason,  should  be 
prevented.  The  large  fiber  area  linear  converter  with  50  micro¬ 
meter  cores  spaced  at  100-125  micrometer  intervals  will  be 
reduced  4  to  5  times  onto  the  array.  The  resolution  of  the 
electron  optics  will  cause  this  image  to  spread.  Fortunately, 
the  resolution  capability  of  electron  optics  is  very  good,  par¬ 
ticularly  in  the  small  6  mm  long  fields  considered  here,  so  that 
the  image  spots  should  not  spread  appreciably.  (Poor  performance 
of  intensifiers  is  due  to  the  phosphors  as  much  as  anything  else.) 
Cross  talk  and  pulse  width  increases  should  be  minimized. 

The  linear  CCD's  might  be  treated  in  the  same  way  as 
the  Reticon,  using  a  thick  external  shield  to  protect  the  trans¬ 
fer  lines.  The  structure  of  the  CCD  which  is  plowed  by  the  high 
energy  electrons,  in  the  second  sketch  of  Figure  44,  is  composed 
of  a  sequence  of  layers  starting  with  a  1  micrometer  Si02  insul¬ 
ator,  then  a  0.3  micrometer  polysilicon  layer,  a  0.05  micrometer 
silicon  nitride  electrode,  another  SiC>2  insulator,  and  then  the 
depleted  region  in  which  the  usable  electron  hole  pairs  are 
generated . 

This  last  layer  is  about  7  micrometers  deep,  limiting 
the  gain  to  7800  and  the  optimum  electron  beam  energy  to  32  kV, 

The  failure  mechanism  is  surface  trapping  of  carriers 
between  the  layers. 

The  first  of  the  sel f -scanned  linear  devices  was 
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obtained  by  Tull  and  Nathan  at  the  University  of  Texas  in  1974 
and  is  still  operating  in  the  photon  counting  mode. 

The  efficiency  of  the  two  designs  are  basically  dic¬ 
tated  by  the  photocathode  which  will  be  20%  to  30%  in  our  case 
and  40%  in  the  other.  The  difference  in  optical  bandwidth  will 
make  the  opaque  photocathode  three  times  as  efficient  as  the 
other  design.  This  is  worth  the  cost  of  providing  a  complex 
input  optical  lens. 


* 

"The  Self  Scanned  Digicon,"  Tull,  Choisser, 
Optics ,  14 ,  #5,  May  1975. 
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G • 6  Photon  Coupled  Systems 

A  variety  of  image  intensifier  tubes  and  optically 
coupled  solid  state  image  device  configurations  are  shown  in 
Figures  32a-e .  The  basic  device,  Figure  32;l  >  consists  of  a  sim¬ 
ple  intensifier  with  a  fiber  optic  plate  on  the  output.  This 
is  coupled  directly  to  the  solid  state  array  by  a  fiber  optic 
plug.  The  array  is  assumed  to  be  a  buried  channel  device  be¬ 
cause  that  is  the  only  configuration  that  can  be  made  to  operate 
like  a  small  array  at  high  speed  by  the  technique  described  in 
Section  6. 4. 

The  resolution  of  intensifiers  is  often  tin-  result  of 
a  trade-off.  The  efficiency  of  the  phosphor  increases  with 
thickness  because  higher  energy  electrons  can  be  used.  On  the 
other  hand,  the  resolution  decreases  with  thickness  because  of 
multiple  scatter.  The  relatively  low  resolution  that  results 
can  be  used  to  advantage  in  the  AWS  as  a  means  of  creating  photon 
initiated  blobs  which  can  be  large  enough  to  bridge  over  the 
opaque  areas  above  the  CCD  transport  registers.  In  consequence, 
if  the  resolution  is  matched  to  the  CCD  subaperture  cell  size, 
at  least  50%  of  the  photon  induced  signal  will  be  detected.  In 
consequence,  the  quantum  efficiency  will  not  be  reduced  but  the 
noise  factor  due  to  a  wide  pulse  height  distribution  will  in¬ 
crease  . 

The  size  of  the  Fairchild  CCD  elements  in  the  190  x 
2'14  array  are  18  x  14  micrometers.  The  added  pair  plus  the 
opaque  region  yields  a  subaperture  size  of  36  x  30  micrometers. 
The  resolution  of  typical  image  intensifiers,  ranging  from  20  to 
50  1 /mm ,  is  compatible  with  this  subaperture  area. 

The  fraction  of  photons  which  is  emitted  by  a  phosphor 
and  which  emerges  from  the  tubes  is  quite  small  because  of  frus¬ 
trated  internal  reflection  at  the  final  glass/air  surface.  Loss¬ 
es  from  the  electron  beam  side  are  prevented  by  thin  aluminum 
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coating  on  the  phosphor.  If  n  is  the  index  of  the  phosphor,  the 
fraction  of  photons  that  can  escape  at  all  angles  from  an  alumi¬ 
nized  tube  is  l/n^. 

Since  the  index  of  refraction  for  most  phosphors  is 
quite  high,  only  a  small  fraction  can  escape.  This  is  avoided 
by  fabricating  the  phosphor  as  a  powder.  The  particles  then 
scatter  the  emitted  photons  until  they  are  scattered  into  the 
face  plate.  As  a  result,  up  to  45%  of  the  light  in  an  aluminized 
tube  can  be  emitted  over  a  full  hemisphere.  If  this  output  were 
coupled  to  the  array  by  f/1  optics,  the  total  efficiency  would 
only  be  4.5%.  The  use  of  fiber  optics  to  optically  couple  the 
arrays  to  the  tube  and  avoid  internal  reflection  losses  increases 
the  efficiency  dramatically.  The  face  plate  must  be  an  Extra 
Mural  Absorption  (EMA)  type  (fibers  imbedded  in  a  black  matrix) 
to  preserve  resolution;  the  transmittance  of  a  D-14  type  with  6 
micrometer  fibers  from  a  Lambertian  source  is  45%.  The  spectral 
transmittance  to  the  array  cuts  off  at  0.4  micrometers  (u-1 2'-') 
and  is  gone  at  \  =  0,35  micrometer.  This  is  shown  in  Figure  45. 

The  simple  structures  of  the  photon  coupled  arrays  are 
then  fairly  well  represented  by  Figuvs  32a  and  32b.  In  the-  lir^t 
described  above,  there  is  a  simple  cathode,  then  a  fiber  face 
plate  to  an  area  array.  The  second  is  identical  but  uses  a  fi¬ 
ber  area  linear  transformer  to  a  1  inear  array.  Tin-  area  linear 
transformer  must  be  made  of  very  fine  G  micrometer  or  less  li¬ 
bers,  however.  Manufacturers  that  were  contacted  did  not  feel 
that  this  was  reasonable.  The  feasible  area  linear  t runs  termer 
would  be  made  of  large  fibers  (50  micrometer  core)  w  i  ' h  sub¬ 
stantial  cladding  creating  100  micrometers  center  to  center 
spacing.  Such  a  largo  bundle  is  incompatible  with  direct  con 
nection  to  an  array.  The  third  design  with  the  transformer  up 
front  with  a  lenticular  lens  array  for  efficient  input,  is  pro- 


ferred.  The  intensifier  electron  optics  provides  the  necessary 
change  of  scale  to  the  small  linear  array. 

None  of  these  designs  provide  much  gain;  however,  30  to 
50  luminous  gains  are  typical  for  single  stage  tubes;  conse¬ 
quently  second  generation  intensifiers  or  cascaded  tubes  in  Fig¬ 
ures  32c  or  32d  are  desirable. 

The  basic  problem  of  the  second  generation  tube  is  the 
fact  that  the  microchannel  plate  has  a  fairly  coarse  structure 
which  has  opaque  regions  and  only  65^  of  the  primary  photoelec¬ 
trons  get  amplified.  This  direct  loss  of  quantum  efficiency 
should  be  avoided,  and  thus  a  cascade  arrangement  with  a  low 
gain  Gen  I  tube  preceding  the  Gen  II  tube  to  provide  a  fat  pho¬ 
ton  pulse  is  desirable.  There  are  severe  problems  in  the  use  of 
microchannel  plates  which  will  be  reviewed  in  the  next  section, 
and  this  use  is  not  recommended. 

The  preferred  photon  coupled  design  uses  an  array,  and 
two  stages  of  gain.  The  low  resolution  precludes  the  use  of  de¬ 
tectors  with  small  elements,  30  to  40  micrometers  being  as  small 
as  pract ieal . 

The  last  question  to  be  studied  is  whether  phosphor 
decay  affects  the  measurements  of  tilt  ir  a  nutating  image. 

In  the  case  of  time  delayed  sampling,  which  is  a  con¬ 
sequence  of  the  use  of  self  scanned  diode  arrays,  the  resulting 
effect  was  equivalent  to  a  rotation  of  the  coordinate"  system. 

This  deviation  was  corrected  by  the  use  of  a  weighting  function 
that  included  the  rotation  of  time  delay  angle". 

There  is  a  difference  between  time"  delay  and  phosphor 
decay  because  a  delay  is  well  defined,  while  decay  introduces  a 
smear . 

The  effect  on  the  signal  can  be  derive'd  simply  by  us- 
-  ing  the  square  uniform  object,  and  by  separating  AC  and  DC  terms 


in  numerator  and  denominator,  as  before.  Assuming  a  simple  de¬ 
cay  so  that  the  output  of  the  processor  using  the  simple  integra¬ 
tor  weighting  functions  is  given  by  the  convolution 
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where  a  is  the  decay  constant.  The  convolution  yields  the  I  sum 
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If  this  is  evaluated  by  the  I  weighting  and  applying 
the  appropriate  integrals 
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for  the  x  deflection  we  obtain 
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The  coordinate  axis  rotation  is  given  by  the  ratio  of 
nutation  frequency  to  decay  frequency. 

The  denominators  for  ,'e  and  /s  weightings  become 
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The  values  for  Ax  and  Ay  are  given  by 
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The  variance  on  the  measurement  of  tilt  is 
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where  N  is  the  total  number  of  photons.  In  terms  of  N, 
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The  consequences  of  the  smear  caused  by  decay  is  an  in¬ 
crease  in  the  variance.  It  simply  means  that  the  contrast  of  the 
signal  is  lower.  There  is  one  problem  here;  the  persistence  of 
cathode-luminescence  is  not  generally  simple.  The  decay  curves 
for  these  fast  phosphors  are  shown  in  Figure  46.  These  are  the 
P16-P15  and  P24.  All  of  these  are  best  described  by  two  time 
constants  so  that  the  intensity  is  described  by 

I  =  Ae-  tl  +  (1  -  A)e_t2t  (98) 


1  56 


For  the  P16, 


So  that 


5  =  0.028  ps  and  T2  5  =  °-09  gs 


=  24.8  x  106  and  a2  =  0.77  x  106 


and  A  -  0.75 


(  100) 


The  P15  and  P24  are  both  zinc  oxide.  The  P15  has  a  UV  compo¬ 
nents,  but  the  P24  does  not  and  emits  entirely  in  the  blue  green. 
Its  time  constants  are 


a ^  =  6.93  x  10 


a g  =  0.77  x  10 
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A  -  0.75 


Most  of  the  emitted  energy  is  in  the  second  term,  since 
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'totalg  ^  1  A)/a2  0.33 


The  magnitude  of  the  rotation  term  for  the  P24  —  for  10  kHz  nu¬ 
tation  is  0.082.  This  is  small  enough  so  that  a  modification  to 
the  weight  ng  function  could  be  used  to  remove  the  mean  rotation 


angle;.  The  fact  that  the  proportion  of  light  from  each  decay 
mode  is  a  function  of  temperature:-  and  Learn  intensity  is  not 
crucial  because  the  rotation  angle  is  so  small. 


The  PI  ft  decay  is  faster  but  half  the  L'Y  light  will  be 
lost  in  the  fiber  optic  plug,  This  is  illustrated  in  Figure  do 
where  the  relative  brightness  of  the  P16  and  P24  phosphors  are 
shown , 

The  second  curve  for  the  P16  shows  the  effect  of  ab¬ 
sorption  in  the  fiber  optic  face  plate.  The  relative  losses  in 
the  UV  at  the  CCD  further  reduce  the  value  of  the  P16  output. 

In  consequence,  the  P2-4  is  preferred,  in  a  two  tube  cascade, 
coupled  to  an  area  CCD, 

ft , 7  Electron  Coupled  Mode 

The  development  of  the  mi crochannel  plate,  which  is  a 
honeycomb  of  electron  multiplier  tubes,  increased  the  gain  of 
image  intensifiers  immensely  and  permitted  the  development  of 
adequate  single  stage'  and  proximity  focus  devices. 

The  device  provides  our  •vn  t  rather  than  voltage  gain. 
Thus,  a  sensor  can  simply  collect  tin;  current  rather  than  provide 
the  means  of  converting  the  high  voltage  to  current  and  then  col¬ 
lecting  it  as  well.  The  configurations  are  shown  in  Figures  dd 
,nd  B.  The  first  of  these  consist  s  of  a  parallel  array  of 
anodes  in  a  proximity  focus  mi crochannel  plate  intensi f ier . 

This  tube  is  manufactured  by  ITT  and  is  probably  the  onlv  de¬ 
vice  that  can  be  obtained  without  any  more  special  offer’  than 
changing  the  pattern  of  anodes.  In  this  device,  the  anodes 
would  be  0.9  mm  square  plates  with  0.1  mm  separation.  l'ach  anode 
is  led  through  the  ceramic  header  to  a  pin. 

The  electronics  for  each  anode  would  consist  of  a  1(H) 
MHz  amplifier  discriminator  and  a  counter.  The  counters  would 
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be  buffered  and  read  into  the  processor  in  serial  stream. 

The  second  version  uses  a  self  scanned  diode  array  in 
which  each  diode  is  covered  by  an  aluminum  plate.  This  would  be 
a  Ret  icon  device  similar  to  the  SSANACON^  that  flew  on  Voyager  as 
a  part  of  a  UV  spectrometer. 

The  anodes  in  this  second  device  are  very  small,  0.1  mm 
square  with  0.01  mm  gaps  between  them. 

The  final  electrons  out  of  the  MCP  with  about  100  eV 
have  too  little  energy  to  penetrate  this  device  and  there  is 
little  possibility  of  damage. 

5  7 

The  gain  that  can  be  obtained  is  high,  10  to  10  ,  and 
there  is  little  question  about  the  value  as  a  photon  counting  de¬ 
tector.  The  possibility  of  obtaining  an  array  was  explored  with 
EG&G  Reticon,  and  they  proposed  a  16  x  16  element  device  using 
16  parallel  channels  of  16  anodes  in  each.  The  channels  would  be 
read  at  1  MHz  in  parallel,  providing  a  16  MHz  overall  data  rate. 
The  budgetary  estimate  for  a  dozen  of  these  arrays  with  nonrecur¬ 
ring  costs  included  would  be  100K  and  with  delivery  9  months 
after  receipt  of  order. 

Closer  analysis  of  the  characteristics  of  microchannel 
plates  revealed  some  critical  problems.  These  are 

1.  Life 

2.  Slow  response 

A  microchannel  plate  consists  of  a  honeycomb  glass 
structure  with  a  secondary  emissive  coating  on  the  channel  sur¬ 
faces.  The  structure  is  formed  by  three  glasses  with  different 
viscosi ty . 

During  the  drawing  and  fusing  steps  the  outer  glass,  or 


*  The  SSANACON  -  a  self  scanned  anode  array  with  a  mocrochan- 
no  1  plate  electron  multiplier.  A,  l.  Uroauioot,  B.  K.  Handel, 
APPLIED  OPTICS. 
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cladding  glass,  glows  in  such  a  way  that  external  force  and  sur¬ 
face  tension  deform  the  glass  matrix  to  the  final  hexagonal 
shape.  The  process  of  the  different  manufacturers  leads  to 
slightly  different  shapes  and  sizes. 

The  Galileo  MCP ' s  have  12  to  38  pm  diameters  and  open 
area  ratios  of  55%.  ITT’s  MCP's  have  channels  12  ym  in  diameter 
and  an  open  area  ratio  of  63%.  Only  electrons  striking  the  open 
area  are  multiplied;  thus,  if  no  other  factors  modify  this,  the 
quantum  efficiency  is  55%  (or  63%)  of  the  cathode  efficiency. 
Secondary  electron  emission  from  the  web  faces  can  increase  this 
efficiency  up  to  85%. 

One  of  the  problems  with  the  structure  is  that  the  pos¬ 
itive  ions  can  be  accelerated  in  the  opposite  direction.  These 
will  poison  the  cathode  if  they  are  permitted  to  leave  the  chan¬ 
nel  plate. 

2 

The  life  end  point  will  be  between  0.01  to  0.1  C/cm  . 
Under  CIS  conditions,  with  10°  photons/sec  and  gains  of  10°  for 
photon  counting  with  1  mm  anodes,  this  would  correspond  to  a  life 
range  of  1.5  hours  to  15  hours.  In  the  Reticon  anode  array,  the 

4 

lower  noise  might  permit  one  to  run  at  a  10  gain,  but  the  life 
would  only  be  10  minutes  to  100  minutes. 

To  cure  this,  a  thin  film  is  placed  over  the  input  end 
of  the  channel  plate  which  prevents  the  ions  from  escaping,  but 
permits  the  higher  energy  electrons  to  pass  through.  Unfortu¬ 
nately,  this  decreases  the  secondary  electron  flux,  reducing  the 
quantum  efficiency.  The  life  is  extended  by  a  factor  of  about 
1000.  This  is  marginally  acceptable  for  the  large  anode  array, 
but  not  for  the  small  array. 

The  second  problem,  slow  response,  is  even  worse.  The 
microchannel  plate  behaves  just  like  a  photomultiplier,  which 
uses  discrete  dynodes.  The  potential  on  each  dynode  is  usually 


provided  by  a  resister  chain,  In  such  tubes,  the  resistance  is 
chosen  to  provide  sufficient  current  to  the  dynodes  to  handle  the 
signal  level.  In  the  MCP,  the  same  situation  exists,  but  the  re¬ 
sistance  chain  is  the  channel  wall  itself. 

14 

This  impedance  is  high,  about  5  x  10  ohms/channe] . 

In  a  single  plate  intensifier  with  1000  V  across  it  and  a  gain 

5 

requirement  of  10  ,  the  number  of  pulses  that  could  be  generated 

would  be  only  125/sec,  The  effective  dead  time  is  8  millisecs. 

5  2 

There  are  a  lot  of  channels,  however,  5  x  10  per  cm  , 

5 

so  that  the  maximum  flux  at  10  gain  for  the  1  mm  anode  is  6.3 

5 

x  10  ,  63  events  per  nutation  cycle,  or  S/N  £8.  With  the  small- 

4 

er  anode,  even  if  lower  gain  of  10  is  usable,  the  smaller  area 
limits  the  flux  to  6.3  x  104  events,  or  S/N  £2.5, 

Operation  near  saturation  is  not  reasonable  because  the 
responses  will  be  nonlinear.  The  noise  factors  will  go  up,  and 
linearity  will  be  lost. 

The  MCP  is  thus  seen  to  be  less  efficient  than  other 
schemes,  and  cannot  take  the  flux  levels  anticipated  for  compen- 

4 

sated  imaging.  These  levels,  varying  from  10  events/sec  to  5  x 

7 

10  events/sec,  will  generally  overwhelm  the  electron  coupled  de¬ 
tector  using  microchannel  plates. 

The  only  electron  coupled  mode  that  will  work  and,  in¬ 
deed,  the  only  tried  sensing  mechanism  that  can  be  purchased  off- 
the-shelf,  is  a  photomultiplier  array.  This  is  not  recommended 
for  the  full  256  channel  system,  but  it  is  recommended  for  a  21 
channel  breadboard  system. 

The  basic  arrangement  is  shown  in  Figure  34,  A  digital 
photon  counter  will  be  used  instead  of  an  analog  integrator.  The 
amplifier  discriminator  will  be  a  single  integrated  circuit  chip. 
This  is  followed  by  a  12  bit  counter  and  a  12  bit  buffer.  This 
parallel  arrangement  will  be  read  sequentially  into  the  pipeline 


processor . 


The 

lower  quantum 


only  drawback  in  the  use  of  PMT’s  is  the  slightly 
efficiency  that  is  a  result  of  mass  production. 


t 


6.8  Recommended  Detector  Design 

There  are  three  recommended  designs;  two  designs  use 
electron  bombardment  and  differ  only  in  the  nature  of  the  optical 
system.  The  third  design  is  a  photon  coupled  arrangement. 

In  addition,  the  recommended  experimental  system  for 
twenty-one  channels  uses  photon  multipliers  because  it  presents 
the  least  expensive  way  to  prove  the  performance  of  the  sensor 
and  processing  techniques. 

The  three  designs  are  shown  in  Figures  47,  4!*  and  50. 

The  most  efficient  photocathodes  that  have  been  devel¬ 
oped  are  opaque.  Normally,  this  has  prevented  their  use  in  image 
tubes  because  it  is  difficult  to  arrange  for  parallel  conjugate 
planes,  although  one  device  with  a  Schmidt  telescope  objective 
inside  an  intensifier  was  used  in  a  UV  spectrometer  on  an  Apollo 
mission.  The  occluded  center  image  plane  was  the  cathode. 

The  fact  that  the  final  image  on  the  cathode  is  linear 
because  the  detector  arrays  are  linear  permits  simplifications, 
since  the  arrangement  need  only  provide  for  parallel  conjugate 
lines.  The  cathode  can  then  be  narrow  and  tilted  at  an  angle  so 
that  the  optical  axis  and  electron  imaging  axis  can  be  at  right 
angles . 


b , 8 , 1  Design  1  Internal  Optics  KB  Mode.  The  first 
and  preferred  design  uses  an  elliptical  mirror  inside  of  the 
image  intensifier  to  collect  the  photons  from  the  fiber  optics 
linear  array  and  to  image  these  back  on  the  closely  adjacent 
photocathode.  The  optical  element  for  a  very  large  numerical 
aperture  may  be  an  ellipsoid  of  revolution  about  an  axis  between 


the  line  of  fiber  optic  and  the  GaAs  cathode. 

The  value  of  this  design  can  lie  judged  by  t  lie  ef¬ 


ficiency  of  the  optical  system  to  collect  all  of  the  photons 


from  the  fibers  and  reimage  these  on  the  cat  bode . 


Thus,  in  order 


<>xt  I'rnal  ont  i  rs 


to  evaluate  the  design,  we  must  obtain  an  estimate  of  the  numeri¬ 
cal  aperture  the  optics  can  handle. 

The  calculation  for  the  focal  plane  assumes  a 
curved  field,  and  for  simplicity,  a  spherical  element.  If  h  is 
the  distance  off-axis  of  the  source  point,  displaced  toward  the 
mirror  by-a  distance  o,  and  if  h'  is  the  image  distance  or  ray 
intercept  from  the  axis  in  the  plane  normal  to  the  axis  including 
the  source  for  unit  magnification  and  if  9  is  the  angle  of  a  ray 
from  the  source  relative  to  the  axis,  the  relation  for  h  and  h' 
is 


(h  +  h'  ) 


— 

cos-? 


+  sin 


hcos9 


(  103) 


=  2hcos9  +  25sin9 


At  minimum  aberration,  h  +  h’  -  2h ,  so  that  one 
obtains  the  curvature  of  the  field 

2 

o  =  h  / R  (104) 

where  R  is  the  radius  of  the  sphere.  The  radius  of  the  focal 
surface  is  Rf  =  -R/2,  which  is  typical  of  spherical  mirror  sys- 
t  eras . 

With  this  field  curvature  assumed,  the  aberration 
for  a  spherical  mirror  is  given  by 


Ah 


(2h“/R) 


t  a  n  °  (  1 


-  -sinn)  -  sin?) 


(  I  Of. ) 


The  evaluation  is  given  in  Figure  )g.  If  the  al¬ 
lowed  deviation  is  not  to  exceed  ±  50  micrometers,  which  is  ha  1 f 
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SPHERE 


Ahor  rat  i  nns  of  spherical  and  ellipsoid  mirrors 
at  1:1  magnification.  Field  ±  12.8  mm, 
radius  =  100  mm. 
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the  spacing  between  fibers,  then  the  numerical  apertures  at  the 
extreme  elements  h  =  12.5  mm  and  R  =  100  mm  is  limited  to  NA  = 
0.34. 

At  other  field  positions  closer  to  the  center,  the 
numerical  aperture  increases  so  that  at  h  =  8.3,  NA  =  0.40.  At 
h  =  6.25,  NA  =  0.47;  at  h  =  3.125,  NA  =  0.69;  and  at  the  center, 
of  course,  NA  =  1 ,  since  the  point  is  at  the  center  of  curvature. 
This  point  is  not  achieved,  however,  because  the  fibers  and 
cathode  are  separated  by  a  few  millimeters  so  that  the  greatest 
numerical  aperture  w-ill  probably  be  about  0.7. 

Obviously,  better  performance  can  be  obtained  by 
using  an  ellipsoid  which  provides  two  fully  corrected  field 
points  instead  of  the  one  that  is  provided  by  the  sphere.  If 
these  points  are  located  at  ±  1/4  of  the  field,  then  the  extreme 
point  aberrations  are  given  by  curve  3,  and  the  best  points  by 
curve  4.  Assume,  then,  the  spherical  case  gives  a  good  estimate 
of  small  deviation  errors  about  a  design  point.  The  overall  av¬ 
erage  NA  will  probably  be  about  0.60. 

This  is  completely  compatible  with  the  optics  at 
the  input  to  the  fiber  optics,  and  with  the  numerical  aperture  of 

single  fiber  elements.  If  the  input  to  the  fibers  is  constrained 

to  f/2,  NA  =  0.24;  if  the  fiber  NA  -  0.5  and  tight  bends  are 
avoided,  all  of  the  light  from  the  fibers,  except  for  that  lost 
in  the  open  center  of  the  ellipsoid  (NA  <_0.06  (a  27  loss)  where 
the  electrons  go  through  the  mirror),  will  reach  the  cathode. 

The  efficiency  of  the  system,  as  long  as  the  nu¬ 
merical  aperture  is  ^0.5,  is  dictated  by  the  antireflection 
coatings  that  can  be  obtained  for  the  wide  bandwidth  0.42  to 
0.85  urn  sensed  by  the  GaAs  cathode,  which  if  fully  utilized  will 
provide  photoelectrons  for  40/.  of  the  solar  photons.  The  wide 

bandwidth  will  prevent  the  use  of  very  efficient  coatings  so  that 

one  can  probably  only  achieve  987  average*  transmission  per  sur- 
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face.  Thus,  with  four  lenticular  lens  surfaces,  two  fiber  sur¬ 
faces,  one  mirror  with  protected  silver  coating,  and  a  2%  loss 
through  the  central  hole,  the  system  transmission  will  be  about 
84%. 

The  properties  of  the  design  are  summarized  in 
Table  11.  The  noise  factor  is  low  because  the  variance  on  the 
gain,  which  is  given  by  the  reciprocal  number  of  photoelectron 
pairs  generated,  is  small.  This  is  unlike  a  PMT,  where  the  vari¬ 
ance  on  gain  is  given  by  the  reciprocal  number  of  electrons  gen¬ 
erated  in  the  first  dynode. 

With  gains  >1000,  the  noise  factor  is  1  + 

( 1  /  /1000)  . 

The  biggest  problem  in  the  design  may  be  the  de¬ 
sign  of  the  electron  optics,  which  must  provide  a  4X  reduction 
onto  the  array  and  a  worst  case  resolution  of  40  1/mm.  This 
problem  has  been  discussed  with  Electron  Vision  and  ITT,  and 
little  concern  over  the  problem  was  expressed--but  whether  they 
realized  the  full  extent  of  the  internal  optics  or  not  is  uncer¬ 
tain.  The  present  effort  has  not  studied  the  electron  imaging, 
and  consequently,  the  design  of  the  electron  optics  must  be  con¬ 
sidered  a  cost  problem. 

The  last  element  in  the  comparative  table  is  a 
relative  net  S/N  value  given  by  >  T ■ n / ( NF )  where  T  is  the  trans¬ 
mission,  n  the  solar  spectral  efficiency  and  (NF)  the  noise  fac¬ 
tor  . 


0,8.2  Design  2  External  optics,  EB  mode.  The  second 
design  listed  in  the  table  avoids  the  problem  of  optics  inside  of 
the  image  tube,  but  introduces  other  problems.  The  obvious  first 
solution  is  to  use  fast  relay  optics  and  this  solution  is  shown 
in  Figure  40.  Kawa  relay  lenses  are  shown  here  because  the  de¬ 
signs  were  available.  The  picture  tells  all;  there  are  .80  op- 


TABLE  1 1 .  RECOMMENDED  DETECTORS 


DES 1  ON : 

1 

2 

3 

4 

MODE 

EB 

EB 

PC 

EC 

OPTICS 

Internal 

External 

None 

None 

AREA-LINEAR 

CONVERTER 

ARRAY 

Re  1 1 con 

1  i  near 
SSDA 

Ret  icon 

SSDA 

Fa i rchl  Id 

CCD 

Paral lei 
detector 

CATHODE 

Opaque 

Ga  As 

Opaque 

GaAs 

Trans . 

GaAs 

Bi- 

alkali 

,'S(  •  UI(  ■  )d- 

Refrac-  Clam- 
t  ors  shel 1 

/S(OE4)B(  '■  'fP 

0.40 

0.40  0.40 

0.31 

0.11 

TRANSMISSION 

0 . 84 

0.55  0.79 

0.99 

0.97 

DETECTOR  RISK 

None 

None  (all 
off  shelf) 

None 

None 

1 NTENSIEI KR 

DESIGN 

DI  EFI  C'LTiTY 

Greatest 

Med i um 

None,  off 
the  shelf 

None 

NOISE  FACTOR 
( FAT  PHOTON 
PROBLEM  k 

PIT.SE  HEIGHT 
DISTRIBUTION) 

1  .03 

1 .03 

1  .  2 

1.19 

COST  PROBLEMS 

Design  of 
elect  run 
opt i rs 

Co  1 ] er t i en 
efficiency 
of  lenses, 
c  h  roma t i c 
alierrat  ions 

Reso 1 u 1 i on 
i f  3  t  ubes 
requ  i  rc'd 

No  IK1 

NET  SNR 

0. 563 

0.455  0.546 

0.4  62 

0.27 

the 


tical  surfaces  in  all,  and  if  each  surface  reflects  2%, 
transmission  would  be  only  55%.  The  design  of  these  lenses  is 
probably  not  compatible  with  the  very  wide  spectral  bandwidths 
used  to  obtain  full  advantage  of  the  GaAs  cathode. 

The  one  advantage,  though,  is  that  the  back 
focal  length  is  great  enough  so  that  the  designs  of  the  relay 
and  the  electron  optics  are  independent. 

Another  possibility  shown  in  Figure  50  is  the 
clamshell  lens.  The  device  is  completely  reflective  and  chro¬ 
matic  aberrations  are  not  present.  The  numerical  aperture  of  a 
design  which  has  the  input  and  output  in  contact  with  the  para¬ 
bolic  surfaces  is  0,94.  If  a  space  is  created  so  that  the  out¬ 
put  can  get  through  the  face  plate  and  be  reflected  onto  the 
cathode  that  is  tilted,  the  numerical  aperture  will  be  consider¬ 
ably  less. 

The  biggest  problem  is  getting  from  the  clamshell 
to  the  cathode.  A  configuration  with  the  clamshell  axis  at  right 
angles  to  the  electron  tube  axis,  with  the  cathode  at  45°,  would 
place  the  large  conductive  surface  at  the  edge  of  the  electron 
trajectory,  and  it  is  doubtful  that  the  arrangement  could  be  made 
to  work.  The  other  arrangement,  as  shown,  would  image  directly 
onto  the  end  of  the  tube,  where  a  mirror  element  would  direct  the 
light  sideways  onto  the  long  and  narrow  photocathode.  In  order 
to  work,  the  focused  line  on  the  cathode  will  need  to  be  flat. 

In  consequence,  it  may  be  difficult  to  design  the  electron  op¬ 
tics.  Finally,  the  last  right  angled  bend  will  restrict  the  nu¬ 
merical  aperture  in  one  dimension.  This  restriction,  as  seen  in 
the  insert  to  Figure  50,  can  reduce  the  numerical  aperture  to 
0.6.  In  consequence,  although  the  clamshell  has  the  capability 
of  NA  =  0.94,  the  realization  will  be  no  better  than  the  inter¬ 
nal  ellipsoid  configuration. 
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If  this  were  not  enough  of  a  problem,  the  making 
of  the  cathode  and  the  application  of  the  cesium  material  on  the 
GaAs  with  the  mirror  present,  masking  the  areas  that  need  to  be 
clear,  may  pose  practical  difficulties.  For  these  reasons,  the 
clamshell  is  less  attractive  than  the  internal  mirror  design. 

6.8.3  Photon  coupled  system.  The  third  choice,  the 
photon  coupled  system,  has  one  major  advantage:  it  does  not  use 
an  area  linear  converter  and  thus  has  no  loss  elements  before 
the  cathode,  except  for  the  tube  face  itself.  The  photocathode, 
however,  must  be  transparent  and  will  consequently  be  less  sen¬ 
sitive,  providing  only  31£  of  the  possible  number  of  photoelec¬ 
trons.  Thus,  in  Table  11  the  expected  flux  number  is  0.31  • 

0.39  -  0.307,  which  must  be  compared  to  the  0.4  •  0.84  =  0.336 
value  obtained  from  the  best  EB  system. 

The  big  differences  between  the  PC  and  EB  systems 
is  the  gain  mechanism.  The  PC  system,  with  a  gain  of  about  25 
per  stage,  will  require  three  cascaded  image  tubes  to  produce  a 
quantum  limited  signal.  With  this  low  gain  in  the  first  stage, 
the  noise  factor  becomes  significant  and  will  be  at  least  1.2. 

In  consequence,  the  net  SNR  is  the  least  of  all  of  the  recom¬ 
mended  systems,  0.462  vs.  0.563.  Except  for  the  problem  of  reso¬ 
lution,  the  system  has  few  problems,  and  can  be  put  together  en- 
>  rely  using  off-the-shelf  devices. 

In  the  last  column  of  Table  11  an  evaluation  of 
photomultipliers  with  a  bialkali  cathode  has  been  included.  Few 
photoelectrons  are  obtained,  but  the  input  efficiency  can  be 
good.  The  noise  factor  will  be  the  worst  because  of  the  low 
voltage  to  the  first  dynode;  100  V  usually  producing  27  second¬ 
aries  yields  a  factor  of  1.19,  thus  the  net  SNR  is  0.27. 

The  conclusion  is  that  the  preferred  designs  all 


represent  a  factor  of  two  gain  in  SNK  and  in  tilt  or  phase  ac¬ 
curacy  over  wavefront  sensors  that  must  use  PMT’s  to  provide  an 
optimum  response. 


7  0  THE  21  ELEMENT  BREADBOARD  SYSTEM 


The  purpose  of  a  breadboard  system  is  to  test  the  concepts 
and  design  of  a  large  scale  system.  The  characteristics  of 
the  AWS  that  should  be  tested  are  as  follows. 

3 

(1)  The  LSI  configuration  accuracy  as  a  function  of 
photon  count,  target  size  and  wavelength. 

(2)  The  ability  of  a  single  serial  processor  to  produce 
phase  error  signal  from  256  elements,  in  100  ps. 

In  order  to  achieve  these  goals,  it  is  necessary  to 
process  that  data  from  only  one  subaperture  if  the  subaperture 
data  are  processed  with  the  same  algorithms  that  are  applied 
to  the  full  array  and  at  the  same  speed.  A  21  element  bread¬ 
board  would  not  be  any  more  complex,  except  for  more  detectors, 
but  it  would  yield  wavefront  data  that  could  be  applied  to  a 
closed-loop  system. 

The  AWS  consists  of  three  elements:  an  optical  system, 
a  detection  system  and  the  electronic  processor.  Of  these, 
the  design  of  this  optical  system  is  not  dependent  on  the 
number  of  subapertures. 

The  nature  of  the  detection  system  is  the  major  variable 
in  the  design  of  a  breadboard.  If  the  breadboard  is  to 
represent  a  full  system  and  demonstrate  its  performance,  then 
it  should  be  quantum  noise  limited.  This  level  of  operation 
can  only  be  obtained  by  using  photomultiplier  tubes  with  noise- 
free  gain,  such  as  those  discussed  in  Section  6. 

All  of  the  detectors  which  are  recommended  for  the  full 
system  require  some  degree  of  development,  and  it  would  be 
wasteful  to  implement  any  of  these  as  a  21  element  system, 
although  the  crossover  might  be  at  approximately  64  subaper¬ 
tures.  If  a  special  intensified  self  scanned  disk  array  is  to 
be  designed  and  built,  a  256  element  chip  may  as  well  be  used. 
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Because  of  this  and  because  of  the  cost,  it  is  recommended  that 
photomultipliers  be  used  in  the  breadboard  The  cost  per  channel 
will  be  relatively  high.  At  $200  per  tube,  including  photon 
counting  circuit,  with  four  tubes  per  subaperture,  the  detector 
system  cost  would  be  $16,800  for  21  channels.  Since*  no  engineer¬ 
ing  time  would  be  required  to  build  the  tubes,  this  would 
probably  be  the  least  expensive  system 

The  major  difference  between  the  use  of  an  intensifier 
and  an  array  of  PMT's  is  size.  The  smallest  photomultipliers 
are  the  Hamamatsu  R6471  13  mm  diameter  tubes.  These  art*  13 
times  larger  than  the  lenticular  array  lens  element?;  which 
serve  to  divide  the  pupil  in  the  proposed  optical  system. 

An  enlarging  optical  system  could  be  used  instead  of  the  fiber 
optical  area-linear  transducers.  The  subaperture  size  ior  a  21 
channel  system  covering  the  same  pupil  area  as  a  200  channel 
system  is  3.2  times  larger.  The  size  of  a  detector  array  is 
only  75  mm  in  diameter  and  the  length,  including  socket, 
resistance  ladder  network  for  dynode  voltage,  preamplifier  in 
discriminator  chip  and  resettable  counter  would  be  about  150  mm 
This  is  about  the  size  of  the  intensi f iers  discussed  in  Section 
6.  A  straightforward  enlarging  optical  system  using  a  conven¬ 
tional  eyepiece  would  not  be  overly  large  and  would  be  cost 
e f feet ive . 

An  alternate  method,  providing  a  better  test  of  tin*  pupil 
dividing  optics  of  the  full  system,  would  be  to  use  a  fiber 
optic  bundle  with  lenslets  at  the  input,  and  then,  instead  ol 
bringing  the  fibers  together  into  a  line,  spread  then  out, 
taking  one  fiber  to  each  photomultiplier.  These  would  be 
optically  coupled  to  the  end  windows 

If  the  fiber  is  coupled  to  the  tube  through  a  small  ■lf>n 
prism,  then  a  majority  of  the  photon.-;  that  would  otherwise  pass 
through  the  semitransparent  photorut hode  would  would  he  inter- 


nally  reflected,  and  pass  through  the  cathode  material  many  times 
as  it  reflects  off  both  surfaces.  This  will  enhance  the  real 
response  of  the  multi-alkali  cathodes  that  can  be  obtained  in 
the  1/2  inch  tubes,  so  that  there  is  some  spectral  similarity 
to  the  GaAs  cathodes  to  bo  used  in  the  linal  system 

The  rest  of  t  tie  detection  system  would  consist  of  an 
amplifier-discriminator,  which  can  now  be  obtained  in  a  single 
chip  and  provide  a  1 00  MU/,  maximum  count  rate,  a  12  bit  counter 
and  some  kind  of  buffer  memory  for  each  tube.  These  simulate 
the  integrators,  the  oil-chip  analog  sampling  circuits,  and  the 
A/D  converter.  After  these  circuits,  the  electronics  become 
identical  to  the  AWS  circuit. 

The  AWS  circuitry,  described  in  Section  5,  consists  of 
three  major  elements:  the  Front.  Knd  Processor  <  FKP ) ,  the 
Programmable  Algorithm  Processor  (PAP),  and  tin-  System 
Supervisor  Processor  (SSP). 

The  FEP  is  largely  replaced  bj  the  parallel  circuitry  of 
the  detector  system.  The  normal  buffer  in  the  FKP  is  a  random 
access  memory  which  can  be  used  because  the  data  is  obtained 
serially  regardless  of  the  type  <»,  «i.  •  1 1  >  r  --  whether  a  self- 

scanned  diode  array  or  a  CCD.  The  pat. ».  lei  structure  of  tin  de¬ 
tectors  permits  simulation  of  tin  sampling  process  oi  the  self- 
scanned  arrays  or  the  CCD's  b>  u.-.ing  the  approximate  liming 
i'. -cause  CCD  detectors  are  gated  mu.,  transfer  lines  in  parallel, 
simulation  of  the  process  will  required  Duller  registers  winch 
would  then  be  read  out  : n  serial  fashion  direct  1\  into  another 
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can  be  smaller  because  the  amount  of  data  is  only  a  tenth  as 
great.  In  the  AW'S  PAP,  four  arithmetic  units  were  required  in 
each  processor  unit  to  handle  256  subapertures  at  a  10  kHz  field 
rate.  In  the  breadboard,  only  one  All  will  be  required.  Further 
the  instruction  set  in  each  microprogrammed  processor  can  be 
larger  so  that  the  number  of  processors  can  be  proportionally 
fewer.  The  validity  of  the  PAP  is  not  reflected  by  the  size. 

The  speed  of  the  PAP  will,  in  all  cases,  be  the  same. 

The  System  Supervisor  Processor  will  be  unchanged  in 
building  the  breadboard.  Once  the  breadboard  has  been  completed 
a  large  amount  of  the  full  system  will  have  been  completed.  The 
optics  will  be  done1.  The  detector  modules  and  fiber  optics 
interface  elements  will  simply  not'd  to  be  replaced  to  turn  the 
breadboard  into  a  256  subaperture  system. 

The  electronics  for  the  AWS  will  have'  boon  at  least  half 
completed  at  the  breadboard  stage,  and  the  critical  element 
designs  will  have  been  done*.  The?  changes  in  upgrading  to  the 
AWS  will  mainly  consist  of  upgrading  the  FKP. 

The  relationship  of  the,*  breadboard  to  the  lull  AWS 
is  summarized  in  Table  12.  The  major  element  in  building  a 
full  256  element  system  will  be  the  procurement  of  the  detectors 
Because  of  the  detector  interchangeability,  we  propose  that 
the  breadboard  be  a  little  more  than  a  model  so  that  the  optics 
need  be  done  only  oner. 


TABLE  12 


ITEM  BREADBOARD  AWS 

OPTICS  Full  System  Complete 

Detectors  4  sets  of  21  PMT ' s  4  intensified  arra 

( interchangeable  sizes) 


Electronics 

FEP  Special  for  parallel  serial  processor 

PMT ’ s  for  arrays 

PEP  ^  full  system  full  system 

SSP  full  system  complete 
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8.0  CONCLUSION 


The  objective  of  this  study  was  the  evolution  of  a  new  sen- 
3 

sor  concept,  the  LSI  sensor,  for  potential  cost  and  engineering 
improvements  over  the  CIS  system  presently  under  development. 

The  results  of  analysis  have  shown  that  the  wavefront  sensor  will 
not  only  significantly  reduce  the  cost  of  the  sensor  system  and 
provide  a  greatly  simplified  data  processing,  but  will  also  im¬ 
prove  the  performance  by  a  factor  of  two  for  low  target  signal 
levels.  The  features  of  the  system  include  the  following. 

Detectors:  The  recommended  detector  system  designed  for  the 

256  subaperture  application  uses  a  photocathode  electron  multi¬ 
plier  tube  with  direct  electron  bombardment  of  a  linear  self 
scanned  diode  array.  The  detector  configuration  employs  internal 
reflective  optics  so  that  an  opaque,  wide  bended  GaAs  photocath¬ 
ode  can  be  utilized  so  that  full  utilization  of  available  solar 
illumination  in  the  0.45  to  0.85  pm  band  is  achieved. 

Optics:  The  optical  system  as  designed  is  a  highly  com¬ 

pact,  all-reflective  system  to  take  advantage  of  the  wide  band 
detector  capability.  Two  small  galvanometer  nutation  devices  are 
the  only  moving  elements  in  the  sensor  system.  Using  a  single 
prism  divider  as  the  critical  alignment  element,  and  a  local  ref¬ 
erence  input,  auto-alignment  of  the  system  is  achievable. 
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Sensor  concept:  The  LSI  sensor  makes  use  of  the  best  fea¬ 
tures  of  both  DC  Hartmann  and  AC  shearing,  i  n!  or  femmr  or 
wavefront  measurement  systems.  The  Hartmann  technique  allows 
complete  utilization  of  photon  for  highest  plane  measurement  ac¬ 
curacy  for  low  SNU  levels,  while  the  AC  shearing  inter  IVrninHi  r 

avoids  many  of  the  object  shape  dependencies  of  the  d.c.  Hartmann. 
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In  the  LSI  ,  AC  modulation  is  obtained  by  nutation  of  the  full 
aperture  image  on  a  quad  division  prism  while  t  ht*  measurement 
uses  both  AC  and  DC  components  in  a  llartmann-1  ike  quad-cel! 
tilt  detection  algorithm.  And  optimizing  this  algorithm  allows 


3 

the  LSI  to  achieve  a  factor  of  two  improvement  over  the  shear¬ 
ing  interferometer  system. 

Processing:  The  optimized  algorithm  along  with  the  use  of 

quantum  efficient  self  scanned  detector  array  system  is  compati¬ 
ble  with  high  speed  serial  data  processing.  Thus,  instead  of  a 
processor  per  channel,  as  in  the  CIS  system,  only  4  flexible  ser¬ 
ial  processors  are  required  for  256  subaperture  measurements. 

The  measurement  update  rate  is  at  the  nutation  frequency  of 
10  kHz. 

3 

LSI  sensor  design  therefore  successfully  achieves  the 
goals  of  the  study  and  is  recommended  for  consideration  in  all 
visible  light,  large  scale  image  compensation  systems. 
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